Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,508,654 T→G 33.2% L178R (CTG→CGG)  ftsI → Peptidoglycan D,D‑transpeptidase FtsI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,508,6540TG33.2% 20.7 / 5.3 37L178R (CTG→CGG) ftsIPeptidoglycan D,D‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base T (10/13);  new base G (9/4);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGCTCGACCTCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGCGTTGTCGATAGACC  >  KT2440/1508523‑1508788
                                                                                                                                    |                                                                                                                                      
cGTGCGCGACCTCAAAGTGCCTGGTGGGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGGGACTGCGCACAGGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGc                                                                                                                                     <  2:306790/135‑1 (MQ=255)
     tCGACCTCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACCTGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGGAGGGGTGGAGCGGGCCTAt                                                                                                                                 >  5:279913/1‑134 (MQ=255)
     tCGACCTCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATg                                                                                                                                >  8:445937/1‑135 (MQ=255)
       gACCTCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGGCGGGGGGGGGGTGGCCttt                                                                                                                                 >  5:99055/1‑130 (MQ=255)
           tCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGGAGGGGTGGAGCGGGGCGATGAAGAGt                                                                                                                          >  3:95407/1‑135 (MQ=255)
                gTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGGACGGGTGGGGGTGGGCTATTAATAGttggtt                                                                                                                     >  3:281994/1‑130 (MQ=255)
                     tGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGGAGGGGGGGGGGGGGGCGATGAAGAGTGGCTGGGCgg                                                                                                                >  1:117361/1‑135 (MQ=255)
                      ggTGTGTACGGACTGGAAGACTTCCGCCGCTGCTATCCGGCCGGTGATGTGCCTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGGATGGGGGGGGGGGGGCGATGAAGAGTGGCTGGCCGGcc                                                                                                              >  7:184844/1‑135 (MQ=255)
                                    ggAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGGTGGGGTGGGGCGGGGCGATGTAGCGTGGCTGGCCGGGGGGCCCCGccaccg                                                                                                 >  5:349210/1‑130 (MQ=255)
                                    ggAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGCAGGGGGGGAGCGGGGCGATGAAGACTGGCTGGCCGGCGGGCCCCGCCAgcg                                                                                                 >  4:369167/1‑135 (MQ=255)
                                    ggAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGAAGGGGTGGCGCTGGCCAATGCACACTGGCTGGCCCGCGGGCCCCGccaccc                                                                                                 >  4:345281/1‑130 (MQ=255)
                                    ggAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGAAGCGGTCGAGCAGGCCTATGAAGACTGGCTGGCCGGCCTGCCCGGCAAacg                                                                                                 >  6:104677/1‑135 (MQ=255)
                                     tattaaTTCCGCGTCTCATAGCGGGCGGGTGATGTGACTCCGACCATGGGGGGCTTTA‑CCGACGGCCACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATCAAGAGTGGCTGGCCGGCGTGCCCGGAAAGCGcc                                                                                               <  8:96327/132‑1 (MQ=255)
                                      aataaTGCCGCCGCTTCAATCCGGCCGGTTCGGTCACGGCGAACAGGGGGGGGTTAC‑CCCACCGCGACGACCACGGGCGTGAAGGGGTGGAGCTGGCCTATGAAGAGTGGCGGGCCGGCGTGCCCGGCAAGCGcc                                                                                               <  6:279913/132‑1 (MQ=255)
                                      aagaaTTCCGCCGTTTCTAACCGGCCGGTGATGGGCCTGGGCACAGGGGCGGCTTTC‑CCGACCGCGACGACCACGGTCGTGATGGGGTGGAGCTGGCCTATGACGAGGGGCTGGCCGGCGTGCCCGGCAAGCGcc                                                                                               <  6:122510/135‑1 (MQ=255)
                                         aaTTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGGAGGGGTGGCGCTGGGCGATGTAGAGTGGCTGGCCGGCGTGGCCGGCCAGCGCCCgg                                                                                            >  1:92985/1‑135 (MQ=255)
                                          tttCCGCCGCTTCTTCCGGGCCGGGGAGGTGACTGCGCACTTGGTGGGCTTTC‑CCGACCGCGACGACCACGGTGGTGATGGGGTGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTg                                                                                          <  7:445937/135‑1 (MQ=255)
                                              cgccgcTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGCAGGGGTGGGGCGGGCCTATGAAGAGTGGCTGGCCCGCGGGGCCGGCACGCGCCAGGTGAt                                                                                        >  2:342950/1‑134 (MQ=255)
                                                             gccggTGATGTGCCTGCGCACTTGGGCGGCTTAC‑CCGACGTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCAATGAAGAGTGGCTGGCCGGGGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGt                                                                        <  1:83290/135‑1 (MQ=255)
                                                                      gttcgtGCGCACATGGT‑GGTTTTACCCTACGTGGACGCCCTCGGGCGGGGCGGGGGGGAGCTGGCCTATGAAGAGGGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGAt                                                                <  1:57690/129‑1 (MQ=255)
                                                                            gcgcACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTg                                                                                          >  8:323015/1‑102 (MQ=255)
                                                                            gcgcACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTg                                                                                          <  7:323015/102‑1 (MQ=255)
                                                                                        ggCTTTC‑CCGACGTCGAGGACCACGGTGGTTATGGGTGGGTGGTGGCGTTTGATGTGTGGGTGGCCGGGGTGCCGGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAAcc                                             <  2:241634/135‑1 (MQ=255)
                                                                                          cTTTA‑CCGACCTCGACGACCACGGTCGTGCAGGGGTGGGGCGGGCCTATGAAGAGTGGCTGGCCCGCGTGCCCGGCCAGCGCCAGGTGAGCAAGGCCCGGCGCGGGCGGCTGATGAAGGGCATTGAGGTAAc                                              >  4:319154/1‑132 (MQ=255)
                                                                                                ccGACCTCCACGCCCCCGGTCGTGAAGGGGTGCCGCTGGCCTATGAAGCGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACg                                      <  3:117569/135‑1 (MQ=255)
                                                                                                 cgacggcgacgaCCACGGTCGTTAAGGGGGGGAGCGGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACg                                      <  4:95407/128‑1 (MQ=255)
                                                                                                        tagccccAGGGTCGTTAAGGGGGGGAGCTGGCCTATGACTTGGGGCGGGCGGGGGGGCCCGGAAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCACGGCCg                             <  7:294709/131‑1 (MQ=255)
                                                                                                        gacgacCACGGTCGTGAAGGGGGGGAGCGGGCCGATGAACTGTGTCTGGGCCGCGGGGCCGGCCCGCGCCACGTGATGAACGCCCCGCGCGGGGGGGGGAGCACGGCTCTTGAGGGAACCAAGAACGCCacagcc                              >  6:32000/1‑130 (MQ=255)
                                                                                                         acgacCACGGTCGTGACGGGGTGGAGCGGGCCTATGAAGAGGGGCTGGCCGGCGGGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCg                             <  3:241816/135‑1 (MQ=255)
                                                                                                                      tGATGGGGGGGCGCGGGCTTATTACGAGTGGCGGCCGGGCGTGCCCGGCAAGCGCCAGGTGAGCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGGACCCAAGAACGCCAAGGCCGGGAAGACCTTGGc                <  3:72434/135‑1 (MQ=255)
                                                                                                                      tGACGGGGGGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGc                <  6:220986/135‑1 (MQ=255)
                                                                                                                      ggacggggTGGAGCGGGCCTATGAAGAGTGGCTGCCCGGCGGGCCCGTCAAGCGCCAGGTGATCACGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGc                <  4:31449/131‑1 (MQ=255)
                                                                                                                       tatggggTGGCGCTGGCCTATGACGAGTGGCTGGCCGGGGGGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGc                <  7:78741/131‑1 (MQ=255)
                                                                                                                                gAGCTGGCCTATGAAGATTGGCTTGCAGACGTGCCCGGAAAGCGACAGGTGATCAAGGACCGGCGCGGTCGGCTGATAAAGGACATTCAGGTAACCAAGAAAGCCAATGCAGGGAAGACCTTGTCGTTGTCGAt       >  6:420067/1‑134 (MQ=255)
                                                                                                                                gAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGCGTTGTCg         >  4:361118/1‑132 (MQ=255)
                                                                                                                                gAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGCGTTGTCGa        >  7:3760/1‑133 (MQ=255)
                                                                                                                                   cTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGCGTTGTCGATAGAcc  >  8:353251/1‑136 (MQ=255)
                                                                                                                                    tGGCCTATGAAGAGGGGCGGGCCGGCGGGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGCGTTGTCGATAGAcc  <  4:98428/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
CGTGCTCGACCTCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTA‑CCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGCGTTGTCGATAGACC  >  KT2440/1508523‑1508788

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