Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,716,832 A→C 6.7% intergenic (+90/+108) 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,8320AC 6.7% 193.3 / 9.5 75intergenic (+90/+108)yphH/glyApredicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (29/41);  new base C (2/3);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

CTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCA  >  minE/1716764‑1716901
                                                                    |                                                                      
cTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACg                                                                      >  1:306367/1‑71 (MQ=255)
   cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:1661947/70‑1 (MQ=255)
   cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:3790146/70‑1 (MQ=255)
   cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:3585983/70‑1 (MQ=255)
   cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:1290254/70‑1 (MQ=255)
    ccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGgag                                                                   <  1:3879586/70‑1 (MQ=255)
     cGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:3261819/68‑1 (MQ=255)
     cGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑At                                                                <  1:1044909/71‑1 (MQ=255)
            agCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAata                                                         >  1:2977423/1‑71 (MQ=255)
             gCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:4078590/60‑1 (MQ=255)
              ccTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAatata                                                       >  1:3051359/1‑71 (MQ=255)
               cTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTaca                                                                        >  1:3955819/1‑54 (MQ=255)
                ttGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:2090217/57‑1 (MQ=255)
                 tgatgaTCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATAtt                                                     >  1:1846672/1‑70 (MQ=255)
                  gatgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTaca                                                                        >  1:2982090/1‑51 (MQ=255)
                  gatgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTg                                                    <  1:3964185/70‑1 (MQ=255)
                   atgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCaa                                                           >  1:3674266/1‑62 (MQ=255)
                   atgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCaa                                                           >  1:1685762/1‑62 (MQ=255)
                    tgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAAt                                                 >  1:1039740/1‑71 (MQ=255)
                     gatCATAAACGAAAAATTGCCTGATGCGCTGCGCTTATCAGGCCTACACGgag                                                                   <  1:301944/53‑1 (MQ=255)
                       tCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTg                                              <  1:1006996/71‑1 (MQ=255)
                        cATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTg                                              >  1:2701873/1‑70 (MQ=255)
                        cATAAAAGAAAAATTGCCTGATGCTCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTg                                              >  1:1212670/1‑70 (MQ=255)
                         aTAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACTga                                                                    <  1:3571552/48‑1 (MQ=255)
                         aTAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAata                                                         >  1:1481500/1‑58 (MQ=255)
                          tAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAatata                                                       <  1:3940519/59‑1 (MQ=255)
                           aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCa                                            >  1:2118127/1‑69 (MQ=255)
                           aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCa                                            >  1:2461888/1‑69 (MQ=255)
                            aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGgag‑a                                                                 >  1:194831/1‑47 (MQ=255)
                            aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAg                                         >  1:3265183/1‑71 (MQ=255)
                            aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAg                                         >  1:1661596/1‑71 (MQ=255)
                             aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGt                                        <  1:1419646/71‑1 (MQ=255)
                             aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGt                                        <  1:1174803/71‑1 (MQ=255)
                             aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGt                                        <  1:2285852/71‑1 (MQ=255)
                               gAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGtt                                       <  1:1172683/70‑1 (MQ=255)
                                 caaaTTGCCTGATGCGCTACGCTTATCGGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTg                                    <  1:585976/70‑1 (MQ=11)
                                 caaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTg                                    <  1:1180270/70‑1 (MQ=11)
                                 caaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTg                                    <  1:2980459/70‑1 (MQ=11)
                                 aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg                                    <  1:2501861/71‑1 (MQ=255)
                                 aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg                                    <  1:2767673/71‑1 (MQ=255)
                                 aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg                                    <  1:3799106/71‑1 (MQ=255)
                                  aaaTTGCCTGATGCGCTGCGCTTATCAGGCCTACAAGTTGAATTGCAATTTATTGAATTTGCACATTTTTg                                    <  1:1206041/71‑1 (MQ=2)
                                  aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTTCAATATATTGAATTTGCAAAGTTTTGt                                   <  1:3970829/71‑1 (MQ=255)
                                  aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg                                    >  1:1864276/1‑70 (MQ=255)
                                  aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGt                                   <  1:811366/71‑1 (MQ=255)
                                       gCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGcc                              >  1:911080/1‑71 (MQ=255)
                                        ccTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCg                             <  1:61138/71‑1 (MQ=255)
                                         cTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCgg                            <  1:1305034/71‑1 (MQ=255)
                                           gATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                           >  1:1815686/1‑70 (MQ=255)
                                           gATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                           >  1:3293141/1‑70 (MQ=255)
                                           gATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGa                           >  1:834745/1‑70 (MQ=11)
                                            aTGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                         <  1:1412157/71‑1 (MQ=255)
                                            aTGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGAt                          <  1:1877370/70‑1 (MQ=11)
                                                 cTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCg                    <  1:346244/71‑1 (MQ=255)
                                                   aCGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCaaa                                          >  1:1564827/1‑47 (MQ=255)
                                                   aCGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGtttt                                     >  1:1861805/1‑52 (MQ=255)
                                                   aCGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                  >  1:3359635/1‑71 (MQ=255)
                                                   aCGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                  >  1:2876386/1‑71 (MQ=255)
                                                    tgCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTg                                    <  1:3181163/51‑1 (MQ=11)
                                                    tgCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTg                                    <  1:1662940/51‑1 (MQ=11)
                                                    cGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTc                 <  1:843926/71‑1 (MQ=255)
                                                      cTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGttt                                      >  1:2490392/1‑48 (MQ=255)
                                                      cTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAc               >  1:2571736/1‑71 (MQ=255)
                                                        tATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCact                                          >  1:3949887/1‑40 (MQ=11)
                                                         aTCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgc             <  1:2465711/70‑1 (MQ=255)
                                                         aTCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTcgcg              <  1:2009239/69‑1 (MQ=2)
                                                            aGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgccg           >  1:181541/1‑69 (MQ=255)
                                                              gCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc       <  1:2128009/71‑1 (MQ=255)
                                                              gCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc       >  1:528223/1‑71 (MQ=255)
                                                                cTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGttt                                      <  1:1282849/38‑1 (MQ=255)
                                                                cTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                  >  1:939185/1‑58 (MQ=255)
                                                                cTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg     >  1:1275493/1‑71 (MQ=255)
                                                                 tACACGGAG‑ATTGCAATATATTGAATTTGCAAAGttt                                      <  1:3633273/37‑1 (MQ=255)
                                                                 tACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCgg                            <  1:1101214/47‑1 (MQ=255)
                                                                 tACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCg                             <  1:3525913/46‑1 (MQ=255)
                                                                 tACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                          <  1:1665877/49‑1 (MQ=255)
                                                                 tACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                       <  1:798626/52‑1 (MQ=255)
                                                                  acacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                          <  1:575918/48‑1 (MQ=255)
                                                                  acacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg    >  1:2866669/1‑70 (MQ=255)
                                                                   cacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa                >  1:2585801/1‑57 (MQ=255)
                                                                   cacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgccgc          >  1:3393738/1‑63 (MQ=255)
                                                                   cacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg    <  1:307968/69‑1 (MQ=255)
                                                                   cacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa  <  1:413119/71‑1 (MQ=255)
                                                                   cacGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa  >  1:829181/1‑71 (MQ=255)
                                                                    |                                                                      
CTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCA  >  minE/1716764‑1716901

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: ‑