Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,915,303 A→G 48.5% F410L (TTT→CTT)  prc ← carboxy‑terminal protease for penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,915,3030AG48.5% ‑5.3 / 46.2 33F410L (TTT→CTT) prccarboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base A (4/13);  new base G (8/8);  total (12/21)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

ATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGAT  >  W3110S.gb/1915238‑1915366
                                                                 |                                                               
aTCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGAc                                                            >  1:2121064/1‑71 (MQ=255)
   gCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAGCAGACCg                                                          >  1:3318240/1‑70 (MQ=255)
           cACGATCCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGa                                                  <  1:1858097/70‑1 (MQ=255)
           cACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGa                                                  <  1:1524324/70‑1 (MQ=255)
             cGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATAc                                               <  1:3727707/71‑1 (MQ=255)
             cGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATAc                                               <  1:3202854/71‑1 (MQ=255)
                     gCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAGCAGACCGGAGAGCGATACGGCTTCAg                                       >  1:3902526/1‑71 (MQ=255)
                       cGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAGCAGACCGGAGAGCGATACGGCTTCAGtt                                     <  1:1184017/71‑1 (MQ=255)
                       cGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAGCAGACCGGAGAGCGATACGGCTTCAGtt                                     <  1:4046265/71‑1 (MQ=255)
                       cGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAGCAGACCGGAGAGCGATACGGCTTCAGtt                                     <  1:3903202/71‑1 (MQ=255)
                         ttgttATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGTCTTCAGt                                      <  1:968468/68‑1 (MQ=255)
                         ttgttATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGTCTTCAGt                                      <  1:963610/68‑1 (MQ=255)
                          tgttATCGCGGACCTGAACAATGGGACCCGCAGGAATAAGCAGACCGGAGAGCGATACGGCTTCAGTTAAc                                  >  1:211207/1‑71 (MQ=255)
                              aTCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGccc                              >  1:816808/1‑71 (MQ=255)
                                     aCCTGAACAATGGGACCCGCAGGAATAAGCAGACCGGAGAGCGa                                                  <  1:950953/44‑1 (MQ=255)
                                     aCCTGAACAATGGGACCCGCAGGAATAAGCAGACCGGAGAGCGa                                                  <  1:599123/44‑1 (MQ=255)
                                      ccTGAACAATGGGACCCGCAGGAATAAACAGACCGTAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCa                      <  1:2051019/71‑1 (MQ=255)
                                      ccTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCa                      <  1:1272802/71‑1 (MQ=255)
                                         gAACAATGGGACCCGCAGGAATAAGCAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTg                   >  1:3554/1‑71 (MQ=255)
                                          aaCAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACTGCCAtt                    <  1:527478/69‑1 (MQ=255)
                                          aaCAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGc                  <  1:1997421/71‑1 (MQ=255)
                                          aaCAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGc                  <  1:1577391/71‑1 (MQ=255)
                                            cAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAg                                       >  1:2044046/1‑48 (MQ=255)
                                             aaTGGGACCCGCAGGAATAAGCAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGc                        <  1:2613303/62‑1 (MQ=255)
                                             aaTGGGACCCGCAGGAATAAGCAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGc                        <  1:1603366/62‑1 (MQ=255)
                                               tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACg              >  1:3900420/1‑70 (MQ=255)
                                                       gCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCgat     <  1:1154940/71‑1 (MQ=255)
                                                       gCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCgat     <  1:2851431/71‑1 (MQ=255)
                                                          ggAATAAGCAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCgatgat  <  1:1900517/71‑1 (MQ=255)
                                                             aTAAGCAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCgatga   >  1:55840/1‑67 (MQ=255)
                                                             aTAAGCAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCgatga   >  1:664609/1‑67 (MQ=255)
                                                             aTAAGCAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCgatga   >  1:1761109/1‑67 (MQ=255)
                                                             aTAAGCAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCgatga   >  1:2309445/1‑67 (MQ=255)
                                                                 |                                                               
ATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGAT  >  W3110S.gb/1915238‑1915366

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