Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,312,395 A→G 48.9% S341G (AGC→GGC)  hisD → bifunctional histidinal dehydrogenase and histidinol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,312,3950AG48.9% ‑5.9 / 39.5 45S341G (AGC→GGC) hisDbifunctional histidinal dehydrogenase and histidinol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (11/12);  new base G (10/12);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGACTGGTCACCGGAATCGGCAGGTGATTACGCCTCC  >  minE/1312331‑1312463
                                                                |                                                                    
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCAt                                                                  >  1:311046/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCAt                                                                  >  1:1219115/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCAt                                                                  >  1:3004423/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCAt                                                                  >  1:3094935/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCAt                                                                  >  1:1383171/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCAt                                                                  >  1:1424219/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCAt                                                                  >  1:3159761/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCAt                                                                  >  1:1588603/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCAt                                                                  >  1:2500662/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCAt                                                                  >  1:1711087/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCAt                                                                  >  1:261642/1‑69 (MQ=255)
ctcCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCAGTGAACTGGTCGATAGCAt                                                                  >  1:1701413/1‑69 (MQ=255)
     aCCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATc                                                                 <  1:465296/65‑1 (MQ=255)
         gTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:3253324/59‑1 (MQ=255)
         gTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:2474099/59‑1 (MQ=255)
         gTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:2621967/59‑1 (MQ=255)
         gTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:2099685/59‑1 (MQ=255)
         gTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:3085517/59‑1 (MQ=255)
         gTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:2129595/59‑1 (MQ=255)
         gTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:1921264/59‑1 (MQ=255)
         gTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:474404/59‑1 (MQ=255)
         gTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:1522628/59‑1 (MQ=255)
         gTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:109831/59‑1 (MQ=255)
           aCGGCCCGGAGGACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:818062/57‑1 (MQ=255)
              gCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGCGAACTGGTCGATGGCa                                                                   <  1:3259077/54‑1 (MQ=39)
                 cGGAGCACCTGATCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCgg                                                 <  1:2191234/69‑1 (MQ=255)
                     gCACCTGATCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGt                                                <  1:711575/66‑1 (MQ=255)
                     gCACCTGATCATTCAGACCCGCAACGCCCGTAAACTGGTCGATAGCATCACCAGCGCCGGTTCGGt                                                <  1:706480/66‑1 (MQ=255)
                      cACCTGATCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGt                                                <  1:2707932/65‑1 (MQ=255)
                           gATCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCgg                                                 <  1:1922651/59‑1 (MQ=255)
                            aTCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGt                                      >  1:944814/1‑69 (MQ=255)
                            aTCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGt                                      >  1:823329/1‑69 (MQ=255)
                            aTCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGt                                      >  1:1235110/1‑69 (MQ=255)
                               aTTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGAc                                   >  1:2954298/1‑69 (MQ=255)
                               aTTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGAc                                   >  1:2179190/1‑69 (MQ=255)
                               aTTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGAc                                   >  1:3349730/1‑69 (MQ=255)
                               aTTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGAc                                   >  1:1366784/1‑69 (MQ=255)
                               aTTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGAc                                   >  1:1187069/1‑69 (MQ=255)
                                          cAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGACTGGTCACCGGa                        <  1:2692190/69‑1 (MQ=255)
                                          cAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGACTGGTCACCGGa                        <  1:2138563/69‑1 (MQ=255)
                                            aCGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGACTGGTCACCGGAAt                      >  1:2777495/1‑69 (MQ=255)
                                                         ggTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGt                                      <  1:472437/40‑1 (MQ=255)
                                                         ggTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGt                                      <  1:1645654/40‑1 (MQ=255)
                                                         ggTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGt                                      <  1:1307343/40‑1 (MQ=255)
                                                         ggTCGATAGCATCACCAGCGCCGGTTCGGTATTTCCTGGt                                      <  1:905311/40‑1 (MQ=255)
                                                                aGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGACTGGTCACCGGAATCGGCAGGt              <  1:1518099/57‑1 (MQ=255)
                                                                aGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGACTGGTCACCGGAATCGGCAGGTGATTACGCCTcc  >  1:2846417/1‑69 (MQ=255)
                                                                |                                                                    
CTCCAACCAGTACGGCCCGGAGCACCTGATCATTCAGACCCGCAACGCCCGTGAACTGGTCGATAGCATCACCAGCGCCGGTTCGGTATTTCTTGGTGACTGGTCACCGGAATCGGCAGGTGATTACGCCTCC  >  minE/1312331‑1312463

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: