Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 935,676 C→T 10.0% T818I (ACC→ATC)  ftsK → DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913935,6760CT10.0% 113.9 / 4.6 40T818I (ACC→ATC) ftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base C (16/20);  new base T (2/2);  total (18/22)
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.70e-01

TTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACAA  >  NC_000913/935538‑935812
                                                                                                                                          |                                                                                                                                        
ttGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGGTATCCTGcggttg                                                                                                                                <  3:182155/149‑6 (MQ=255)
             acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTc                                                                                                                                  >  1:274144/1‑134 (MQ=255)
             acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTc                                                                                                                                  <  2:274144/134‑1 (MQ=255)
             acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGc                                                                                                                   >  1:226271/1‑149 (MQ=255)
             acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGc                                                                                                                   >  3:150173/1‑149 (MQ=255)
                      gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCg                                                                                                          >  3:288918/1‑149 (MQ=255)
                      gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCg                                                                                                          >  3:214119/1‑149 (MQ=255)
                            gCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGcc                                                                                                    >  3:204008/1‑149 (MQ=255)
                                 aacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccg                                                                                               >  1:116505/1‑149 (MQ=255)
                                 aacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccg                                                                                               >  3:400818/1‑149 (MQ=255)
                                  acCGGTTGCGCCGCAGCAGCAATATCAGCAGCCGCAACAACCCGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccgt                                                                                              <  2:219792/149‑1 (MQ=255)
                                   ccGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGtt                                                                                             <  3:222062/149‑1 (MQ=255)
                                      gTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa                                                                                          <  3:293360/149‑1 (MQ=255)
                                           gcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAc                                                                                     >  3:464220/1‑149 (MQ=255)
                                           gcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAc                                                                                     <  2:116505/149‑1 (MQ=255)
                                                    gcaATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGAcgc                                                                            <  2:429963/149‑1 (MQ=255)
                                                             gcagcCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAc                                                                                     >  3:303735/1‑131 (MQ=255)
                                                             gcagcCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAAGCGTAATGGCGACAGCCGTCCGTTGCATAAAc                                                                                     <  4:303735/131‑1 (MQ=255)
                                                                 ccGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTctct                                                               >  1:267551/1‑149 (MQ=255)
                                                                   gCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCCCAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTgg                                                             <  4:400818/149‑1 (MQ=255)
                                                                    caacaaCCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGa                                                            <  2:137692/149‑1 (MQ=255)
                                                                         acCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGc                                                       <  2:272345/149‑1 (MQ=255)
                                                                         acCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGc                                                       <  2:203459/149‑1 (MQ=255)
                                                                                       agcagcagTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACAccgccgccg                                         >  2:88383/1‑149 (MQ=255)
                                                                                        gcagcagTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGa                                        <  4:150173/149‑1 (MQ=255)
                                                                                                   agcagcCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAg                             <  4:214119/149‑1 (MQ=255)
                                                                                                    gcagcCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGc                            >  3:111029/1‑149 (MQ=255)
                                                                                                    gcagcCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGc                            <  3:207767/149‑1 (MQ=255)
                                                                                                    gcagcCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGc                            >  1:277260/1‑149 (MQ=255)
                                                                                                       gcCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCgg                         <  3:1397/149‑1 (MQ=255)
                                                                                                           caacaGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGa                     <  2:277260/149‑1 (MQ=255)
                                                                                                              caGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATAc                  <  4:288918/149‑1 (MQ=255)
                                                                                                                gCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGCGCCGGTAGATACCt                <  2:267551/149‑1 (MQ=255)
                                                                                                                gCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCt                <  3:457641/149‑1 (MQ=255)
                                                                                                                        cgcCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCtt        <  4:464220/149‑1 (MQ=255)
                                                                                                                         gcCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTg       >  4:472917/1‑149 (MQ=255)
                                                                                                                            acaACCGCAGGATATCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTctct                                                               <  2:285819/90‑1 (MQ=255)
                                                                                                                            acaACCGCAGGATATCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTctct                                                               >  1:285819/1‑90 (MQ=255)
                                                                                                                            acaACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAAc    >  1:305450/1‑149 (MQ=255)
                                                                                                                              tgaagCAGGGTATCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACaa  >  1:477429/5‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                        
TTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCGACGACGCCGCTGCCTTCTCTGGATTTGCTGACACCGCCGCCGAGCGAAGTGGAGCCGGTAGATACCTTTGCGCTTGAACAA  >  NC_000913/935538‑935812

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