Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 514,757 T→C 15.7% S35G (AGC→GGC)  ybbJ ← inner membrane protein; stimulator of the QmcA suppressor of ftsH‑htpX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913514,7570TC15.7% 95.4 / 2.1 36S35G (AGC→GGC) ybbJinner membrane protein; stimulator of the QmcA suppressor of ftsH‑htpX
Reads supporting (aligned to +/- strand):  ref base T (27/3);  new base C (2/4);  total (29/7)
Fisher's exact test for biased strand distribution p-value = 7.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

AAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGA  >  NC_000913/514557‑514977
                                                                                                                                                                                                           |                                                                                                                                                                                                                            
aaaaCGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCAACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAACCACCGCTGCCACGCCGACCCACAACAAACAACCATTTCCGCCCa                                                                                                                                                                                                    >  1:1435455/1‑227 (MQ=255)
                       gcCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGCAATCACCGCTGCCCCGCCGCTCCACAACCAACAACCATTTCCGCCCAGCATCTCGCATGCCcgcagc                                                                                                                                                                                >  1:1619944/1‑224 (MQ=255)
                                           cTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAACTAACCATTTCCGCCCAGCATCTCGGCCGCCACCAACAAACCGCCGCGACTGAGCCAGaaa                                                                                                                                                        >  1:453458/1‑228 (MQ=255)
                                           cTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAACTAACCATTTCCGCCCAGCATCCCGGCTGCCAGCAGCCAACCGCCCAGACCGATCCa                                                                                                                                                            >  1:1774775/1‑224 (MQ=255)
                                           cTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGaaa                                                                                                                                                        >  1:1046801/1‑228 (MQ=255)
                                                    tgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCCGAAAATATGGGGa                                                                                                                                               >  1:1048244/1‑228 (MQ=255)
                                                    tgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGa                                                                                                                                               >  1:1547181/1‑228 (MQ=255)
                                                         ttGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGGTGaa                                                                                                                                           >  1:1182232/1‑227 (MQ=255)
                                                                cgcACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCCCCACAACAAATAACCATTTCCGc                                                                                                                                                                                                       >  1:1295599/1‑160 (MQ=255)
                                                                                        caccacAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAAt                                                                                                                                                                                                                   >  1:1940862/1‑124 (MQ=255)
                                                                                        caccacAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTc                                                                                                               >  1:1202226/1‑224 (MQ=255)
                                                                                         accacAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATa                                                                                                                                                                                                                  <  2:1937421/124‑1 (MQ=255)
                                                                                                        cgAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCCGAAAATATGTGGATGAACGACCATTAACCCCATCCTGGCTGAGTCCGTTTGTTGGCGCTGTCtt                                                                                            >  1:231970/1‑227 (MQ=255)
                                                                                                          agcagcGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCacaaca                                                                                                                                                                                                                      >  2:409907/1‑103 (MQ=255)
                                                                                                           gcagcGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCAcaacaa                                                                                                                                                                                                                     <  1:410589/103‑1 (MQ=255)
                                                                                                            gcgcGTCAGGAAGGCAAACAATACCCCTTGCCCCCCCCAACCCGGCCGGCCAGCGCCAACCACCCGCGCCGGAAACCCCCGCTCCCCCGCCGACCCAAAACAAAGACCCATTCCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCAc                                                                                        <  2:925388/226‑1 (MQ=255)
                                                                                                                      aTGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGGCCGTTTGTTGGTGCTGTCTTTTACCAGCTCGGCaa                                                                             >  2:1045011/1‑228 (MQ=255)
                                                                                                                                   ccccTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGgcg                                                                                                   >  1:56209/1‑193 (MQ=255)
                                                                                                                                        tGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGa                                                                                                                                            >  1:106453/1‑147 (MQ=255)
                                                                                                                                         cccccccccACCCCAGCGCCCCCCCCCAACCCCCCAGCCCGGTAACCCCCGCTCCCACGCCGCCCCACAACAAATAACCATTCCCCCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGaa                                                                                                                                           <  2:106263/141‑1 (MQ=255)
                                                                                                                                              tCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGAc                                                               >  1:196379/1‑218 (MQ=255)
                                                                                                                                                    cccAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTgg                                               >  2:1128373/1‑228 (MQ=255)
                                                                                                                                                       ggcgcCACCACCAAGCCCCACAGGCAGGTACCCACCGCTCCCCCGCCGCCCCAAAACAAACAACCATGCCCCCCAAGCACCCCGGCGGCAAGCAGCAACCCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCct                                            <  2:1086697/227‑1 (MQ=255)
                                                                                                                                                       aGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGCCTGCTGGcgt                                            >  1:289455/1‑226 (MQ=255)
                                                                                                                                                        gCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGcc                                             >  1:1088644/1‑226 (MQ=255)
                                                                                                                                                           cgaCCAGCCCTACCCCCAGGGCCGATACCCCCCGCGGCCAGCCCGTCCCAAAAAAAATACCCATTCCCCCCCAGACTCCGGGCTGCCAGCAACAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTa                                         <  2:231576/226‑1 (MQ=255)
                                                                                                                                                                                   cGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCt                   >  2:735565/1‑228 (MQ=255)
                                                                                                                                                                                   cGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACACCATCAGGCTGCTGGCCTCTCATGGCCTCAGCACTACCTTGCt                   >  2:1876797/1‑228 (MQ=255)
                                                                                                                                                                                            ccGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCa                                                                                      >  2:1951744/1‑152 (MQ=255)
                                                                                                                                                                                             cGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAg                                                                                     <  1:1955226/152‑1 (MQ=255)
                                                                                                                                                                                               cTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTgga          >  2:1564799/1‑225 (MQ=255)
                                                                                                                                                                                                tGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGGCTCTAATGGCATCATCACTACTTTGCTGTTACTgga          >  2:1334420/1‑224 (MQ=255)
                                                                                                                                                                                                tGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTggt          >  2:682229/1‑222 (MQ=255)
                                                                                                                                                                                                    aCGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGa  >  2:793846/1‑228 (MQ=255)
                                                                                                                                                                                                          cTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATcc                                                                          >  1:921234/1‑150 (MQ=255)
                                                                                                                                                                                                          cTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGgctgc                                                  >  1:1513699/1‑174 (MQ=255)
                                                                                                                                                                                                           tCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATccc                                                                         <  2:919617/150‑1 (MQ=255)
                                                                                                                                                                                                           tCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGgctgct                                                 <  2:1510990/174‑1 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                            
AAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGA  >  NC_000913/514557‑514977

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