Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,916,448 C→T 11.1% R28H (CGT→CAT)  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,916,4480CT11.1% 131.9 / 7.2 45R28H (CGT→CAT) prccarboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base T (3/2);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GCGACGTTACGCGCTCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGC  >  W3110S.gb/1916380‑1916513
                                                                    |                                                                 
gcgACGTTACGCGCTCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCAcgc                                                                 <  1:4268299/71‑1 (MQ=255)
   aCGTTACGCGCTCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTg                                                              <  1:539828/71‑1 (MQ=255)
     gTTACGCGCTCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGa                                                             >  1:4268015/1‑70 (MQ=255)
      ttACGCGCTCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCAc                                                                   <  1:3082019/63‑1 (MQ=255)
        aCGCGCTCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATAtc                                                         >  1:2804523/1‑71 (MQ=255)
        aCGCGCTCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATAtc                                                         >  1:690994/1‑71 (MQ=255)
         cgcgCTCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATAtc                                                         <  1:114013/70‑1 (MQ=255)
             cTCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTa                                                    <  1:3431725/71‑1 (MQ=255)
              tCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTa                                                    >  1:2917702/1‑70 (MQ=255)
               cACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTa                                                    >  1:4694234/1‑69 (MQ=255)
                aCTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCATGCGTGATATCTTCTACAg                                                 >  1:3263003/1‑71 (MQ=255)
                aCTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCATGCGTGATATCTTCTACAg                                                 >  1:1347619/1‑71 (MQ=255)
                aCTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCATGCGTGATATCTTCTACAg                                                 >  1:3410728/1‑71 (MQ=255)
                 cTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGc                                                >  1:299918/1‑71 (MQ=255)
                    aCCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCg                                               >  1:1811919/1‑69 (MQ=255)
                    aCCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCg                                               <  1:2582263/69‑1 (MQ=255)
                    aCCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCg                                               >  1:2752898/1‑69 (MQ=255)
                    aCCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGa                                              <  1:1724857/70‑1 (MQ=255)
                      cGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTa                                                    >  1:1521145/1‑62 (MQ=255)
                             tgctgcGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGa                                                             >  1:257289/1‑46 (MQ=255)
                              gctgcGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTg                                       <  1:3562521/67‑1 (MQ=255)
                                 gcGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGatat                                                          >  1:4276467/1‑45 (MQ=255)
                                 gcGTCTCTTCCTTTAATACAGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCt                                        >  1:720354/1‑63 (MQ=255)
                                   gTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTAt                              <  1:705524/71‑1 (MQ=255)
                                     ctctTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTAt                              >  1:210797/1‑69 (MQ=255)
                                        ttCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGc                           <  1:2931003/69‑1 (MQ=255)
                                          ccTTTAATACCGGAATTTGATCAGCATGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCa                          <  1:82415/68‑1 (MQ=255)
                                          ccTTTAATACCGGAATTTGATCAGCATGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCa                          <  1:765251/68‑1 (MQ=255)
                                          ccTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAg                        <  1:3498809/70‑1 (MQ=255)
                                          ccTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAg                        <  1:2938355/70‑1 (MQ=255)
                                           cTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGc                       <  1:3369005/70‑1 (MQ=255)
                                           cTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGc                       <  1:2728246/70‑1 (MQ=255)
                                           cTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGc                       <  1:2722589/70‑1 (MQ=255)
                                              tAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTAt                              >  1:1917493/1‑60 (MQ=255)
                                              tAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGc                   >  1:2725926/1‑71 (MQ=255)
                                                   ccGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAg                                            >  1:580685/1‑41 (MQ=255)
                                                   ccGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCa                      >  1:326834/1‑63 (MQ=255)
                                                   ccGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCt              <  1:3541874/71‑1 (MQ=255)
                                                   ccGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCt              <  1:2422159/71‑1 (MQ=255)
                                                         tttGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGc                       <  1:3874581/56‑1 (MQ=255)
                                                          ttGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGc                       <  1:3815322/55‑1 (MQ=255)
                                                           tGATCAGCACGCGTGATGTCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCgg       >  1:1907795/1‑70 (MQ=255)
                                                            gATCAGCACGCGTGATATCTTCTACAGCGAAGGTCt                                        >  1:929541/1‑36 (MQ=255)
                                                              tCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTaa    <  1:1484893/70‑1 (MQ=255)
                                                              tCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAg   >  1:1160670/1‑71 (MQ=255)
                                                               cAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTAtt                             >  1:2135048/1‑44 (MQ=255)
                                                               cAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGc  >  1:649901/1‑71 (MQ=255)
                                                                    |                                                                 
GCGACGTTACGCGCTCACTTACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGC  >  W3110S.gb/1916380‑1916513

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: