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

CGTGCTCGACCTCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGCGTTGTCGATAGACC  >  KT2440/1508523‑1508788
                                                                                                                                    |                                                                                                                                      
cGTGCGCGACCTCAAAGTGCCTGGTGGGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGGGACTGCGCACAGGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGc                                                                                                                                     <  2:306790/135‑1 (MQ=255)
     tCGACCTCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACCTGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGGAGGGGTGGAGCGGGCCTAt                                                                                                                                 >  5:279913/1‑134 (MQ=255)
     tCGACCTCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATg                                                                                                                                >  8:445937/1‑135 (MQ=255)
       gACCTCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGGCGGGGGGGGGGTGGCCttt                                                                                                                                 >  5:99055/1‑130 (MQ=255)
           tCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGGAGGGGTGGAGCGGGGCGATGAAGAGt                                                                                                                          >  3:95407/1‑135 (MQ=255)
                gTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGGACGGGTGGGGGTGGGCTATTAATAGttggtt                                                                                                                     >  3:281994/1‑130 (MQ=255)
                     tGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGGAGGGGGGGGGGGGGGCGATGAAGAGTGGCTGGGCgg                                                                                                                >  1:117361/1‑135 (MQ=255)
                      ggTGTGTACGGACTGGAAGACTTCCGCCGCTGCTATCCGGCCGGTGATGTGCCTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGGATGGGGGGGGGGGGGCGATGAAGAGTGGCTGGCCGGcc                                                                                                              >  7:184844/1‑135 (MQ=255)
                                    ggAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGGTGGGGTGGGGCGGGGCGATGTAGCGTGGCTGGCCGGGGGGCCCCGccaccg                                                                                                 >  5:349210/1‑130 (MQ=255)
                                    ggAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGCAGGGGGGGAGCGGGGCGATGAAGACTGGCTGGCCGGCGGGCCCCGCCAgcg                                                                                                 >  4:369167/1‑135 (MQ=255)
                                    ggAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGAAGGGGTGGCGCTGGCCAATGCACACTGGCTGGCCCGCGGGCCCCGccaccc                                                                                                 >  4:345281/1‑130 (MQ=255)
                                    ggAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGAAGCGGTCGAGCAGGCCTATGAAGACTGGCTGGCCGGCCTGCCCGGCAAacg                                                                                                 >  6:104677/1‑135 (MQ=255)
                                     tattaaTTCCGCGTCTCATAGCGGGCGGGTGATGTGACTCCGACCATGGGGGGCTTTACCGACGGCCACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATCAAGAGTGGCTGGCCGGCGTGCCCGGAAAGCGcc                                                                                               <  8:96327/132‑1 (MQ=255)
                                      aataaTGCCGCCGCTTCAATCCGGCCGGTTCGGTCACGGCGAACAGGGGGGGGTTACCCCACCGCGACGACCACGGGCGTGAAGGGGTGGAGCTGGCCTATGAAGAGTGGCGGGCCGGCGTGCCCGGCAAGCGcc                                                                                               <  6:279913/132‑1 (MQ=255)
                                      aagaaTTCCGCCGTTTCTAACCGGCCGGTGATGGGCCTGGGCACAGGGGCGGCTTTCCCGACCGCGACGACCACGGTCGTGATGGGGTGGAGCTGGCCTATGACGAGGGGCTGGCCGGCGTGCCCGGCAAGCGcc                                                                                               <  6:122510/135‑1 (MQ=255)
                                         aaTTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGGAGGGGTGGCGCTGGGCGATGTAGAGTGGCTGGCCGGCGTGGCCGGCCAGCGCCCgg                                                                                            >  1:92985/1‑135 (MQ=255)
                                          tttCCGCCGCTTCTTCCGGGCCGGGGAGGTGACTGCGCACTTGGTGGGCTTTCCCGACCGCGACGACCACGGTGGTGATGGGGTGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTg                                                                                          <  7:445937/135‑1 (MQ=255)
                                              cgccgcTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGCAGGGGTGGGGCGGGCCTATGAAGAGTGGCTGGCCCGCGGGGCCGGCACGCGCCAGGTGAt                                                                                        >  2:342950/1‑134 (MQ=255)
                                                             gccggTGATGTGCCTGCGCACTTGGGCGGCTTACCCGACGTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCAATGAAGAGTGGCTGGCCGGGGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGt                                                                        <  1:83290/135‑1 (MQ=255)
                                                                      gttcgtGCGCACATGGTGGTTTTACCCTACGTGGACGCCCTCGGGCGGGGCGGGGGGGAGCTGGCCTATGAAGAGGGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGAt                                                                <  1:57690/129‑1 (MQ=255)
                                                                            gcgcACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTg                                                                                          >  8:323015/1‑102 (MQ=255)
                                                                            gcgcACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTg                                                                                          <  7:323015/102‑1 (MQ=255)
                                                                                        ggCTTTCCCGACGTCGAGGACCACGGTGGTTATGGGTGGGTGGTGGCGTTTGATGTGTGGGTGGCCGGGGTGCCGGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAAcc                                             <  2:241634/135‑1 (MQ=255)
                                                                                          cTTTACCGACCTCGACGACCACGGTCGTGCAGGGGTGGGGCGGGCCTATGAAGAGTGGCTGGCCCGCGTGCCCGGCCAGCGCCAGGTGAGCAAGGCCCGGCGCGGGCGGCTGATGAAGGGCATTGAGGTAAc                                              >  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)
                                                                                                                                    |                                                                                                                                      
CGTGCTCGACCTCAAAGTGCCTGGTGTGTACGGACTGGAAGAATTCCGCCGCTTCTATCCGGCCGGTGATGTGACTGCGCACATGGTCGGCTTTACCGACCTCGACGACCACGGTCGTGAAGGGGTGGAGCTGGCCTATGAAGAGTGGCTGGCCGGCGTGCCCGGCAAGCGCCAGGTGATCAAGGACCGGCGCGGTCGGCTGATCAAGGACATTCAGGTAACCAAGAACGCCAAGGCCGGGAAGACCTTGGCGTTGTCGATAGACC  >  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: