Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,444,462 G→T 15.8% L157L (CTC→CTA) ‡ sulA ← cell division inhibitor SulA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,444,4620GT15.8% 89.5 / 9.6 38L157L (CTC→CTA) ‡sulAcell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base G (15/17);  new base T (3/3);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

CATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTT  >  NC_002947/2444321‑2444609
                                                                                                                                              |                                                                                                                                                   
cATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGAt                                                                                                                                             >  2:833583/1‑150 (MQ=255)
cATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCGGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGAGCGGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGAt                                                                                                                                             <  2:1122126/150‑1 (MQ=255)
      cAGTGTGGGCGTGCAGATGGCCGGTCGCGTGATGGATGGCGAGCTCGGGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGTAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGTCGAGGCGGATGTTCAg                                                                                                                                       <  2:55315/150‑1 (MQ=255)
           tgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCttt                                                                                                                                  >  1:469222/1‑150 (MQ=255)
             ggCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACCTCCGCCGAGGCGGATGTTCAGGCTTtgt                                                                                                                                <  1:996783/150‑1 (MQ=255)
             ggCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTtgt                                                                                                                                <  2:418498/150‑1 (MQ=255)
                                   tGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCAc                                                                                                          >  2:809028/1‑150 (MQ=255)
                                        gatggCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGc                                                                                                     <  1:487761/150‑1 (MQ=255)
                                          tggCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAg                                                                                                   >  2:767512/1‑150 (MQ=255)
                                                 cTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCt                                                                                                              >  1:651510/1‑132 (MQ=255)
                                                 cTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCt                                                                                                              <  2:651510/132‑1 (MQ=255)
                                                  tCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGc                                                                                           <  1:809028/150‑1 (MQ=255)
                                                    gAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCa                                                                                                                                        >  2:652288/1‑103 (MQ=38)
                                                    gAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCa                                                                                                                                        <  1:652288/103‑1 (MQ=38)
                                                        gTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAg                                                                                     >  1:1002740/1‑150 (MQ=255)
                                                          gatgaAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCg                                                                                   <  2:1118738/150‑1 (MQ=255)
                                                               aGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCAtt                                                                              <  2:823717/150‑1 (MQ=255)
                                                                 ggCCGCGTGTGCTGCACTGAGTGTGCAGCAGCTTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGa                                                                            <  1:964951/150‑1 (MQ=255)
                                                                 ggCCGCGTGTGCTGCACGGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCGTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGa                                                                            <  2:332283/150‑1 (MQ=255)
                                                                          tgctgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCcca                                                                   <  2:1002740/150‑1 (MQ=255)
                                                                                aCTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCa                                                                                                                                        <  2:893844/75‑1 (MQ=21)
                                                                                aCTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCa                                                                                                                                        >  1:893844/1‑75 (MQ=21)
                                                                                    aGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTcac                                                         >  2:990494/1‑150 (MQ=255)
                                                                                          cagcagCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGgt                                                   >  2:311682/1‑150 (MQ=255)
                                                                                               gcGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGATCATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGc                                              >  2:649096/1‑150 (MQ=255)
                                                                                                   tATCAGGCTATTTGAAGTGACAAGCAAAGCCTGAACATCAGCCGAGCCGCATGTTCAGGCTTTGTGCTTTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGa                                           <  1:860479/150‑1 (MQ=255)
                                                                                                        ggCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGccc                                      >  2:6639/1‑150 (MQ=255)
                                                                                                         gCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCa                                     >  2:367521/1‑150 (MQ=255)
                                                                                                          cTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAg                                    <  1:502177/150‑1 (MQ=255)
                                                                                                                      aCAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAg                                          >  2:564798/1‑132 (MQ=255)
                                                                                                                      aCAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAg                                          <  1:564798/132‑1 (MQ=255)
                                                                                                                                ccTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGt              >  2:370640/1‑150 (MQ=255)
                                                                                                                                 cTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGtt             >  2:979384/1‑150 (MQ=255)
                                                                                                                                  tGAACATCCGCCTCGGCTGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCATGCAGCGGCATTGACGTTGCcca                                                                   <  2:249299/95‑1 (MQ=25)
                                                                                                                                  tGAACATCCGCCTCGGCTGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCATGCAGCGGCATTGACGTTGCcca                                                                   >  1:249299/1‑95 (MQ=37)
                                                                                                                                      cATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAgc        <  2:175959/150‑1 (MQ=255)
                                                                                                                                            ccGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCtt  >  2:13525/1‑150 (MQ=255)
                                                                                                                                            ccGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCtt  <  1:6639/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCGCTT  >  NC_002947/2444321‑2444609

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