Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,643,819 C→T 12.8% G721G (GGG→GGA pbpC ← fused transglycosylase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,643,8190CT12.8% 100.5 / 5.9 39G721G (GGG→GGApbpCfused transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (15/19);  new base T (3/2);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

TCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTAAACGTTTAATAATCGCGCCATCGC  >  W3110S.gb/2643758‑2643889
                                                             |                                                                      
tCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACtt                                                               <  1:314595/71‑1 (MQ=255)
            cgcgTTCAGTTAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACtt                                                               <  1:2270161/59‑1 (MQ=255)
            cgcgTTCAGTTAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGc                                                   <  1:338432/71‑1 (MQ=255)
             gcgTTCAGTTAACGGTTCGCGATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCa                                                  <  1:1105888/71‑1 (MQ=255)
             gcgTTCAGTTAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGc                                                   >  1:4018977/1‑70 (MQ=255)
             gcgTTCAGTTAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCa                                                  <  1:678311/71‑1 (MQ=255)
                ttCAGTTAACGGTTCGCCGTTCAGAAACCACCAGCGTTCACCTGCTCCTCCACTTGa                                                             <  1:2913066/57‑1 (MQ=255)
                ttCAGTTAACGGTTCGCCGTTCAGAAACCACCAGCGTTCACCTGCTCCTCCACTTGa                                                             <  1:4410419/57‑1 (MQ=255)
                ttCAGTTAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAg                                               <  1:2068613/71‑1 (MQ=255)
                ttCAGTTAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAg                                               <  1:2543242/71‑1 (MQ=255)
                     ttAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCt                                          <  1:4198281/71‑1 (MQ=255)
                      tAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCaa                                                       >  1:1744620/1‑57 (MQ=255)
                       aaCGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGa                                                             >  1:3254629/1‑50 (MQ=255)
                       aaCGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGa                                                             >  1:4071139/1‑50 (MQ=255)
                       aaCGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTc                                        >  1:3140789/1‑71 (MQ=255)
                        aCGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTcc                                       <  1:775865/71‑1 (MQ=255)
                        aCGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTc                                        <  1:2972807/70‑1 (MQ=255)
                           gTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCt                                    <  1:1392754/71‑1 (MQ=255)
                               gCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTg                                                              <  1:3682498/41‑1 (MQ=255)
                                 cATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGccc                              >  1:3589316/1‑71 (MQ=255)
                                   ttCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCgg                            >  1:1660159/1‑71 (MQ=255)
                                    tCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCa                                                  <  1:1479226/48‑1 (MQ=255)
                                    tCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGc                                     >  1:109899/1‑61 (MQ=255)
                                     cAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGt                           >  1:3031961/1‑70 (MQ=255)
                                     cAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGt                           >  1:3994913/1‑70 (MQ=255)
                                     cAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGt                           >  1:3082431/1‑70 (MQ=255)
                                     cAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCAGCGCCCGGt                           >  1:3518856/1‑70 (MQ=255)
                                      aGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTa                          <  1:3466671/70‑1 (MQ=255)
                                           ccaccaGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTAAACGtt                    <  1:3793817/71‑1 (MQ=255)
                                             accaGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTAAACGTTta                  <  1:4154499/71‑1 (MQ=255)
                                               caGCGTTCACCTGCTCCTCCACTTGATTGCAACGGCaaa                                                >  1:223982/1‑39 (MQ=255)
                                               caGCGTTCACCTGCTCCTCCACTTGATTGCAACGGCaaa                                                >  1:1536388/1‑39 (MQ=255)
                                               caGCGTTCACCTGCTCCTCCACTTGATTGCAACGGCaaa                                                >  1:1056105/1‑39 (MQ=255)
                                                  cGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTAAACGTTtaataa              <  1:2097189/70‑1 (MQ=255)
                                                   gTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTAAACGTTTAATAATc            <  1:361605/71‑1 (MQ=255)
                                                     tCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGc                                           >  1:2277382/1‑38 (MQ=255)
                                                        ccTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTAAACGTTTAATAATCgcgc        >  1:1696422/1‑70 (MQ=255)
                                                         cTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTAAAc                       <  1:213349/54‑1 (MQ=255)
                                                           gCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTAAACGTTTAATAATCGCGcc       >  1:4172666/1‑68 (MQ=255)
                                                             cccTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTTAACGTTTAATAATCGCGCCATcgc  <  1:2649228/71‑1 (MQ=255)
                                                             cccTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTAAACGTTTAATAATCGCGCCATcgc  >  1:2092694/1‑71 (MQ=255)
                                                             |                                                                      
TCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCACCAGCGTTCACCTGCCCCTCCACTTGATTGCAACGGCAAAGTCGCTTCCGCTGCGCCCGGTAAACGTTTAATAATCGCGCCATCGC  >  W3110S.gb/2643758‑2643889

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