Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 933,387 C→T A55V (GCC→GTC)  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_000913933,3870CT98.0% 171.8 / ‑4.1 51A55V (GCC→GTC) ftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (22/28);  total (22/29)
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.25e-01

TGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACAC  >  NC_000913/933291‑933477
                                                                                                |                                                                                          
tGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTg                                                                                         <  1:867022/100‑1 (MQ=255)
  gAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGc                                                                                       <  2:736521/100‑1 (MQ=255)
   aGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCa                                                                                      <  1:598045/100‑1 (MQ=255)
   aGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCa                                                                                      <  2:1509673/100‑1 (MQ=255)
   aGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCa                                                                                      <  1:318658/100‑1 (MQ=255)
       gTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGaa                                                                                  >  2:220171/1‑100 (MQ=255)
       gTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGaa                                                                                  <  1:470008/100‑1 (MQ=255)
             gATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATc                                                                            >  2:1552005/1‑100 (MQ=255)
                 cTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATa                                                                        >  2:628113/1‑100 (MQ=255)
                      tGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTa                                                                   <  2:1493657/100‑1 (MQ=255)
                      tGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTa                                                                   >  2:1173197/1‑100 (MQ=255)
                       gTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAg                                                                  <  1:1445364/100‑1 (MQ=255)
                       gTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAg                                                                  <  2:431708/100‑1 (MQ=255)
                            gTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTggg                                                             >  1:1534703/1‑100 (MQ=255)
                            gTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTggg                                                             <  1:901316/100‑1 (MQ=255)
                            gTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTggg                                                             >  2:791001/1‑100 (MQ=255)
                                gCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGc                                                         >  1:1565525/1‑100 (MQ=255)
                                   gTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCg                                                      >  1:1394174/1‑100 (MQ=255)
                                   gTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCg                                                      <  1:1262063/100‑1 (MQ=255)
                                      tGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTg                                                   >  1:56493/1‑100 (MQ=255)
                                        gTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCg                                                 <  1:1006586/100‑1 (MQ=255)
                                          tGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTg                                               >  2:1228902/1‑100 (MQ=255)
                                           gATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTgg                                              >  1:739957/1‑100 (MQ=255)
                                                cTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGTGTGGTTGGc                                         <  1:1332344/100‑1 (MQ=255)
                                                   ccTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGa                                      <  1:372317/100‑1 (MQ=255)
                                                    cTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGAt                                     >  2:428469/1‑100 (MQ=255)
                                                     ttACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATa                                    >  2:1026514/1‑100 (MQ=255)
                                                        cTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGc                                 <  1:1450735/100‑1 (MQ=255)
                                                          aaGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTg                               <  1:630486/100‑1 (MQ=255)
                                                          aaGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTg                               >  1:13020/1‑100 (MQ=255)
                                                          aaGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTg                               >  2:121326/1‑100 (MQ=255)
                                                           aGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGt                              <  1:473839/100‑1 (MQ=255)
                                                           aGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGt                              >  1:1003032/1‑100 (MQ=255)
                                                             cTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTc                            >  2:1326776/1‑100 (MQ=255)
                                                              tttAACCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTCGGTGGGATGCCCGGTGCGTGGTTGGCGGATACGCTGTTCt                           >  2:1251204/1‑100 (MQ=255)
                                                                ctaCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCttt                         <  1:1175727/98‑1 (MQ=255)
                                                                  aCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTag                       >  2:1175722/1‑99 (MQ=255)
                                                                  aCCCTTCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTAt                       >  1:832227/1‑100 (MQ=255)
                                                                       tCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTg                  <  1:551140/100‑1 (MQ=255)
                                                                       tCGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTg                  <  1:1173202/100‑1 (MQ=255)
                                                                        cGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTgg                 <  1:361294/100‑1 (MQ=255)
                                                                        cGGACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTgg                 <  1:1355916/100‑1 (MQ=255)
                                                                         ggACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGc                <  2:13020/100‑1 (MQ=255)
                                                                          gACCCCAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCg               <  2:407406/100‑1 (MQ=255)
                                                                             cccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGa            >  1:1174303/1‑100 (MQ=255)
                                                                              ccAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGAt           <  1:283397/100‑1 (MQ=255)
                                                                               cAGCTGGTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATg          <  1:758276/100‑1 (MQ=255)
                                                                                     gTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTac    <  2:441409/100‑1 (MQ=255)
                                                                                     gTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTac    <  2:1565515/100‑1 (MQ=255)
                                                                                     gTCGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTac    <  1:791006/100‑1 (MQ=255)
                                                                                       cGCAAACGGTCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTacac  >  2:823475/1‑100 (MQ=255)
                                                                                                |                                                                                          
TGGAGGCGTTGCTGATCCTTATTGTCCTGTTTGCCGTCTGGTTGATGGCTGCCTTACTAAGCTTTAACCCTTCGGACCCCAGCTGGTCGCAAACGGCCTGGCATGAACCTATCCATAATTTAGGTGGGATGCCCGGTGCGTGGTTGGCAGATACGCTGTTCTTTATTTTTGGCGTGATGGCTTACAC  >  NC_000913/933291‑933477

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