Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,602,498 T→C 16.0% I84V (ATT→GTT)  ftsE ← cell division protein FtsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,602,4980TC16.0% 44.4 / 4.0 25I84V (ATT→GTT) ftsEcell division protein FtsE
Reads supporting (aligned to +/- strand):  ref base T (10/11);  new base C (2/2);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

GCTTTGTCCAGTAGCCCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCG  >  NC_000913/3602356‑3602634
                                                                                                                                              |                                                                                                                                        
gCTTTGTCCAGCAGCCCGACTGCATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAAt                                                                                                                                          <  2:378112/143‑1 (MQ=255)
    ggTCCAGTAGCCCGACTGTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATg                                                                                                                                      <  1:256788/142‑1 (MQ=255)
          gTAGCCCGACTTTATCCGGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATc                                                                                                                                <  1:559359/143‑1 (MQ=255)
               ccGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTggcg                                                                                                                           >  1:383451/1‑143 (MQ=255)
                    tttATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGCa                                                                                                                      <  1:683063/143‑1 (MQ=255)
                          cAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTCGATGGTCATCCTGGAAAATCATGCCTATCTAGCGACAAAGaaaaa                                                                                                                >  1:387451/1‑141 (MQ=255)
                          cAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACg                                                                                                                <  1:577443/143‑1 (MQ=255)
                                  ccgACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCTGTCAGTAGATGGTGATCCTGGAAAACCATGCCAATCTGGCGGCGCAGAAACGGAACTTCa                                                                                                        >  2:493524/1‑143 (MQ=255)
                                  ccgACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAACCATGCCAATCTGGCGGCGCAGAAACGGAACTTCa                                                                                                        >  2:124727/1‑143 (MQ=255)
                                         ccGGCGACGAATATCGTCACCGCAGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACGGTTCGGTCCCTGAGTAGATGGTGATCCTGGAAAATCATGCCAATAACGCGTCGCGGAAACGCAACATCACGGtata                                                                                                 >  2:817152/1‑140 (MQ=255)
                                            gccacgaATACCATCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCa                                                                                              <  1:629451/140‑1 (MQ=255)
                                                        gTCACCGCTGGCACCGGCGATAATCAACGGGATCGCCACGGGATCGTCGACAGTACGGCCCATCAGTAGATGGTGATCCTGGAAAATCATGCAAATCTAGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGt                                                                                  >  1:35781/1‑143 (MQ=255)
                                                                      cGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTCATGGCCGCTAAAc                                                                    <  2:26121/143‑1 (MQ=255)
                                                                      cGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGCACGGTCCATCCGTAGATGGTGATCCTGGAAAATCATACCAATCTGGCGGCGCAGAAACGGGACTTCACGGTTGCTCAGACGCCTGCGGTCATGGCCGcgcaac                                                                    >  2:115740/1‑138 (MQ=255)
                                                                         cgtgaaTCAGCGGGATCCTCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAACTCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTCATGGCCGCTAAACCAg                                                                 <  2:650301/139‑1 (MQ=255)
                                                                                        tCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTGCACGGGTTTTCAGACGCGTGATGTCATGGCCGCTAAACCAGATTTTCCCGGCGCTg                                                  >  1:689594/1‑143 (MQ=255)
                                                                                              cGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGc                                            >  1:48643/1‑143 (MQ=255)
                                                                                                                         gTAGATGGTGATCCTGGAAAATCATGGCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTTCTGTCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCaggagg                 <  1:518800/143‑1 (MQ=255)
                                                                                                                          tAGATGGTGATCCTGGAAACCCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAggg                <  1:124727/143‑1 (MQ=255)
                                                                                                                          tAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTGTTACAGACGCGTGATGGCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTGAATCCCACAGAACAACAACAGGAggg                >  2:435631/1‑143 (MQ=255)
                                                                                                                          tAGATGGTGATCCTGGAAAACCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAggg                <  1:493524/143‑1 (MQ=255)
                                                                                                                              tGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCCCTCAATCCCACACATCAGCTTCAgggggggcct            >  1:811078/1‑139 (MQ=255)
                                                                                                                                  gATCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTc        <  1:382254/143‑1 (MQ=255)
                                                                                                                                   aTCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGTGTGATGTCATGGCCGCTAAACGAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTcc       <  1:625485/143‑1 (MQ=255)
                                                                                                                                       tGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTCATACCCCCTCAACCAGTTTTTGGCACAACTTGGGCTCTCAATGCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTgc   >  1:153893/1‑143 (MQ=255)
                                                                                                                                        tgagaaTCATGCCAATCTGGCGGCGCGGAAACGGAACTTCACGGTTTTGCAGACGCGTGATGTCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTgcg  <  2:387451/139‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                        
GCTTTGTCCAGTAGCCCGACTTTATCCAGCGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCGGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCATCAGTAGATGGTGATCCTGGAAAATCATGCCAATCTGGCGGCGCAGAAACGGAACTTCACGGTTTTTCAGACGCGTGATGTCATGGCCGCTAAACCAGATTTTCCCGGCGCTGGGCCGCTCAATCCCACAGATCAGCTTCAGGAGGGTACTTTTCCCTGCG  >  NC_000913/3602356‑3602634

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: