Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 64,573 G→A 15.1% I344I (ATC→ATT leuC ← 3‑isopropylmalate isomerase subunit, dehydratase component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE64,5730GA15.1% 125.0 / 9.6 72I344I (ATC→ATTleuC3‑isopropylmalate isomerase subunit, dehydratase component
Reads supporting (aligned to +/- strand):  ref base G (19/42);  new base A (9/2);  total (28/44)
Fisher's exact test for biased strand distribution p-value = 2.36e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TTTCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAG  >  minE/64508‑64638
                                                                 |                                                                 
tttCGCCCTTTTGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                <  1:2787736/69‑1 (MQ=255)
tttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                <  1:658168/69‑1 (MQ=255)
tttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                <  1:764481/69‑1 (MQ=255)
tttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                <  1:2448277/69‑1 (MQ=255)
tttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                <  1:604970/69‑1 (MQ=255)
tttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                <  1:260158/69‑1 (MQ=255)
 ttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGATCCGATa                                                                >  1:2381568/1‑68 (MQ=255)
 ttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                >  1:3426296/1‑68 (MQ=255)
 ttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                >  1:2032505/1‑68 (MQ=255)
 ttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                <  1:134610/68‑1 (MQ=255)
 ttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                >  1:2793775/1‑68 (MQ=255)
 ttCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATa                                                                >  1:1132176/1‑68 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:2664104/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:1023497/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:2552601/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:2513944/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:2809646/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:2876660/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:2241994/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:2996994/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:2089880/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:420232/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:1105205/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:1248787/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:1817512/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:1364279/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:1025600/68‑1 (MQ=255)
            gCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:1659348/68‑1 (MQ=255)
             cGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTc                                                     <  1:538892/67‑1 (MQ=255)
                 ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATaa                                                               >  1:964297/1‑53 (MQ=255)
                 ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAg                                                >  1:1097119/1‑68 (MQ=255)
                 ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGc                                               >  1:2756496/1‑69 (MQ=255)
                 ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGc                                               >  1:261998/1‑69 (MQ=255)
                 ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGc                                               >  1:1254817/1‑69 (MQ=255)
                 ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGc                                               >  1:3127967/1‑69 (MQ=255)
                 ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGc                                               >  1:211331/1‑69 (MQ=255)
                 ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGc                                               >  1:423244/1‑69 (MQ=255)
                 ctcCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGc                                               >  1:281341/1‑69 (MQ=255)
                 ctcCGCTGCCGCGCGTAAATCTTCAATGCGCCAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGc                                               >  1:2845919/1‑69 (MQ=255)
                                  aaTCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGaa                              <  1:2708499/69‑1 (MQ=255)
                                  aaTCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGaa                              <  1:823226/69‑1 (MQ=255)
                                  aaTCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGaa                              <  1:473721/69‑1 (MQ=255)
                                  aaTCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGaa                              <  1:617921/69‑1 (MQ=255)
                                            gcgcGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGat                              <  1:2038585/59‑2 (MQ=255)
                                             cgcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:1053688/56‑1 (MQ=255)
                                             cgcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:3089378/56‑1 (MQ=255)
                                             cgcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:352623/56‑1 (MQ=255)
                                             cgcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:207779/56‑1 (MQ=255)
                                             cgcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:2823813/56‑1 (MQ=255)
                                             cgcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:271187/56‑1 (MQ=255)
                                             cgcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:1350160/56‑1 (MQ=255)
                                             cgcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:2501226/56‑1 (MQ=255)
                                             cgcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:2484070/56‑1 (MQ=255)
                                             cgcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:677895/56‑1 (MQ=255)
                                             cgcgAGTTAGTACAGGAACCAATAAACACTTTGTc                                                     <  1:2105367/35‑1 (MQ=255)
                                             cgcgAGTTAGTACAGGAACCAATAAACACTTTGTc                                                     <  1:1921292/35‑1 (MQ=255)
                                              gcgAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCgg                                <  1:838394/55‑1 (MQ=255)
                                                    tGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCACCGGAATACCCGGTTTCAGCCCCa            <  1:1160215/69‑1 (MQ=255)
                                                       tACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGAATACCCGGTTTCAGCCCCata          <  1:256907/68‑1 (MQ=255)
                                                         cAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAg       >  1:731690/1‑69 (MQ=255)
                                                           ggAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAg       >  1:2101661/1‑67 (MQ=255)
                                                           ggAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAg       >  1:1768964/1‑67 (MQ=255)
                                                           ggAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAg       >  1:1141145/1‑67 (MQ=255)
                                                              aCCAATAAACACTTTGTCGATAGCCACTTCTGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCa   >  1:379482/1‑68 (MQ=255)
                                                              aCCAATAAACACTTTGTCGATAGCCACTTCTGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCa   >  1:1905648/1‑68 (MQ=255)
                                                              aCCAATAAACACTTTGTCGATAGCCACTTCTGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCa   >  1:1777248/1‑68 (MQ=255)
                                                              aCCAATAAACACTTTGTCGATAGCCACTTCTGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAg  >  1:289588/1‑69 (MQ=255)
                                                              aCCAATAAACACTTTGTCGATAGCCACTTCTGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAg  >  1:2493611/1‑69 (MQ=255)
                                                              aCCAATAAACACTTTGTCGATAGCCACTTCTGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAg  >  1:1687546/1‑69 (MQ=255)
                                                              aCCAATAAACACTTTGTCGATAGCCACTTCTGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAg  >  1:1470651/1‑69 (MQ=255)
                                                              aCCAATAAACACTTTGTCGATAGCCACTTCTGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAg  >  1:985445/1‑69 (MQ=255)
                                                              aCCAATAAACACTTTGTCGATAGCCACTTCTGTCAGCGGAATACCCGGTTTCAGCACCa            >  1:375345/1‑59 (MQ=255)
                                                                 |                                                                 
TTTCGCCCTTTGGCGATCTCCGCTGCCGCGCGTAAATCTTCAATGCGCGAGTTGGTACAGGAACCGATAAACACTTTGTCGATAGCCACTTCGGTCAGCGGAATACCCGGTTTCAGCCCCATATAGGCCAG  >  minE/64508‑64638

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: