Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,825 G→C 11.2% P88A (CCC→GCC)  ftsX ← cell division protein FtsX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,8250GC11.2% 106.1 / 4.4 43P88A (CCC→GCC) ftsXcell division protein FtsX
Reads supporting (aligned to +/- strand):  ref base G (14/24);  new base C (3/2);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

TCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAG  >  NC_000913/3601686‑3601967
                                                                                                                                           |                                                                                                                                              
tCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAg                                                                                                                                             >  2:151812/1‑143 (MQ=255)
  gtccTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGGCTTTGTACACCATATCACCGACGCTGTGCAGCg                                                                                                                                           >  2:899092/3‑143 (MQ=255)
      tccAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCag                                                                                                                                       <  1:747432/141‑1 (MQ=255)
        agaCTATGCCACCACGCCAGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGCCCAGCGTctgtt                                                                                                                                     <  1:274324/141‑5 (MQ=255)
         aaCTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAgaga                                                                                                                                    <  1:900363/143‑1 (MQ=255)
           cTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTACCGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAgcgaca                                                                                                                                  >  1:110575/1‑141 (MQ=255)
              tgCCACCGCGCCCGCAGCAGCGGCATCGGCCAGCGGTTTTTGCAGAGAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATgg                                                                                                                               <  1:252346/143‑1 (MQ=255)
              tgCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGATCTTTGCCGCTAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCGGCGTCAGAGAAACgg                                                                                                                               >  2:498496/1‑143 (MQ=255)
                ccaccaCGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGCTTTGCGGTGACGGATAATGCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGAAGCTGGGCAGCGTCAGAGAAATGGCg                                                                                                                             >  1:888806/1‑143 (MQ=255)
                 cacgacGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGa                                                                                                                            <  2:670861/143‑1 (MQ=255)
                 caccacGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGa                                                                                                                            >  1:435050/1‑143 (MQ=255)
                      cGCCCGCAGCAGCGTCATCGTCCAGCGTTTTGTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAAcc                                                                                                                       <  2:523942/143‑1 (MQ=255)
                                    tCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTaaaaa                                                                                                         >  1:81708/1‑143 (MQ=255)
                                     cATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAAc                                                                                                        >  2:733263/1‑143 (MQ=255)
                                          tCCAGCGTTTTTTGCAGATAAGCAGTGATTTGCTGTGACGGATCATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGc                                                                                                   <  2:117548/143‑1 (MQ=255)
                                             aGCGTTTTTTGCAGATAAACCGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCgcc                                                                                                >  1:316992/1‑141 (MQ=255)
                                              gCGGTTTTTGCAGATAGACAGTGATTTGCGGTGACGGATAATTCTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACATACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATTACCGTTAAAAACGTGGCGAAc                                                                                               <  2:417054/143‑1 (MQ=255)
                                                  tttttGCAGATAGACAGTGATTTGCGGTGACGGATAATACTGGGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGtt                                                                                           <  1:887098/143‑1 (MQ=255)
                                                           aTAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAACCGTGGCGAACGGTTTGCTTTTCa                                                                                  <  1:28043/143‑1 (MQ=255)
                                                           aTAAACAGTGATTTCCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTGGTACACCATATAGCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCa                                                                                  <  2:519789/143‑1 (MQ=255)
                                                                 aGTGATTTGCGGTGACGGATAAGACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCt                                                                            <  1:727355/143‑1 (MQ=255)
                                                                  gTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTg                                                                           <  2:421589/143‑1 (MQ=255)
                                                                  gTGAGTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTg                                                                           <  1:878537/143‑1 (MQ=255)
                                                                   tGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGCCGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGc                                                                          <  1:572565/143‑1 (MQ=255)
                                                                     cgttGCGGGCACGGATAATCCTTCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGGTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCaa                                                                        <  1:512230/141‑1 (MQ=255)
                                                                             tGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGt                                                                <  1:733263/143‑1 (MQ=255)
                                                                             tGACGCATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCCTATAACAGACGCTGGCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGt                                                                <  2:672063/143‑1 (MQ=255)
                                                                                  gATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAg                                                           <  1:151812/143‑1 (MQ=255)
                                                                                                     cgggTTCACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGt                                        <  2:428974/141‑1 (MQ=255)
                                                                                                                ttttGTACACCATATAGCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGTAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTg                             <  2:84928/143‑1 (MQ=255)
                                                                                                                ttttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTg                             <  1:303975/143‑1 (MQ=255)
                                                                                                                ttttGTACACCAGATAACAGCCGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCCCCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGCTCCTGCGATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTg                             <  2:84175/143‑1 (MQ=255)
                                                                                                                  ttGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACggtggt                           >  1:233454/1‑143 (MQ=255)
                                                                                                                        acCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCg                     >  1:45921/1‑143 (MQ=255)
                                                                                                                             ataACAGTCGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCTCCTGTCAGTTGAAAACGTTGGTTTTGCGATTTa                <  2:110575/143‑1 (MQ=255)
                                                                                                                               aaCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGCTTTAcc              >  1:484337/1‑143 (MQ=255)
                                                                                                                                tctgACGCTGCCCAGGGTCAGAGAAATGGCGATAACCATCACCGTTCAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCACTGCGCCGTGGAAGGCCTAGCGCACCTGTTCGTTGAAACCGTTGGTCTTCCGAATTAcag             >  2:410705/4‑141 (MQ=255)
                                                                                                                                tctgACGCTGCCCAGCGTCTGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGCGGAAGCCATAGCGCCCCAGTACGTTGAAAACGTTGGTTTTGCGAtttcacg             >  2:100579/4‑139 (MQ=255)
                                                                                                                                tctgACGCTGCCCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGCTGCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGTTTTACCg             >  1:247036/4‑143 (MQ=255)
                                                                                                                                aCAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTAACg             >  1:242353/1‑143 (MQ=255)
                                                                                                                                     cGCTGGGCGGCGTCAAAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGtttt        <  1:899092/143‑1 (MQ=255)
                                                                                                                                      gCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGGTTGCTTTTCAGATCCTGCAATGCGCCGTGGACGGCATAGCGCACCTGTTCGTTGAGAACGTTGGTTTTGCGATTTACCGGTtgtg       >  1:348080/1‑143 (MQ=255)
                                                                                                                                      cgtGGGCAGCGCCAGCTAATTGGCCATAACGATCACCGTTAACAACGTGGCGAACGGTTTTCTTTTCATATTCTGCAATGCGCCTTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTg       <  2:84348/141‑1 (MQ=255)
                                                                                                                                           gCAGCGTTCGATAAAGGGCGATCACCATCACCGTTAAAAACGTGGCGAACGGGCTGCCTTTCAGATCCTGCAATGCTCCGTGGAAGGCATAGCTCACCTGTTCGTTGAAAACGTCGGTTTTGCGATTTACCGGTTTTGGCgcg  <  1:70031/143‑3 (MQ=255)
                                                                                                                                           gCAGCGTCATGGTAATGTCGATACCCATCACCTTTGAAAAGGTGGAGAAGGGTTTGGTGTGCAGATCCTGCAATGCGCCGGGGAAGGCATAGCGCACCGGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAg  <  1:147482/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
TCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGTCGCCGCCTGGTTAACGTTTTTGTACACCATATAACAGACGCTGGGCAGCGTCAGAGAAATGGCGATAACCATCACCGTTAAAAACGTGGCGAACGGTTTGCTTTTCAGATCCTGCAATGCGCCGTGGAAGGCATAGCGCACCTGTTCGTTGAAAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAG  >  NC_000913/3601686‑3601967

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