Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 935,464 A→G 11.2% A747A (GCA→GCG 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,4640AG11.2% 135.2 / 4.7 54A747A (GCA→GCGftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base A (17/30);  new base G (4/2);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGC  >  NC_000913/935321‑935591
                                                                                                                                               |                                                                                                                               
ggCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCgc                                                                                                                             <  5:41859/148‑1 (MQ=255)
 ctgAACAACCGGCTTGGGCGAATCCGTTCTCGCTGGTTGATTTTTAATTTTCGCCAATGAAAGCGTTGCTTGATTATGGTCCACACGAACCGTTGTTTAAGCCAATTTTTTAACCTGTACAGCAGCCTCAACAACCGGTTGCACCgc                                                                                                                             <  8:410679/145‑1 (MQ=255)
    aacaacCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCgcagca                                                                                                                         <  4:121166/148‑1 (MQ=255)
    aacaacCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGGTGTTTACGCCAATTGTTGAACCGGTACAGCAGCCGCAACAACCCGTTGCACCgcagcc                                                                                                                         >  3:229098/1‑147 (MQ=255)
     acaacCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGTATGATGGTCAAAACGAACCGTTGTTTACGCCAATTGTTGAACATGTACAGCAGCCGCAACAACCGCTTGCACCgca                                                                                                                            >  7:410679/1‑144 (MQ=255)
       aacCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCgcagcagca                                                                                                                      <  1:273372/148‑1 (MQ=255)
        acCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCCCCCCAGCAGCaa                                                                                                                     >  3:264474/1‑148 (MQ=255)
          cGGCTGGGGCGAATGCGTTCTCTCTGGATGATTTTGAATTTTAGCCAATGAAAGCGTTGCTGGATGATGGTCCACTCGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAata                                                                                                                   <  6:102707/148‑1 (MQ=255)
          cGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAata                                                                                                                   <  8:417117/148‑1 (MQ=255)
          cGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAata                                                                                                                   <  7:159517/148‑1 (MQ=255)
            gCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATc                                                                                                                 <  7:214348/148‑1 (MQ=255)
            gCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTGGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATc                                                                                                                 <  7:168003/148‑1 (MQ=255)
                gggCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATcagca                                                                                                             <  7:356823/148‑1 (MQ=255)
                gggCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATcagca                                                                                                             <  5:403477/148‑1 (MQ=255)
                gggCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATcag                                                                                                               >  4:134291/1‑146 (MQ=255)
                    gAATCCGTTCTCGCTGGATGATTTTGAATTTGCGCCAATGAAAGCGGTGCGGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCg                                                                                                         <  6:412547/148‑1 (MQ=255)
                      aTCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCCTCAATATCAGCcgccgc                                                                                                        >  7:316851/1‑147 (MQ=255)
                        ccGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAATATCCGCAGCCGcaa                                                                                                      >  5:274768/1‑147 (MQ=255)
                           ttCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAGTATCAGCAGCCGCaacaac                                                                                                  >  3:40219/1‑148 (MQ=255)
                            tctcGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCaacaac                                                                                                  >  7:127090/1‑147 (MQ=255)
                                gCTGGATGTTTTTGTACTTTCGCCTCTGAATCCGTTGCGGGAGGATGGTCCTCACGAAGCGCTGTTGACGGCTATTGTGGAACCTGTACAGCAGGCGCATCAACCGGTTGCGCCTCAGCAGCAATATCAGCAGCCGCAACAAccg                                                                                                <  4:264474/145‑2 (MQ=16)
                                gCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGt                                                                                              <  3:394672/147‑1 (MQ=255)
                                  tGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTcc                                                                                           >  5:42591/1‑148 (MQ=255)
                                     atgatTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTccgcc                                                                                        >  1:68086/1‑148 (MQ=255)
                                     atgatTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCCGCAATATCAGCAGACGCAACAACCACTTccgcc                                                                                        >  8:435731/1‑148 (MQ=255)
                                        atTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTccgccg                                                                                       >  7:294476/1‑146 (MQ=255)
                                               aTTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGccgcagccgcag                                                                              <  6:42591/148‑1 (MQ=255)
                                                  ttCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTa                                                                            <  8:329419/147‑1 (MQ=255)
                                                    cGCCTATGAAAGCGTTGCTGGATGATGGTCGACACGAACCGTTGTTTACGGCAATTGTTGTACCTGTACAGCAGCCGCAACAACCTGTTGCACCGCAGCAGCAATATCAGCAGCGGCAACAACCAGTTCCGCCGCAGCCGCTGTATc                                                                          <  1:216831/147‑1 (MQ=255)
                                                    cGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATc                                                                          <  3:184640/147‑1 (MQ=255)
                                                     gCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATcag                                                                        <  8:294476/148‑1 (MQ=255)
                                                      ccAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATcag                                                                        >  8:15630/1‑147 (MQ=255)
                                                          tGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGcc                                                                   <  8:127090/148‑1 (MQ=255)
                                                           gAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAATATCAGCAGCCGCACCCACCAGTTCCGCCGCAGCCGCAGTATCAGCCGCca                                                                  >  3:313449/1‑148 (MQ=255)
                                                                        gatgatGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCca                                                                  >  1:196395/1‑135 (MQ=255)
                                                                        gatgatGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCca                                                                  <  2:196395/135‑1 (MQ=255)
                                                                        gatgatGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGt                                                       <  3:134291/146‑1 (MQ=255)
                                                                            atgGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGcgcc                                                 >  5:161916/1‑148 (MQ=255)
                                                                               gTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTAt                                                                           >  5:235351/1‑119 (MQ=255)
                                                                               gTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTAt                                                                           <  6:235351/119‑1 (MQ=255)
                                                                                  cacaCGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGcc                                           <  7:15630/148‑1 (MQ=255)
                                                                                    cacGAACCGTTGGTTACGCCAATTGTTGAACATGTCCAGCAGCAGCAACAACCGTTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCACCTCCTCCGCAGCCTCATTATCAGCAGCCACAACATCCGGTTGCGCCGCAGCcac                                         <  4:313449/148‑1 (MQ=21)
                                                                                     acGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCTCAACAGCCAGTTGCGCCGCAGCCGCAGGATCAGCAGCAACAACCGCCGGTTACGCCGCAGAcaca                                        >  3:76157/1‑148 (MQ=21)
                                                                                       gAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAat                                      <  6:161916/148‑1 (MQ=255)
                                                                                          ccGGTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCTCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATc                                   <  7:435731/148‑1 (MQ=255)
                                                                                                tttACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTCCGCCCCAGCCGCAGTAACAGCAGCCACAACAGCCGGTTGCGCCTCAGCCCCAAGATcaccag                              >  4:398519/1‑147 (MQ=21)
                                                                                                        aaTAGTTCAACCTGACCATCACCCGCAACAACAGGATGCACCGCAGCATCAATATCAACAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGcaacaa                     <  1:67183/148‑1 (MQ=21)
                                                                                                              tgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATc                                                                          <  3:253742/89‑1 (MQ=255)
                                                                                                              tgAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATc                                                                          >  4:253742/1‑89 (MQ=255)
                                                                                                                    tGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAg                                                            >  2:252291/1‑97 (MQ=255)
                                                                                                                    tGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAg                                                            <  1:252291/97‑1 (MQ=255)
                                                                                                                        cagcagCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagc     <  4:178755/148‑1 (MQ=255)
                                                                                                                            agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCCGCCGCACCAACCGGTTGCGCCACagcagc  >  5:387903/1‑147 (MQ=255)
                                                                                                                                     aacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACATCCTGTTGCGCCGCAGCCGCAGTATCAGCAGCGACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGt                <  4:195925/124‑1 (MQ=255)
                                                                                                                                               |                                                                                                                               
GGCGAACAACCGGCTGGGGCGAATCCGTTCTCGCTGGATGATTTTGAATTTTCGCCAATGAAAGCGTTGCTGGATGATGGTCCACACGAACCGTTGTTTACGCCAATTGTTGAACCTGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGC  >  NC_000913/935321‑935591

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