Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,602,701 T→C 27.5% Q16R (CAG→CGG)  ftsE ← cell division protein FtsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,602,7010TC27.5% 37.3 / 2.6 27Q16R (CAG→CGG) ftsEcell division protein FtsE
Reads supporting (aligned to +/- strand):  ref base T (5/14);  new base C (3/5);  total (8/19)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTA  >  NC_000913/3602559‑3602836
                                                                                                                                              |                                                                                                                                       
gCCGCTGAATCAGATTTTCCCGGCGCTGTGCCGTTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCt                                                                                                                                         <  2:400539/143‑1 (MQ=255)
    cTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCCGCGCCTGtct                                                                                                                                     >  1:106109/1‑143 (MQ=255)
         ccAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCAc                                                                                                                                <  1:519394/143‑1 (MQ=255)
          cAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGCCCGGTCAGCAACGCCATCTAAGCTGGCCCCACATGGAACGTAACACCCTGCAGCGCCCGTCTcccccc                                                                                                                               >  1:639293/1‑143 (MQ=255)
                              ccGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACAt                                                                                                           <  1:57873/143‑1 (MQ=255)
                                 cTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAAGGTTACGGGGAGAAGCCCCCGTCTCCCACCGACACCACCCTTGCTGCCATgcc                                                                                                        >  1:144072/1‑141 (MQ=255)
                                    aaccccACAGACCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCaa                                                                                                     <  1:346849/140‑1 (MQ=255)
                                           cataTCAGCTTCAGGAGGCTACTTTTCTCTTCTCTATAATGTCCGGTCAGATACGCCATCTCCCCCGGCTGCATATGGTACGTTTTTTTTTTAATGATACGGCGACCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAAt                                                                                              <  2:102495/140‑1 (MQ=255)
                                              accAGCGTCAGGAGGGGACTTTGCCCGGCGCCGGAATGACCGGTCAGAACCGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCAt                                                                                           <  1:502560/141‑1 (MQ=255)
                                               tgaGCGTCAGGAGGGTACTTTGGCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTATCGTTCTGAAGGGCACGGCGACCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCAtt                                                                                          <  2:351516/141‑1 (MQ=255)
                                               ccAGCGTCAGGAGGGTACTTTTCCCTGCGCCGGAATGATCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTTTCGTTTTGCAGGGTCCGGCGACCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCAtt                                                                                          <  2:272282/142‑1 (MQ=255)
                                                 atcTTTAGGAGGGTATTTTTCCCTGCGTCGGAATGTCCTTTCTGATTCTCCATCTCACCCGGCTGCATATGGAACGTTTTTTTTTTAATGGTACGGCGACCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGt                                                                                        <  2:39343/141‑1 (MQ=255)
                                                   cTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGGTa                                                                                      >  1:61575/1‑143 (MQ=255)
                                                               aCTTTTCCCGGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGc                                                                          <  1:277895/143‑1 (MQ=255)
                                                                        tGCGCCGGAATGACCGGTCAGAAGCGCCAGGTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGc                                                                 <  1:116707/143‑1 (MQ=255)
                                                                        tGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGc                                                                 >  1:493451/1‑143 (MQ=255)
                                                                        tGCGCCGGAATGACCGGTCAAAAACGCCATCTGACCCGGCTGCATATGGACCGTAACGCACTGCAGCGCCTGTCTGCCACCGAGATAAGCCTTGCTGACATGGTCAAAGCGAACCATTCGTAATCATCTCGGCCaaaacgagc                                                                 >  2:336513/1‑138 (MQ=255)
                                                                                 aggACCGGTCAGGAACGCCAGCCCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGGCTCCCACCGAGATAAGTCGTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATaaa                                                        <  1:624499/141‑1 (MQ=255)
                                                                                 aTGACCGGTCAGAAACGCCATCTCACCCGGCTGCAGATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATaaa                                                        <  1:328088/143‑1 (MQ=255)
                                                                                  tGACCGGCCAGAAACGCCCTCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAg                                                       <  1:265712/143‑1 (MQ=255)
                                                                                              aaCGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCCGCAGCGCCCGTCTCCCACCTAGCTCAGCCTTGCCGGCCTGATCACAGCCAACCGTTCGCCGTCCTCTAGGGCGAAAAATGGCACTACAACGTCGTCAGccacc                                           >  2:394447/1‑140 (MQ=255)
                                                                                                 gCCATCACACCCGGAAGAATATGGGACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAAc                                        <  1:497405/143‑1 (MQ=255)
                                                                                                            cgtcTGCATATTGAACATAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATc                             <  1:358976/140‑1 (MQ=255)
                                                                                                                   atatGGAACGTAACTCTCTTCAGCGCCCGGCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATa                      <  2:149300/143‑1 (MQ=255)
                                                                                                                                cGCCCCGCAGCGCCTGTCTCCCCCCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGacac         <  1:124132/143‑1 (MQ=255)
                                                                                                                                    cTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAAt     <  1:240698/143‑1 (MQ=255)
                                                                                                                                     tGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCCAAAAGTGCCTCTATAAAGTTGTCCGCCGAAAACCAACGCAGCTCTTTGTCACGATCGCCGACACCAATg    >  1:293717/1‑143 (MQ=255)
                                                                                                                                       cAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTa  <  1:374662/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                       
GCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCGCCGGAATGACCGGTCAGAAACGCCATCTCACCCGGCTGCATATGGAACGTAACGCCCTGCAGCGCCTGTCTCCCACCGAGATAAGCCTTGCTGACATGTTCAAAGCGAATCATTGTTAATCCTCTCGGGCAAAAAGTGCCTCTATAAAGTCGTCCGCCTTAAACGGACGCAAATCCTCAATACGTTCGCCGACACCAATGTA  >  NC_000913/3602559‑3602836

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