Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 91,904 G→A E164E (GAG→GAA)  ftsI → peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091391,9040GA96.5% 87.2 / ‑3.7 29E164E (GAG→GAA) ftsIpeptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (10/18);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

GGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCG  >  NC_000913/91661‑92138
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
ggTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGcc                                                                                                                                                                                                                                     <  2:27065/251‑1 (MQ=255)
    aGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGACAGACCTGCCTGACGCTGGCGGTATGAGCGTCGGTGACCGCTGGAAGGCACTGGATAACGCGCTCAATATTCCGCTGGATCGGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTCTTTATCTGGCGCGCCCGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATATGCGGGAAGAATATCGCCGtt                                                                                                                                                                                                                                  >  1:29094/1‑250 (MQ=255)
     gCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTat                                                                                                                                                                                                                                <  1:233054/251‑2 (MQ=255)
                                            ggCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCttt                                                                                                                                                                                         <  2:99312/250‑1 (MQ=255)
                                              cATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTAc                                                                                                                                                                                       <  1:123769/251‑1 (MQ=255)
                                                      tGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCg                                                                                                                                                                               >  1:246300/1‑251 (MQ=255)
                                                                            cGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCCAGGGAGTGAGGGCGt                                                                                                                                                         <  2:210241/251‑1 (MQ=255)
                                                                                      cgcTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGt                                                                                                                                               >  2:111161/1‑251 (MQ=255)
                                                                                          ggCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCg                                                                                                                                           <  2:269670/251‑1 (MQ=255)
                                                                                              aaCGCGCTCAAGATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGGCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCGCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTTAGAAGAGTTTCGAt                                                                                                                                         <  2:234987/249‑1 (MQ=255)
                                                                                               aCGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCTGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATaaa                                                                                                                                      <  1:1368/251‑1 (MQ=255)
                                                                                                cgcgCTCAATAGTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAACCTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAAt                                                                                                                                     <  2:151406/251‑1 (MQ=255)
                                                                                                                                cccGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATtgtg                                                                                                     <  1:10482/251‑1 (MQ=255)
                                                                                                                                 ccGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGc                                                                                                    <  1:226636/251‑1 (MQ=255)
                                                                                                                                           gCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCg                                                                                          >  1:169205/1‑251 (MQ=255)
                                                                                                                                                  cGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTACCGTCGATAGCCAAGGGATTAGGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGt                                                                                   >  2:236441/1‑251 (MQ=255)
                                                                                                                                                      aGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTcgcg                                                                               <  1:111161/251‑1 (MQ=255)
                                                                                                                                                              ttatttatCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGTGCGCATTGTGCGTAAAGAACGCTATGGTCGCGTAATTGa                                                                        >  2:205715/1‑250 (MQ=255)
                                                                                                                                                               tatttatCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGaa                                                                       <  1:218323/250‑1 (MQ=255)
                                                                                                                                                                atttatCTGGCGCGTCAGGGGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGa                                                                     <  2:144882/251‑1 (MQ=255)
                                                                                                                                                                                       cccTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAgg                                              >  2:76398/1‑251 (MQ=255)
                                                                                                                                                                                                  cGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAAc                                   >  1:216481/1‑251 (MQ=255)
                                                                                                                                                                                                            cAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTg                          >  1:121278/1‑250 (MQ=255)
                                                                                                                                                                                                             aaaaaaCTGAAACTGCCGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAg                        <  2:154409/251‑1 (MQ=255)
                                                                                                                                                                                                                            ccGGGGATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCt         <  1:193357/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggggATTCATCTGCGTGAAGAATCTCGCCGTTCTTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGc       >  2:161744/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAg     <  1:10935/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gATTCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAgg    >  2:218168/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCATCTGCGTGAAGAATCTCGCCGTTATTATCCATCCGGCGAAGTGACTGCTCACCTCATTGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACATCCAGGCTGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGgcg  <  1:248307/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
GGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCTGGCGCGTCAGGTGAACCCTGACATGGCGGACTACATCAAAAAACTGAAACTGCCGGGGATTCATCTGCGTGAAGAGTCTCGCCGTTACTATCCGTCCGGCGAAGTGACTGCTCACCTCATCGGCTTTACTAACGTCGATAGTCAAGGGATTGAGGGCGTTGAGAAGAGTTTCGATAAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCG  >  NC_000913/91661‑92138

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑