Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 934,388 C→T P389S (CCT→TCT)  ftsK → cell division DNA translocase FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913934,3880CT98.2% 175.4 / NA 55P389S (CCT→TCT) ftsKcell division DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/29);  minor base A (0/1);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCC  >  NC_000913/934252‑934522
                                                                                                                                        |                                                                                                                                      
gccgccTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTg                                                                                                                                     <  2:17220/140‑1 (MQ=255)
  cgccTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCa                                                                                                                                   <  1:767060/140‑1 (MQ=255)
  cgccTGTTGCCTCTGGTGATGTTCCTCCTCCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCa                                                                                                                                   <  1:273514/140‑1 (MQ=255)
    ccTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGt                                                                                                                                 <  2:77109/140‑1 (MQ=255)
     cTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTg                                                                                                                                <  1:441841/140‑1 (MQ=255)
       gTTGCCTCTGTTGATGTCCCACTTGCGCAACCTACCGTATCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCa                                                                                                                              <  2:225347/140‑1 (MQ=255)
                  tGATGTTCCACCTGCGCAACCTTCAGTAGCCGGGCAGCCGGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGc                                                                                                                   <  2:351743/140‑1 (MQ=255)
                    aTGTTCCACCTGCGCAACCTACCGTCGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCg                                                                                                                 <  1:374632/140‑1 (MQ=255)
                         ccACCTGCGCAACCTCCAGTAGCCTAGCACCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGca                                                                                                            <  1:240832/140‑1 (MQ=255)
                         ccACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGca                                                                                                            >  2:847587/1‑140 (MQ=255)
                         ccACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGca                                                                                                            >  2:633821/1‑140 (MQ=255)
                         ccACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGca                                                                                                            >  2:126804/1‑140 (MQ=255)
                           aCCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAAGATGAGCAGTCTGCAGTACAATTTAATGAGCCGCTGcaac                                                                                                          >  2:329409/1‑140 (MQ=255)
                            ccTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGcaaca                                                                                                         <  1:846776/140‑1 (MQ=255)
                                 gcAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAATCTGCAGAGCAATGGAATGAGTCGCTGCAACAACCAg                                                                                                    >  1:253285/1‑140 (MQ=255)
                                       tACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGc                                                                                              >  1:48393/1‑140 (MQ=255)
                                               ccTGGCAGTCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGGCGCAGCGACAACCGGTACAGCCGcagcag                                                                                      >  1:267254/1‑140 (MQ=255)
                                               ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGACGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATTAGCCGCTGCGACAACCAGTACAGCCAcagtag                                                                                      >  2:816645/1‑140 (MQ=255)
                                               ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCgc                                                                                           >  2:667512/1‑135 (MQ=255)
                                               ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCgc                                                                                           <  1:667512/135‑1 (MQ=255)
                                               ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGcagcag                                                                                      >  1:144883/1‑140 (MQ=255)
                                               ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGcagcag                                                                                      >  2:680846/1‑140 (MQ=255)
                                               ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGcagcag                                                                                      >  1:687814/1‑140 (MQ=255)
                                               ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGcagcag                                                                                      >  1:149691/1‑140 (MQ=255)
                                               ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGcagcag                                                                                      >  1:162900/1‑140 (MQ=255)
                                               ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGcagcag                                                                                      <  2:149691/140‑1 (MQ=255)
                                                 tGGCAGCCTGTACCGTGTCCACTAATGGGCGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGcc                                                                                    <  1:403007/140‑1 (MQ=255)
                                                 tGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGcc                                                                                    <  1:660489/140‑1 (MQ=255)
                                                   gCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGt                                                                                  >  1:146851/1‑140 (MQ=255)
                                                   gCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAAAGAGCCGCGGCAACAGCCAGTAAAGCCGCAGCAGCCGt                                                                                  >  1:413116/1‑140 (MQ=255)
                                                     aGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGtat                                                                                <  1:733118/140‑1 (MQ=255)
                                                               gggTCTACAAACGGGAGTGCCGCTTATTGCTTCTGAACCGGAAGGTTAGCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTg                                                                      <  2:232361/140‑1 (MQ=255)
                                                                 gTCCACAAACGGGAGGGCCGGTCATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCa                                                                    <  1:329409/140‑1 (MQ=255)
                                                                 gTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTg                                                                      >  1:609196/1‑138 (MQ=255)
                                                                 gTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTg                                                                      <  2:609196/138‑1 (MQ=255)
                                                                 gTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCa                                                                    <  2:162900/140‑1 (MQ=255)
                                                                 gTCCACAAACGGGAGAGCCGGTCATTGCTCCGGCACCGGTAGGTTACCTACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCa                                                                    <  2:267254/140‑1 (MQ=255)
                                                                     acaAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTg                                                                >  2:73413/1‑140 (MQ=255)
                                                                                   cGGTTATTGCTCCTGCCCCGGAAGGTTCCCCGCAGTAGTAACAATATGAGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCaa                                                  <  2:311417/140‑1 (MQ=255)
                                                                                   cGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTgcgc                                                    <  1:33699/138‑1 (MQ=255)
                                                                                   cGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTgcgc                                                    >  2:33699/1‑138 (MQ=255)
                                                                                   acagtATTGCTCCTGCCGCGGAGGGTTACCCTCAGCAGTCACAATATGTGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCaa                                                  <  2:312444/136‑1 (MQ=255)
                                                                                      ttATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAg                                               <  2:144883/140‑1 (MQ=255)
                                                                                      ttATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAg                                               <  2:48393/140‑1 (MQ=255)
                                                                                       tATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGc                                              <  1:126804/140‑1 (MQ=255)
                                                                                        tttGCTCCTGCACCGGACGGTTCCCCACAGCAGTCCCAATATGCGCAGACTGCAGTGCAATGTAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGcc                                             <  2:253285/139‑1 (MQ=255)
                                                                                        aTTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGcc                                             <  1:633821/140‑1 (MQ=255)
                                                                                        aTTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGcc                                             <  1:847587/140‑1 (MQ=255)
                                                                                               cTGCACCGGAAGGATACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCCCCTGCAGCTGAACAACCTGCGCAACAGCCGtattat                                      >  2:847415/1‑140 (MQ=255)
                                                                                                  cACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACGACCTGCGCACCAGCCGtattactgt                                   >  1:312444/1‑136 (MQ=255)
                                                                                                   aCCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGTCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGccc                                  <  1:752455/140‑1 (MQ=255)
                                                                                                            ttACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCATCTGCAGCTGAACGACCGGCGCAACACCGGTATTATGCCTCTGCGACac                         >  2:360291/1‑139 (MQ=255)
                                                                                                             tACCCACAGCAGTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGa                        >  1:44956/1‑140 (MQ=255)
                                                                                                                    agcagTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCg                 >  2:787239/1‑140 (MQ=255)
                                                                                                                    agcagTCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCg                 >  2:802985/1‑140 (MQ=255)
                                                                                                                         tCACAATATGCGCAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGc            >  1:118518/1‑140 (MQ=255)
                                                                                                                                   cgcAGTCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGcc  >  1:745083/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAACCGGTGGCAGGTAACGCC  >  NC_000913/934252‑934522

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