Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 152,487 G→A 52.9% H407H (CAC→CAT hemL ← glutamate‑1‑semialdehyde aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE152,4870GA52.9% 4.0 / 76.3 87H407H (CAC→CAThemLglutamate‑1‑semialdehyde aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (26/15);  new base A (16/30);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 1.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

CACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGCAT  >  minE/152428‑152552
                                                           |                                                                 
cacaACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACg                                                          >  1:1093476/1‑69 (MQ=255)
   aaCTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGa                                                        <  1:1753221/68‑1 (MQ=255)
      ttCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                                    <  1:3106048/69‑1 (MQ=255)
      ttCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                                    <  1:2476897/69‑1 (MQ=255)
                  cGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCa                                        <  1:2130588/69‑1 (MQ=255)
                  cGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCa                                        <  1:2587960/69‑1 (MQ=255)
                  cGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCa                                        <  1:2667335/69‑1 (MQ=255)
                  cGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCa                                        <  1:1251454/69‑1 (MQ=255)
                  cGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCa                                        <  1:1142076/69‑1 (MQ=255)
                     cGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGAACATAAAGCCCGCTTCaaa                                      <  1:191410/68‑1 (MQ=255)
                         cagcaTCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCt                                  >  1:2995113/1‑68 (MQ=255)
                         cagcaTCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCt                                  >  1:1172069/1‑68 (MQ=255)
                         cagcaTCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTg                                 >  1:1227757/1‑69 (MQ=255)
                         cagcaTCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTg                                 >  1:1192505/1‑69 (MQ=255)
                         cagcaTCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTg                                 >  1:523503/1‑69 (MQ=255)
                           gcaTCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGAc                               <  1:3431155/69‑1 (MQ=255)
                           gcaTCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGAc                               <  1:1690819/69‑1 (MQ=255)
                           gcaTCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTAAAACGCTGAc                               <  1:2567738/69‑1 (MQ=255)
                             aTCGATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:871498/68‑1 (MQ=255)
                                gATTGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCg                                             >  1:1571068/1‑50 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCt                                  >  1:2247520/1‑61 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCt                                  >  1:1450389/1‑61 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCTCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:2217592/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:3385216/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:2568637/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:3279690/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:3237146/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:2549335/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:3084089/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:2759100/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:2918491/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:2887105/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:531264/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:691101/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:1030578/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:873318/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:1215640/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:86420/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:757507/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:1482436/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:752140/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:725572/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:2296763/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:1666724/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:605561/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:1744332/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:1965413/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:2022740/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:2061270/65‑1 (MQ=255)
                                gATGGTGTTATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:2068920/65‑1 (MQ=255)
                                  tGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGcc                         <  1:974529/68‑1 (MQ=255)
                                     tgtTATTGATATCATCCTTGCTGTGCGCCACGTACATa                                                    <  1:732719/38‑1 (MQ=25)
                                      gtTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTa                    <  1:362857/69‑1 (MQ=255)
                                       ttATTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACg                              <  1:3142381/58‑1 (MQ=255)
                                         aatGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:2681914/3‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:104763/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:851071/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:1412219/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:1504120/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:151586/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:1649469/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:1235235/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:2327237/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:380702/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:1188129/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:907920/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:2672094/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:3145887/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:94807/1‑63 (MQ=255)
                                         aTTGATATCTTCCATGCTATGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg                       >  1:3071362/1‑63 (MQ=255)
                                                    ccATGCTGTGCGCCACGGACATAAAGCCCGCTTCa                                        >  1:2599906/1‑35 (MQ=255)
                                                        gCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTaaa                  >  1:2709888/1‑53 (MQ=255)
                                                        gCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:2390067/1‑69 (MQ=255)
                                                        gCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:3381199/1‑69 (MQ=255)
                                                        gCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:374415/1‑69 (MQ=255)
                                                        gCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:3308876/1‑69 (MQ=255)
                                                        gCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:2869129/1‑69 (MQ=255)
                                                        gCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:1302646/1‑69 (MQ=255)
                                                        gCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:703016/1‑69 (MQ=255)
                                                         cTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:381353/1‑68 (MQ=255)
                                                         cTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:828766/1‑68 (MQ=255)
                                                         cTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:1741166/1‑68 (MQ=255)
                                                         cTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:1917678/1‑68 (MQ=255)
                                                         cTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:2585512/1‑68 (MQ=255)
                                                         cTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:330597/1‑68 (MQ=255)
                                                         cTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:3073681/1‑68 (MQ=255)
                                                         cTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGcat  >  1:2821301/1‑68 (MQ=255)
                                                           |                                                                 
CACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCTTCGTCCAGCAT  >  minE/152428‑152552

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: