Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,716,884 T→C 16.0% intergenic (+142/+56) yphH → / ← glyA predicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,716,8840TC16.0% 116.4 / 16.0 57intergenic (+142/+56)yphH/glyApredicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (26/21);  new base C (6/3);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

CTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACC  >  minE/1716813‑1716952
                                                                       |                                                                    
cgatGCTTATCAGGCCTACAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGt                                                                      <  1:2714889‑M2/67‑1 (MQ=255)
  atGCTTATCAGGCCTACCAGAAGATAGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGtt                                                                     <  1:4203362‑M2/69‑1 (MQ=255)
   cGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                                                                     <  1:5436240/70‑1 (MQ=255)
     cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa                                                                   >  1:2513460/1‑70 (MQ=255)
     cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa                                                                   >  1:4071410/1‑70 (MQ=255)
     cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAc                                                                  >  1:3248722/1‑71 (MQ=255)
     cTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATATGCGATGTTTTGTAGGACTGATAAGGCGCTcg                                                                   >  1:5395870/1‑69 (MQ=2)
        aTCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTcgcg                                                                 <  1:3798215/69‑1 (MQ=2)
          cAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGcgccg                                                              >  1:934475/1‑70 (MQ=2)
                  cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                       <  1:5050166/69‑1 (MQ=255)
                   acGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                    >  1:5566762/1‑71 (MQ=255)
                    agaagaTTGCAATATATTGAATTTGCACTGTTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATg                                                   >  1:1023432/4‑70 (MQ=11)
                       agaTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaac                                                <  1:4980081/71‑1 (MQ=11)
                       agaTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaac                                                <  1:1869594/71‑1 (MQ=11)
                       agaTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaac                                                <  1:5502856/71‑1 (MQ=11)
                       agaTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaac                                                <  1:1028001/71‑1 (MQ=11)
                        caTTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATCaact                                               <  1:5329792/70‑2 (MQ=11)
                        caTTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATCaac                                                >  1:3731757/2‑70 (MQ=11)
                         aTTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATcaacaa                                              <  1:3874039/71‑1 (MQ=11)
                           tGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                             >  1:5574960/1‑70 (MQ=255)
                           tGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                             >  1:1650007/1‑70 (MQ=255)
                             cAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAg                                          >  1:3828834/1‑71 (MQ=255)
                              aaTATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGc                                         <  1:3284311/71‑1 (MQ=255)
                               atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca                                        <  1:3494598/71‑1 (MQ=255)
                               atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca                                        <  1:5076453/71‑1 (MQ=255)
                               atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca                                        <  1:4537113/71‑1 (MQ=255)
                                tctatTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATcaacaa                                              >  1:882821/3‑64 (MQ=11)
                                tctatTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATcaacaa                                              >  1:2837279/3‑64 (MQ=11)
                                  tatTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATCaac                                                >  1:1107182/1‑60 (MQ=11)
                                  tatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca                                      <  1:4737930/70‑1 (MQ=255)
                                  tatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACAt                                     <  1:575186/71‑1 (MQ=255)
                                  ctttGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCgg                                                       >  1:2539797/3‑52 (MQ=0)
                                  ctttGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCgg                                                       >  1:1595834/3‑52 (MQ=0)
                                  ctttGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCgg                                                       >  1:1909925/3‑52 (MQ=0)
                                  ctttGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGc                                                      <  1:2514547/51‑1 (MQ=0)
                                      gAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGa                                  >  1:3443282/1‑70 (MQ=255)
                                       gaTTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCaa                                                 >  1:3873219/2‑53 (MQ=11)
                                           tGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCa                             >  1:5128176/1‑70 (MQ=255)
                                                aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAc                       >  1:1399563/1‑71 (MQ=255)
                                                aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAc                       >  1:1091837/1‑71 (MQ=255)
                                                    ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGt                    <  1:4691661/70‑1 (MQ=255)
                                                    ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGt                    <  1:1857053/70‑1 (MQ=255)
                                                    ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGt                    <  1:1625196/70‑1 (MQ=255)
                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCt              <  1:810276/71‑1 (MQ=255)
                                                              gATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGc         >  1:1768273/1‑71 (MQ=255)
                                                               aTAAGGCGTTCACGCCGCATCCTGCATGAACAACGAGCACATTGACAGCAAATCAc                       >  1:3559667/1‑56 (MQ=255)
                                                               aTAAGGCGTTCACGCCGCATCCGGCTTGAACAACGAGCACATTGACAGCAAATCAc                       >  1:2560486/1‑56 (MQ=255)
                                                               aTAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAc                       >  1:2227212/1‑56 (MQ=255)
                                                               aTAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAc                       >  1:3092956/1‑56 (MQ=255)
                                                               aTAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAc                       >  1:1876317/1‑56 (MQ=255)
                                                                tAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCtt             >  1:3746483/1‑65 (MQ=255)
                                                                  aGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGt                    <  1:3070448/56‑1 (MQ=255)
                                                                  aGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGt                    <  1:99794/56‑1 (MQ=255)
                                                                  aGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGt                    <  1:5170061/56‑1 (MQ=255)
                                                                  aGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGt                    <  1:175512/56‑1 (MQ=255)
                                                                  aGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTa      >  1:4653903/1‑70 (MQ=255)
                                                                    gCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTa            >  1:5472104/1‑62 (MQ=255)
                                                                    gCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTa            >  1:5010198/1‑62 (MQ=255)
                                                                    gCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAAc   <  1:3648637/71‑1 (MQ=255)
                                                                    gCGGTCACGCCGCATCCGGCATGAACAACGAGCACATTg                                   <  1:3785733/39‑1 (MQ=255)
                                                                       ttCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAAcc  <  1:3534560/69‑1 (MQ=255)
                                                                       |                                                                    
CTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACC  >  minE/1716813‑1716952

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: