Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,314,264 T→C 50.6% R173R (CGT→CGC hisB → fused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,314,2640TC50.6% ‑2.2 / 71.5 79R173R (CGT→CGChisBfused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (29/10);  new base C (21/19);  total (50/29)
Fisher's exact test for biased strand distribution p-value = 6.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

CAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAAC  >  minE/1314212‑1314330
                                                    |                                                                  
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:2531587/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:1767672/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:1808150/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:1850906/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:2022661/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:2200629/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:2372258/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:2428536/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:2457868/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:1176008/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:25853/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:2756915/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:3071039/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:320169/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:405993/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:675360/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:842111/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:993568/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:1004565/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:1294335/3‑68 (MQ=255)
cgaTGATTGGTGAGCAACTCACTAAACGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACg                                                     >  1:1351958/3‑68 (MQ=255)
    gATTGGCGAGCAACTCACCAGACGTGACCGTTATGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                                                 >  1:2625292/1‑68 (MQ=255)
    gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                                                 >  1:844075/1‑68 (MQ=255)
    gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                                                 >  1:506822/1‑68 (MQ=255)
    gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                                                 >  1:467611/1‑68 (MQ=255)
    gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                                                 >  1:1717406/1‑68 (MQ=255)
    gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                                                 >  1:1546485/1‑68 (MQ=255)
    gATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGa                                                 >  1:3278024/1‑68 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGTTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:2157673/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:2827078/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:2890087/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:285864/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:1401765/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:3119884/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:3427525/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:36162/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:384040/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:392819/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:463524/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:836983/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:843167/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:2783626/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:2248628/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:1695059/1‑69 (MQ=255)
     aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAtt                                               >  1:1514384/1‑69 (MQ=255)
       tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                                             >  1:2177280/1‑69 (MQ=255)
       tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                                             >  1:1550899/1‑69 (MQ=255)
       tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                                             >  1:1570672/1‑69 (MQ=255)
       tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                                             >  1:645457/1‑69 (MQ=255)
            aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc                                        <  1:2142160/69‑1 (MQ=255)
            aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc                                        <  1:3198566/69‑1 (MQ=255)
            aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc                                        <  1:220111/69‑1 (MQ=255)
            aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc                                        <  1:899019/69‑1 (MQ=255)
            aGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTCACGTTc                                        <  1:435799/69‑1 (MQ=255)
                   cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATGCCAAAGAGACGCAGATTGACGTTCAGGtgtg                                 <  1:1304822/69‑1 (MQ=255)
                   cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg                                 <  1:2886370/69‑1 (MQ=255)
                   cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg                                 <  1:1397864/69‑1 (MQ=255)
                   cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg                                 <  1:2858445/69‑1 (MQ=255)
                   cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg                                 <  1:2246130/69‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCTCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:1715190/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAATAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:978431/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:2865435/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:2791658/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:1993443/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:2000525/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:2014797/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:2255749/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:596153/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:1575179/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:2946927/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:3038616/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:1542223/66‑1 (MQ=255)
                      taaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:3122365/66‑1 (MQ=255)
                       aaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:1460387/66‑1 (MQ=255)
                       aaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:914324/66‑1 (MQ=255)
                       aaaCGAGACCGTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:429962/66‑1 (MQ=255)
                                 gTTACGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:2562875/58‑1 (MQ=255)
                                    aCGCCCATGTAGTGCGCAACACCAAAGAGACGCAAATTGACGTCCAGGTGTGGCt                              <  1:808087/55‑1 (MQ=255)
                                                  cGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAac  >  1:1407849/1‑69 (MQ=255)
                                                    |                                                                  
CAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAAC  >  minE/1314212‑1314330

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: