Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 648,716 A→T V61E (GTG→GAG)  mrdB ← cell wall shape‑determining protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606648,7160AT97.5% 150.1 / ‑3.2 40V61E (GTG→GAG) mrdBcell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (22/17);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

CAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCG  >  REL606/648643‑648787
                                                                         |                                                                       
cAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCacctcc                                                                       >  1:1565667/1‑76 (MQ=255)
cAAAATAAGACAGAGGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCacctcc                                                                       <  2:1288915/76‑1 (MQ=255)
 aaaaTAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCAcctcca                                                                      >  2:1023002/1‑76 (MQ=255)
     taataCAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCatgat                                                                  >  1:97763/1‑76 (MQ=255)
      aaaaCAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCatgatg                                                                 >  2:46641/4‑76 (MQ=255)
        taCAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGAc                                                               >  1:1709468/1‑76 (MQ=255)
        aaCAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGAc                                                               <  2:1067344/75‑1 (MQ=255)
          cAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCAACTCCATGATGACCa                                                             <  1:988909/76‑1 (MQ=255)
              tgatATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGAcc                                                         <  1:1363458/76‑1 (MQ=255)
                atatAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGAACCa                                                       >  2:358561/1‑76 (MQ=255)
                   taGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATcg                                                    <  2:183474/76‑1 (MQ=255)
                   taGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATcg                                                    >  1:341581/1‑76 (MQ=255)
                    agagATAGGGGGCTCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGAGCCATcgc                                                   <  1:1179806/76‑1 (MQ=255)
                       gaTAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGAt                                                >  1:827286/1‑76 (MQ=255)
                        aTAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGAtt                                               >  1:184587/1‑76 (MQ=255)
                        aTAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGAtt                                               >  1:87513/1‑76 (MQ=255)
                        aTAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGAtt                                               >  1:900476/1‑76 (MQ=255)
                            ggggCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATAGCGATTTGGc                                           >  2:1264606/1‑76 (MQ=255)
                              ggCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCg                                         <  1:1023002/76‑1 (MQ=255)
                                cccAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGAt                                       >  1:288905/1‑76 (MQ=255)
                                 ccAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGAtt                                      >  2:1444849/1‑76 (MQ=255)
                                 ccAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGAtt                                      >  1:1549559/1‑76 (MQ=255)
                                   aGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGAtttt                                    <  2:226305/76‑1 (MQ=255)
                                   aGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGAtttt                                    >  2:1316528/1‑76 (MQ=255)
                                    gCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTa                                   >  1:188304/1‑76 (MQ=255)
                                      cTTCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACg                                 <  2:520929/76‑1 (MQ=255)
                                       ttCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGc                                >  1:1400113/1‑76 (MQ=255)
                                        tCATAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCt                               >  2:1149667/1‑76 (MQ=255)
                                          aTAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTcc                             >  1:240817/1‑76 (MQ=255)
                                          aTAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTcc                             <  1:1776870/76‑1 (MQ=255)
                                          aTAAACGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTcc                             >  2:1462550/1‑76 (MQ=255)
                                            aaaCGCGTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCcat                           >  2:1453386/1‑76 (MQ=255)
                                                  gTGGAGGAATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGcc                     >  2:1396966/1‑76 (MQ=255)
                                                       ggaATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAatat                <  1:413914/76‑1 (MQ=255)
                                                       ggaATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAatat                <  1:75544/76‑1 (MQ=255)
                                                       ggaATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAatat                <  1:617883/76‑1 (MQ=255)
                                                        gaATTTGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATc               <  1:618253/76‑1 (MQ=255)
                                                             tGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGAc          <  2:1311868/76‑1 (MQ=255)
                                                             tGCGCCATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGAc          >  2:251739/1‑76 (MQ=255)
                                                                 ccATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCt      <  1:774554/76‑1 (MQ=255)
                                                                 ccATCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCt      >  2:352981/1‑76 (MQ=255)
                                                                   aTCACGTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTgg    <  2:1806274/76‑1 (MQ=255)
                                                                    tCACCTCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGgc   <  2:123178/76‑1 (MQ=255)
                                                                     cacctcCATGATGNCCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGgcg  >  1:1241199/1‑76 (MQ=255)
                                                                     cacctcCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGgcg  >  2:397041/1‑76 (MQ=255)
                                                                     cacctcCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGgcg  >  2:865337/1‑76 (MQ=255)
                                                                         |                                                                       
CAAAATAATACAGATGATATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCG  >  REL606/648643‑648787

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