Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 1,196,220 C→T H366H (CAC→CAT icd → isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA1,196,2200CT88.2% 27.6 / ‑5.5 17H366H (CAC→CATicdisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/7);  minor base A (1/1);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGG  >  NC_000913_3_vch_tpiA/1196074‑1196314
                                                                                                                                                  |                                                                                              
gggcgtatcggtatcgCCCCTGGTGCAAACATCGGTGTCGAATGCGCCCTGTTTGAAGCCACTCACGGTACTGCGCCGAAAGATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCatat                                                                                              <  1:44784/145‑1 (MQ=255)
   ggtatcggtatcGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCAATTATTCTCTCCGCTGAGATGATGCTGGGCCATATggg                                                                                           >  1:340244/1‑149 (MQ=255)
       tcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTATACGATGAGAATGTGCTGCGACAAATGGGTTgg                                                                                       >  1:143441/1‑149 (MQ=255)
               gCCCCTGGTGCAAACATCGGTGACGCATACCAACGCGTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACCAAGTAAATCCTGGCTCAATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGACCTGAAGc                                                                               >  2:425103/1‑149 (MQ=255)
                       tGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcc                                                                       >  1:217537/1‑148 (MQ=255)
                            aCATCGGTGACGAATGGGCCCTGTTTGAAGGCACGCACGGTACTGCGCCGAAATATGCAGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAtt                                                                  <  2:395268/149‑1 (MQ=38)
                              aTCAGTAACGAATAACCCCTGTTTGAAGCCACCCACGGTACTGCGCCAAAATAAGCCGCTCAGGACAAAGTAAATCCAGGCTCTATTATTCCCTCCGCTGAGATGATGCTGCGCCATATGGGTTAGACTGAAGCGGCTGACCTGATTGt                                                                <  2:10302/149‑1 (MQ=25)
                                 ggTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCACCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaa                                                             >  2:34650/1‑149 (MQ=38)
                                  gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaaa                                                            >  1:374739/1‑149 (MQ=255)
                                  gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaaa                                                            <  2:217537/149‑1 (MQ=255)
                                               ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                               >  1:369807/1‑149 (MQ=255)
                                                       aaGCCAACCACGGTAATGCGCCGAAATATGCCGGCCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaat                                       <  2:196012/149‑2 (MQ=25)
                                                         gCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATTATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTGAAAGGTATGGAAGGCGCAATCAAtgc                                     <  1:34650/149‑1 (MQ=25)
                                                            acccacGGTACTGCGCCGAAATATGCCGGTCACGACAAAGTTAATCCTGGCTGTATTATGCGCTCCGTTGAGATGATGGTGCGCCAAATGGTTTGGACTGAAGCGGCTGACCTGATTCTTAAAGGTATGGAAGGCGACATCAATGCCaa                                  <  2:297459/149‑1 (MQ=18)
                                                                                     ccGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGa                                                    >  1:434246/1‑106 (MQ=21)
                                                                                     ccGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGa                                                    <  2:434246/106‑1 (MQ=21)
                                                                                         tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTga     <  2:374739/149‑1 (MQ=14)
                                                                                            ggACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAGCGTCTgatgg  >  1:117475/1‑149 (MQ=14)
                                                                                                                                                  |                                                                                              
GGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGG  >  NC_000913_3_vch_tpiA/1196074‑1196314

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: