Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,444,458 A→C 18.4% *159G (TGA→GGA)  sulA ← cell division inhibitor SulA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,444,4580AC18.4% 73.4 / 12.2 37*159G (TGA→GGA) sulAcell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base C (3/4);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

GGCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGC  >  NC_002947/2444312‑2444603
                                                                                                                                                   |                                                                                                                                                 
ggCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGcc                                                                                                                                                >  1:418498/1‑150 (MQ=255)
         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                                                                                                                                  >  1:893844/1‑75 (MQ=21)
                                                                                         aCTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCa                                                                                                                                  <  2:893844/75‑1 (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)
                                                                                                                                                   |                                                                                                                                                 
GGCCAAGGCCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGC  >  NC_002947/2444312‑2444603

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: