Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 701,079 T→C 7.7% L148P (CTC→CCC)  U712_03405 → putative RNA methyltransferase yefA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881701,0790TC 7.7% 144.8 / 4.1 53L148P (CTC→CCC) U712_03405putative RNA methyltransferase yefA
Reads supporting (aligned to +/- strand):  ref base T (30/18);  new base C (2/2);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

AACAAGTCAAAGACGTCTTAGAACGCATCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACGGTGTCAAAGCCTACAACG  >  CP006881/700941‑701216
                                                                                                                                          |                                                                                                                                         
aaCAAGTCAAAGACGTCTTAGAACGCATCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACt                                                                                                                                           >  5:360359/1‑139 (MQ=255)
aaCAAGTCAAAGACGTCTTAGAACGCATCAGCAAGCTGGACCTTTCAAAAGTCACCGTGCCCCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCCCAGGTGCCGGTAGGAGAACGGGAAGGCGGACt                                                                                                                                           >  4:269458/1‑139 (MQ=255)
     gTCAAAGACGTCTTAGAACGCATCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACtcgtcg                                                                                                                                      <  2:317871/139‑1 (MQ=255)
         aaGACGTCTTAGAACGCATCGGCAAGCTGGACCTTTCAAAAGTCACCTTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCggg                                                                                                                                  >  7:135501/1‑139 (MQ=255)
         aaGACGTCTTAGAACGCATCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCggg                                                                                                                                  >  7:335574/1‑139 (MQ=255)
         aaGACGTCTTAGAACGCATCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCgg                                                                                                                                   >  6:236847/1‑138 (MQ=255)
         aaGACGTCTTAGAACGCAGCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCggg                                                                                                                                  <  7:129415/139‑1 (MQ=255)
                   aGAACGCATCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTAt                                                                                                                           <  6:360359/136‑1 (MQ=255)
                        gCATCTGCAAGCTGGAACTTTCAATAGTCACGGTGCACCCGCCACTCGGCATGGTAGACCGGTGGAACTACTGAAACATCGGAGATGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCTaag                                                                                                                   <  8:135501/139‑1 (MQ=255)
                        gCATCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAaag                                                                                                                   <  8:335574/139‑1 (MQ=255)
                         cATCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCGATCAGCAaag                                                                                                                   >  3:232814/1‑138 (MQ=255)
                          aTCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAACCAAAGCACAGGTCCCGGTAGGAGAAGGGGAAGGCGGACTCGTCGCGGGCTTCTAACAGGAaagta                                                                                                                 >  6:6654/1‑137 (MQ=255)
                          aTCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAaag                                                                                                                   >  2:2123/1‑137 (MQ=255)
                          aTCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACCCGTCGCGGGATTCTATCAGCAaagaa                                                                                                                 >  2:273539/1‑139 (MQ=255)
                          aTCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACCCGTCGCGGGATTCGATCAGCAaagaa                                                                                                                 >  2:268601/1‑139 (MQ=255)
                             ggCTAGCGGGACCTTTCAACAGGCCCCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGGCGTGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGcc                                                                                                              <  1:85591/139‑1 (MQ=255)
                                 aGCTGGACATATCAAAAGACACCGTCCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGGCGCGGGATTCTATCAGCAAAGAAGCCATGa                                                                                                          <  3:49768/139‑1 (MQ=255)
                                    tGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGAcat                                                                                                       >  1:4291/1‑139 (MQ=255)
                                     ggACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGAcatc                                                                                                      <  5:236847/139‑1 (MQ=255)
                                       aCCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGAcatcat                                                                                                    >  8:14871/1‑139 (MQ=255)
                                       aCCTTACAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGAcatcat                                                                                                    <  8:12084/139‑1 (MQ=255)
                                                     aCCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAaag                                                                                                                   >  6:236944/1‑110 (MQ=255)
                                                     aCCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAaag                                                                                                                   <  5:236944/110‑1 (MQ=255)
                                                                 acacTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGc                                                                          <  2:4291/139‑1 (MQ=255)
                                                                        gCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGCGCGCCTGCCTGCTCCAGCAATCCa                                                                    >  3:26933/1‑138 (MQ=255)
                                                                        gCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCa                                                                    >  3:22022/1‑138 (MQ=255)
                                                                                 aCCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACCTCATCCACATGAGCGCCTGCCTGATCCAGCAATCCAAAAacgacc                                                           >  8:371544/1‑137 (MQ=255)
                                                                                   ccGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAg                                                        >  7:263479/1‑139 (MQ=255)
                                                                                      tGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGGTTCAATCAGCAAAGAAGCCATGACATCAACCACATGAGCGCCAGCCAGATCCAGCCAACCCAAAACGACCAAAcc                                                      >  6:52337/1‑138 (MQ=255)
                                                                                        gAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGCCGAAGCCGt                                                    >  3:115901/1‑138 (MQ=255)
                                                                                             aCAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAGc                                              >  3:351799/1‑139 (MQ=255)
                                                                                             aCAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAGc                                              >  1:22521/1‑139 (MQ=255)
                                                                                             aCAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCACCAATCCAAAAACGACCAACCCGTTCaatc                                              >  6:11427/1‑137 (MQ=255)
                                                                                                   aCAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAg                                               >  2:224214/1‑132 (MQ=255)
                                                                                                   aCAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAg                                               <  1:224214/132‑1 (MQ=255)
                                                                                                      aaGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGCTTCTAAAAGCAAAGAAGCCAAGACATCATGGACATTAGCACCAGCCGGTACCAGCAATCCCAAAAAGACGCAGCCCGTCAAGCTTTaaaatc                                     >  4:15811/1‑137 (MQ=255)
                                                                                                      aaGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCAGCCGGCTCCAGGAATCCAAAAACGCCGAAGCCGTTCAAGCTGTCAAAGa                                     >  8:163809/1‑139 (MQ=255)
                                                                                                      aaGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATAAGCAAAGAAGCCCTGACATCATCGACATGAGCGCCAGCCTGATCCAGCAATCCAAAAACGACGTAGCCGTTCAAGCAGTCAAAGa                                     >  7:137963/1‑139 (MQ=255)
                                                                                                        gCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGAcatcat                                                                                                    <  5:314487/74‑1 (MQ=255)
                                                                                                        gCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGAcatcat                                                                                                    >  6:314487/1‑74 (MQ=255)
                                                                                                           caGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCGAGCTGTCAAAGACAttt                                >  8:68247/1‑139 (MQ=255)
                                                                                                                     tAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTa                      <  7:14871/139‑1 (MQ=255)
                                                                                                                       ggAGAACGGGAAGGCGGCCCCGCCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACg                    <  1:268601/139‑1 (MQ=255)
                                                                                                                       ggAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACg                    <  1:2123/139‑1 (MQ=255)
                                                                                                                       ggAGAACGGGAAGGCGGACCCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACg                    <  1:273539/139‑1 (MQ=255)
                                                                                                                       ggAGAACGGCAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCAGGACATCATCGACATGCGCGCCGGCCTGATCCAGCAAACCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACg                    <  3:269458/139‑1 (MQ=255)
                                                                                                                           aaCGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCCGCAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACGgtgt                >  1:227678/1‑139 (MQ=255)
                                                                                                                             cGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACGGTGTc               >  1:169006/1‑138 (MQ=255)
                                                                                                                             cGGGAAGGCGGACGCTGCGCGGGATTCTAGCAGCAAACAACCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGAAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTATGGTGTCa              <  4:351799/139‑1 (MQ=255)
                                                                                                                                gAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACGGTGTCAAAg           >  6:340990/1‑139 (MQ=255)
                                                                                                                                    gCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAACACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACGGTGTCAAAGCCTa       <  4:22022/139‑1 (MQ=255)
                                                                                                                                    gCGGACTCGGCGCGGGATGCTATCAGCAAAGCAGCCATGACATCATCGACATGAGCGCCTGCCTGAGCCAGCAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACGGTGTCAAAGCCTa       <  4:26933/139‑1 (MQ=255)
                                                                                                                                     cGGACTCGTCGCGGGATTCTTTAAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCGGATCCAGCAATCCAAAAACGACGCAGCCGTTCAAGCGGTCAAAGACATTTGCGCAAACTACGGTGTAAAAGccaa       >  6:241841/1‑136 (MQ=255)
                                                                                                                                     cGGACTCGTCGCGGGATTCTATCAGCAAAGCAGCCATGCCCTCAGCGCCATGCGCGCCTGCCGGCTCCCGCCATCCAAAAACGCCGACGCCGTTCACGCTGGCCAAGCCCTGTGCGCACCCTCCGGTGGCCAAGccgcc      >  2:62893/1‑136 (MQ=255)
                                                                                                                                     cGGACTCGTCGCGGGATTATAACAGCAAAGAACCCATGACAAAAACCACAAGAGCCCCAGCCAGATCCCCCAAACCAAAAACCCCCCACCCGGTCAAGCTGGCAACTGCAATTGAGCAAACTACGGTGTCAAATCCAAc      >  4:280848/1‑139 (MQ=255)
                                                                                                                                         cTCGGCGCGGGATTCTATCAGAAAATACGCCATGCAAACATCCACATGAGCGCCGGCCTGACCCAGCAATCAAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACGGTGTCAAAGCCTACAACg  <  2:22521/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
AACAAGTCAAAGACGTCTTAGAACGCATCGGCAAGCTGGACCTTTCAAAAGTCACCGTGCACCCGACACTCGGCATGGAAGACCCGTGGAACTACAGAAACAAAGCACAGGTGCCGGTAGGAGAACGGGAAGGCGGACTCGTCGCGGGATTCTATCAGCAAAGAAGCCATGACATCATCGACATGAGCGCCTGCCTGATCCAGCAATCCAAAAACGACGAAGCCGTTCAAGCTGTCAAAGACATTTGCGCAAACTACGGTGTCAAAGCCTACAACG  >  CP006881/700941‑701216

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑