Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,175,232 (T)6→7 coding (877/933 nt) pbpG ← D‑alanyl‑D‑alanine endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,175,2261.T97.1% 130.5 / ‑3.4 35coding (883/933 nt)pbpGD‑alanyl‑D‑alanine endopeptidase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (13/21);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCG  >  REL606/2175152‑2175295
                                                                                |                                                                     
cAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCCTTTTGGCGGGTTTTTGt                                                                      >  2:980414/1‑76 (MQ=255)
cAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCGGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGt                                                                      <  1:558840/76‑1 (MQ=255)
    gTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGttttt                                                                  <  2:1728889/76‑1 (MQ=255)
      aTCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGttttttt                                                                <  2:206701/76‑2 (MQ=255)
         gCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCCCCCCTTTGGGGGGGttttttttttttttt                                                              >  1:1646963/1‑73 (MQ=255)
           cATCCGATTTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGc                                                            <  1:1434882/76‑1 (MQ=255)
              ccGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCa                                                         <  2:520747/76‑1 (MQ=255)
              ccGAATTCACTTAATCGTTCTGTGCCGCCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCa                                                         <  1:782233/76‑1 (MQ=255)
                aatttCCCTTATCCGTGCTGTGCCGTCGGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAAc                                                       <  1:575651/73‑1 (MQ=255)
                 aaTTCACTTAGTCGTTCTGTGCCTTCGCCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACg                                                      <  2:1685760/76‑1 (MQ=255)
                    tCACTTAATCTTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACgctg                                                   <  2:629987/76‑1 (MQ=255)
                    tCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCCTTTGGGGGGCTTTTTTTTTTTTTTTTGCTCAACgctg                                                   >  1:1894933/1‑76 (MQ=255)
                       cTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGctgctg                                                <  2:2055809/76‑1 (MQ=255)
                          aaTCTTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCg                                             <  1:733655/76‑1 (MQ=255)
                                  tgtgCCGTCTGCCCCGCCGCCCCCCCTTTGGGGGCTTTTTTTTTTTTTTTTGCTCAACGCTGCTGCCGGCCCAGGc                                     >  2:1611257/1‑76 (MQ=255)
                                  tgtgCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGc                                     >  2:385965/1‑76 (MQ=255)
                                   gggCCGTCGGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCa                                    <  2:1609173/74‑1 (MQ=255)
                                      ccGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGGATCa                                 >  2:1805663/1‑76 (MQ=255)
                                      ccGTCTGCCCCGCCGCCCCCCTTTGGGGGGGTTTTTTTTTTTTTTTAGCTCAACGCTGCTGCCGGCACAGGCATCa                                 >  2:1871884/1‑76 (MQ=255)
                                      ccGTCTGCCCCGCCGCCCCCCCTTTGGCGGCCTTTTTTTTTTTTTTTTTTCAACCCTTCTGCCCGCCCACGCCTCa                                 >  2:317158/1‑76 (MQ=255)
                                      ccGTCTGCCCCGCCCCCCCCATTTTGGCGGGTTTTTGGTTTTTTTTAGCTCAAAGCTGCTGCCTGCACAGGCATCa                                 >  2:183861/1‑76 (MQ=255)
                                          ccgCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACttt                             <  1:454349/74‑1 (MQ=255)
                                           tcccccGCCGCCGCCATTGGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTa                            <  1:940951/74‑1 (MQ=255)
                                             ccccGCCGCCGCCCTTTTGGCGGCTTTTTGTTTTTTTATAATTCTACCAGGCTGCAGGCACAGGGAACACCTTAcc                          >  2:768045/1‑76 (MQ=255)
                                              cccGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCg                         >  2:500802/1‑76 (MQ=255)
                                               cccccgccgccATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCgg                        <  1:1322688/73‑1 (MQ=255)
                                                cgccgccgcCATTTTGGCGGCTTTTTGTTTTTTTATAGCTCAATGCTGCAGCCGGCCCAGGCATCACTTTACCGGt                       >  2:664863/1‑76 (MQ=255)
                                                cgccgccgcCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGt                       >  2:746078/1‑76 (MQ=255)
                                                cgccgccgcCATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGAATAACTTTACCGGt                       >  2:32293/1‑76 (MQ=255)
                                                 gccgccgccATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGtt                      >  2:1191074/1‑76 (MQ=255)
                                                 gccgccgccATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGtt                      >  2:1236041/1‑76 (MQ=255)
                                                  ccgccgccATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCCCAGGCATCACTTTACCGGttt                     >  2:1307950/1‑76 (MQ=255)
                                                   cgccgccATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTc                    <  2:1289028/76‑1 (MQ=255)
                                                    gccgccATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCa                   <  2:274567/76‑1 (MQ=255)
                                                         ccATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATcc               >  2:1379177/1‑76 (MQ=255)
                                                         ccATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATcc               <  1:2004917/76‑1 (MQ=255)
                                                          cATTTGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCa              <  2:1093466/76‑1 (MQ=255)
                                                            tttGGGCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAgg            <  1:1436315/76‑1 (MQ=255)
                                                                gggCGGCTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTAc         >  2:1075354/1‑76 (MQ=255)
                                                                 ggcggcTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACg        <  1:3422/76‑1 (MQ=255)
                                                                 ggcggcTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACg        <  2:461207/76‑1 (MQ=255)
                                                                   cggcTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCa      >  2:654474/1‑76 (MQ=255)
                                                                   cggcTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCa      >  1:1489303/1‑76 (MQ=255)
                                                                    ggcTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAg     <  1:742522/76‑1 (MQ=255)
                                                                     gcTTTTTGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAgg    <  1:565863/76‑1 (MQ=255)
                                                                         tttttGTTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAggcg  >  2:779480/1‑76 (MQ=255)
                                                                                |                                                                     
CAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCG  >  REL606/2175152‑2175295

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: