Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 2,795,396 C→A Q382K (CAG→AAG)  gabP → gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA2,795,3960CA96.7% 83.5 / NA 31Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (14/15);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGAT  >  NC_000913_3_vch_tpiA/2795251‑2795514
                                                                                                                                                 |                                                                                                                      
cccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTATTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGt                                                                                                                     >  1:18885/1‑149 (MQ=255)
    gaCGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGGGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTTCGGTGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCg                                                                                                                 <  1:277369/148‑1 (MQ=255)
         ggCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATgc                                                                                                            <  2:427535/149‑1 (MQ=255)
                ttACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCTGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAAt                                                                                                     >  1:273433/1‑149 (MQ=255)
                ttACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTGGCGGATGCCCAAAAt                                                                                                     >  2:31254/1‑149 (MQ=255)
                   ccctcCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCTTGCTGGTTGATTGAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCt                                                                                                  <  1:174637/147‑1 (MQ=255)
                          ccGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGCATCCGCaacttt                                                                                                    >  2:506780/1‑136 (MQ=255)
                          ccGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGCATCCGCaacttt                                                                                                    <  1:506780/140‑5 (MQ=255)
                             gcgcGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAgaa                                                                                        <  1:177402/149‑1 (MQ=255)
                                  gCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCCAGCAgaagtatg                                                                                   <  2:326381/149‑5 (MQ=255)
                                       cTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGaaa                                                                              >  2:140308/1‑149 (MQ=255)
                                        ttAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCCACAACTCCAGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAAt                                                                             >  1:244651/1‑149 (MQ=255)
                                           aCGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATAAACAGCACCGACGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCg                                                                          >  1:138568/1‑149 (MQ=255)
                                             ggtggtgctGAACTATTACGCGCCCGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCt                                                                        <  2:158781/149‑1 (MQ=255)
                                                    gtgAACTATTACGCGCCGGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAATGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAt                                                                 <  1:147516/149‑1 (MQ=255)
                                                          tATTACGCGCCTGCGGAGGTATGTAAATATCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTGTGTAGTCATCGCCGTTTCAAAGTTGAGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGACATTCGCTTGCGAACGTggtt                                                           >  2:452063/1‑147 (MQ=255)
                                                          tATTACGCGCCTGCGAAGGTATATAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCCAAGTTGCGGATGAGCAAAATTCTGCGGAAAGAAGGATGCCACGTTCGCTTGCAAATTTGGCt                                                           >  2:245990/1‑149 (MQ=255)
                                                                  gcCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGt                                                   <  1:140308/149‑1 (MQ=255)
                                                                    cTGCGAAGGTATTTAAATTTCTGATCCACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTACTCATCTCCGTTTCAAAGTTGCGGATTCTTAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                 <  2:323506/149‑1 (MQ=255)
                                                                             tATTTAAAGTTCTGCGCGACAGCTCCGGAGCTATCGTCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTTTACGAACAGAAGGAAGCGAAATTCGCTTGCGAAGGTGGCTTTATCCGTGGCTCCCCTgg                                        <  2:211687/149‑1 (MQ=255)
                                                                               tttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGCATGTGGCTTTATCCGTGGCTCACCTGGCt                                      >  1:494857/1‑149 (MQ=255)
                                                                               tttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAATCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCt                                      >  1:337541/1‑149 (MQ=255)
                                                                                ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                     >  2:197580/1‑149 (MQ=255)
                                                                                  aaaTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGt                                   >  2:158907/1‑149 (MQ=255)
                                                                                     tttCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTAAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCAt                                >  1:99676/1‑149 (MQ=255)
                                                                                                  gCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTAcc                   <  2:298406/149‑1 (MQ=255)
                                                                                                   cTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCt                  >  1:129703/1‑149 (MQ=255)
                                                                                                     ccGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCttt                <  1:158907/149‑1 (MQ=255)
                                                                                                       gtcgcTTTCGCCCTTCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCATAATTCTGCGAGCAGTAGGAAGCGAAAGTCGCTTGCGAATGTGGCGTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt              <  1:245990/147‑1 (MQ=255)
                                                                                                          gcTATCGCTCTGCTGGTTTATTTAGTCAGTGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGtt           <  2:241993/149‑1 (MQ=255)
                                                                                                                   ctgctgGTTTATTTAGTCATCGCCGTTTCATAGTTGTGGATGCACAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGACTTTATCCGTGTCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTCGtcatcat  >  2:508956/1‑143 (MQ=255)
                                                                                                                                                 |                                                                                                                      
CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGAT  >  NC_000913_3_vch_tpiA/2795251‑2795514

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