Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 264,098 G→T L50L (CTC→CTA ykfH ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA264,0980GT92.9% 72.5 / ‑5.2 29L50L (CTC→CTAykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base T (16/10);  minor base C (1/1);  total (17/12)
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

CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_pae_tpiA/263952‑264242
                                                                                                                                                  |                                                                                                                                                
cGGGCGCCCAGGTGCAGGGTAATGTGCCGCTCCAGGCTCCTGGTAGGGTTGGTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGGTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                                <  2:66647/149‑1 (MQ=38)
       ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                                         >  1:147784/1‑149 (MQ=255)
        cAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGGGGTCCTGGCGGCGTTCGGCGTATACCGCTCGGACTCCgg                                                                                                                                        >  1:358066/1‑149 (MQ=255)
                                     cccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                           >  1:355889/1‑149 (MQ=255)
                                     cccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                           >  2:367705/1‑149 (MQ=255)
                                      ccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGGGTAAt                                                                                                          >  1:382926/1‑149 (MQ=255)
                                             ggggTGCTCATTATGGTGCGCCACGATTGTCAGTGCAGGGTGATGCTCAGCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGAc                                                                                                   <  2:69637/149‑1 (MQ=21)
                                             gggTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGAc                                                                                                   <  1:291102/149‑1 (MQ=255)
                                              ggTTGCTCATTACAGTTCTCCACTATTGTCAGTTAAGCGCGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTCCCGGTGTAATGACCGACg                                                                                                  >  2:87025/1‑149 (MQ=255)
                                                 ttcTCATTATAGTTGTCCAGGAGTGGCATTTCAGGGTGTTGCTCATCAGGCAGGTATAGGGCCCAGTGGGGTCCTGGCGGCGTTCGGGGTATACCGCTAGGACTCCGGCGATATCCGGAAAGTCCCTGCCTGTCTAATGACAGACGCTa                                                                                               <  2:34155/147‑1 (MQ=255)
                                                      aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTg                                                                                          >  1:386171/1‑149 (MQ=255)
                                                               ctcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTa                                                                                 <  2:386171/149‑1 (MQ=255)
                                                                  cggtgatGTCAGTTCCGCGTGATTCTCATCAGGCAGGTAAAGGGCTCTTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGCCTCCGGCGATATGCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAttt                                                                              <  2:476624/143‑1 (MQ=18)
                                                                     tATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGcc                                                                           <  1:449864/149‑1 (MQ=255)
                                                                     tATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCAATGGTATTTGcc                                                                           >  1:67649/1‑149 (MQ=255)
                                                                         gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTg                                                                       >  1:386353/1‑149 (MQ=255)
                                                                                 aGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCg                                                               >  1:473957/1‑149 (MQ=255)
                                                                                   ggTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                           <  2:257023/103‑1 (MQ=255)
                                                                                   ggTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                           >  1:257023/1‑103 (MQ=255)
                                                                                     tGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAg                                                           <  1:405626/149‑1 (MQ=255)
                                                                                     tGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGGTCGGCGTATACCCCCACCACTCTGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAg                                                           >  1:474937/1‑149 (MQ=21)
                                                                                                 ggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                               <  1:367705/149‑1 (MQ=255)
                                                                                                           gggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGc                                     >  1:197242/1‑149 (MQ=255)
                                                                                                                 aTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTCGCGGAAGATTCGAGACTCAGGATGCTGGCGGcatt                               >  2:62130/1‑148 (MQ=17)
                                                                                                                                  gTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGt              >  1:485516/1‑149 (MQ=255)
                                                                                                                                         gTATACCGCCAGGCCTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAGGACAGACGCTACCGTGTCACTGGTATTTGCCGGTGCAGTTGCGAAAGATTCGGGACTCCGCCTGTACCTGGTAGGGCGTCATTGCCCTGCGTTTACTGa       <  1:249087/149‑1 (MQ=255)
                                                                                                                                          tataCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGAt      >  2:355052/1‑149 (MQ=255)
                                                                                                                                              ccGGGAGGGCGCCGGCGATATCCGGAACGTCCCTGCCGGGGTAAGGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGCAGCCAAAGTCTCTGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAtt  <  1:396686/149‑1 (MQ=14)
                                                                                                                                              ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCGGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCGCTGCGTTTGCTGGTAAtt  >  2:163289/1‑149 (MQ=17)
                                                                                                                                                  |                                                                                                                                                
CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_pae_tpiA/263952‑264242

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