Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_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_pae_tpiA2,795,3960CA97.4% 114.5 / NA 38Q382K (CAG→AAG) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (15/22);  minor base T (1/0);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_pae_tpiA/2795251‑2795539
                                                                                                                                                 |                                                                                                                                               
cccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGt                                                                                                                                              >  2:632713/1‑149 (MQ=255)
cccGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGt                                                                                                                                              >  1:192681/1‑149 (MQ=255)
   gTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGc                                                                                                                                           <  1:472504/149‑1 (MQ=255)
         ggCGGTGTTGCTCTCCGCCGGCCCGGCATTCTTGACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTTTTTGAATTTCTGATTGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTTGTCCTCGCCGTTTCTAAGTTGCGGGTgc                                                                                                                                     <  2:221034/149‑1 (MQ=255)
            ggTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCCAAGTTTCGGATGCGCa                                                                                                                                  >  2:407501/1‑149 (MQ=255)
                 tACTCTCCACCGGCGCGGCATTCTTAGGGGTGGGGTTTAACTATTACGCGCCTGCGCAGGTGTTTAAATTTCTGATCGACAGATCCGGCGCTATCGCCCTAGTGGTTTATTTAGTCATCGCCGTTTCAAAAGTGCGGGTGCGCAAAAtt                                                                                                                             >  2:669019/1‑149 (MQ=255)
                  aCTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTc                                                                                                                            <  1:175636/149‑1 (MQ=255)
                                                ggtggtgCACGATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGc                                                                                              >  2:7166/1‑149 (MQ=255)
                                                   ggtgAACTATTACGCGCCTGCGAAGGTATTTCAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGaa                                                                                           <  1:439824/149‑1 (MQ=255)
                                                    gtgAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAAt                                                                                          <  2:146175/149‑1 (MQ=255)
                                                        aCTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTACCGACCTGCTGGTTTATTTCTTCATCGCGGGTTCTAAGTTGCGGATGCGCAAACTTCTGCGACCAGAAGGGAGCGAAATTGGCTTCCACATGTgg                                                                                      >  2:424631/1‑149 (MQ=255)
                                                        aCTATTACGCGCCTGCGAAGGTATTTAAAGTTCTGCTCGTCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTgg                                                                                      <  1:303515/149‑1 (MQ=255)
                                                        aCTATTACGCGCCTGCGAAGGGATTTAAATTTCTGAGCGGAAGCTCCGGCGCTATCGCCCTGATGATTTAAATAATAATCGCCGTTTCGAGGTTGCGGATGCGCAAAATTTTGCGAGCAGAAGGCAGCGAAATTCGCTTGCGACTGTgg                                                                                      >  2:424596/1‑149 (MQ=255)
                                                                    cTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTgg                                                                          <  2:401701/149‑1 (MQ=255)
                                                                     tGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGATTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGc                                                                         >  1:456334/1‑149 (MQ=255)
                                                                               tttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCt                                                               >  1:190447/1‑149 (MQ=255)
                                                                               tttAAATTTCAGCTCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGCCTACCCTGGcc                                                               >  2:257850/1‑148 (MQ=255)
                                                                                ttAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTa                                                              <  2:362663/149‑1 (MQ=255)
                                                                                 tAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAg                                                             <  2:622122/149‑1 (MQ=255)
                                                                                  aaaTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGt                                                            <  1:356843/149‑1 (MQ=255)
                                                                                    aTTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCa                                                          <  2:428511/149‑1 (MQ=255)
                                                                                        cTGATCGACAGCTCCGACGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTCGCTAACATAGCTAGTCATccg                                                      >  2:530588/1‑147 (MQ=255)
                                                                                             cGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTtta                                                 <  1:7166/149‑1 (MQ=255)
                                                                                               aCAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAGGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGCCTCGTCACCGGCTttatt                                               >  2:430055/1‑149 (MQ=255)
                                                                                                cAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTttatta                                              <  1:632713/149‑1 (MQ=255)
                                                                                                       ggCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCCCAACTGGCTAGTCATCGGCTTTATTACCtttct                                       >  2:206403/1‑147 (MQ=255)
                                                                                                       ggCGCTATCGACCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt                                       <  1:688109/149‑1 (MQ=255)
                                                                                                           cTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTg                                   <  1:163163/149‑1 (MQ=255)
                                                                                                                       ttggttatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTATTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTa                       <  1:275034/145‑1 (MQ=255)
                                                                                                                         gataatataGTCAGCGCGGTTTCATTGTTGCGGATTGTCAAAATTCTGGGAGCAGAAGGAGGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTAGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTCGTGATGCTAtt                     >  2:584502/8‑149 (MQ=255)
                                                                                                                         gTTTATTTAGTCATCGCCGTTTCAAAGGTGCTGATTGGAAATATTCTGCGAGCAGAAGGAAGTCAATTTCGCTTTCGAATGTTGTTTGATCCGTGGCTCACCTGGCTAGTCATCGTCTTTATTACCTTTGTGTTGGTAGTGATGCTAtt                     <  2:321287/149‑1 (MQ=255)
                                                                                                                           ttatttaGTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATtcc                   <  1:163145/149‑1 (MQ=255)
                                                                                                                                  gTCATCGCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTAGTCCGTCCGgc            >  1:429562/1‑149 (MQ=255)
                                                                                                                                        gCCGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACa      <  1:404325/149‑1 (MQ=255)
                                                                                                                                         ccGTGTCAAAGTTGCGCATGAGCAAAGTTCTGCGAGCCGAGGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCATCCGGCGCAACAg     <  2:630458/149‑1 (MQ=255)
                                                                                                                                          cGTTTCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGGGGCTCACCTGGCTAGTCACCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGt    <  1:597290/149‑1 (MQ=255)
                                                                                                                                           gcttCAAAGTTGCCGATGCCCAAAATTCTGCGCGCAGAAGGAATCTAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCCCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGtt   <  1:10296/147‑1 (MQ=255)
                                                                                                                                            tttCAAAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTTTTGGTAGTGATGCTATTCCGTCCGGCTCAACAGTTa  >  1:55552/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCCGTACGTGGCGGTGTTACTCTCCACCGGCGCGGCATTCTTAACGGTGGTGGTGAACTATTACGCGCCTGCGAAGGTATTTAAATTTCTGATCGACAGCTCCGGCGCTATCGCCCTGCTGGTTTATTTAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGTGATGCTATTCCGTCCGGCGCAACAGTTA  >  NC_000913_3_pae_tpiA/2795251‑2795539

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: ‑