Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,039,447 T→C 39.8% A156A (GCA→GCG ribB ← 3,4 dihydroxy‑2‑butanone‑4‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,039,4470TC39.8% 48.5 / 80.2 113A156A (GCA→GCGribB3,4 dihydroxy‑2‑butanone‑4‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (29/39);  new base C (35/10);  total (64/49)
Fisher's exact test for biased strand distribution p-value = 2.37e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 113 total aligned reads displayed.

GTCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCC  >  minE/2039381‑2039512
                                                                  |                                                                 
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:1462779/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:6027/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:517905/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:43138/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:423965/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:419976/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:3062168/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:1698005/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:26208/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:2362080/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:1264643/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:1059016/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGACAGCGTCATCAGATCAATAGTTGc                                                                 >  1:2527239/1‑69 (MQ=255)
gtCATTAGTCAGCTCACACAGCACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGc                                                                 >  1:136824/1‑69 (MQ=255)
     tAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAg                                                            >  1:2397273/1‑69 (MQ=255)
     tAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAg                                                            <  1:2483622/69‑1 (MQ=255)
     tAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAg                                                            <  1:2006457/69‑1 (MQ=255)
     tAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAg                                                            <  1:2629039/69‑1 (MQ=255)
     tAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAg                                                            <  1:1526684/69‑1 (MQ=255)
     tAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAg                                                            >  1:942460/1‑69 (MQ=255)
     tAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAg                                                            <  1:1088894/69‑1 (MQ=255)
      aGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAg                                                            <  1:2966667/68‑1 (MQ=255)
      aGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGt                                                           >  1:3075702/1‑69 (MQ=255)
      aGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGt                                                           >  1:720559/1‑69 (MQ=255)
         cAGCTCACACAGTACACCAGCCTGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTAtt                                                        <  1:1746555/69‑2 (MQ=255)
         cAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATg                                                        <  1:1262232/69‑1 (MQ=255)
         cAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATg                                                        <  1:1407863/69‑1 (MQ=255)
               acacaGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGc                                                   <  1:2849459/68‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:2073492/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:2265927/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:2263251/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:2262535/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:1433078/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:342002/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:37254/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:2155154/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:1244698/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:1178563/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:2914886/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:196988/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:795059/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:1083583/69‑1 (MQ=255)
                      acacCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTc                                           <  1:1735919/69‑1 (MQ=255)
                         ccAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATTAATAGTTGCTTCAGTATGACCGCCACGCGTCAg                                         <  1:1423661/68‑1 (MQ=255)
                         ccAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTCAg                                         <  1:2357450/68‑1 (MQ=255)
                         ccAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTCAg                                         <  1:2556275/68‑1 (MQ=255)
                         ccAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTCAg                                         <  1:3089878/68‑1 (MQ=255)
                         ccAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTCAg                                         <  1:1422315/68‑1 (MQ=255)
                         ccAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTCAg                                         <  1:805594/68‑1 (MQ=255)
                         ccAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTCAg                                         <  1:815406/68‑1 (MQ=255)
                                 tttAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAg                                                            <  1:262626/41‑1 (MQ=255)
                                                    atcaGATCAATAGTTGCTTCAGTATGACCGCCAcgc                                              >  1:40483/1‑36 (MQ=255)
                                                    atcaGATCAATAGTTGCTTCAGTATGACCGCCAcgc                                              >  1:427933/1‑36 (MQ=255)
                                                    atcaGATCAATAGTTGCTTCAGTATGACCGCCAcgc                                              >  1:1504813/1‑36 (MQ=255)
                                                    atcaGATCAATAGTTGCTTCAGTATGACCGCCAcgc                                              >  1:2499151/1‑36 (MQ=255)
                                                    atcaGATCAATAGTTGCTTCAGTATGACCGCCAcgc                                              >  1:1540101/1‑36 (MQ=255)
                                                    atcaGATCAATAGTTGCTTCAGTATGACCGCCAcgc                                              >  1:2033998/1‑36 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGCCCGCCACGCGTCAGTACACCACCTGCCTgcgca                      >  1:644066/1‑55 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGa                   >  1:3205240/1‑59 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGt                >  1:512270/1‑62 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:2114579/1‑69 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:697736/1‑69 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:546245/1‑69 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:522846/1‑69 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:1259955/1‑69 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:3217824/1‑69 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:3123109/1‑69 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:2482290/1‑69 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:3051977/1‑69 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:1275433/1‑69 (MQ=255)
                                                        gATCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAACTGGGaaaa         >  1:3054708/1‑69 (MQ=255)
                                                        gATCGAAAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         >  1:435884/1‑69 (MQ=255)
                                                         aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         <  1:44283/68‑1 (MQ=255)
                                                         aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         <  1:466737/68‑1 (MQ=255)
                                                         aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         <  1:1112096/68‑1 (MQ=255)
                                                         aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         <  1:2037603/68‑1 (MQ=255)
                                                         aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         <  1:3015684/68‑1 (MQ=255)
                                                         aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         <  1:2176274/68‑1 (MQ=255)
                                                         aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         <  1:501995/68‑1 (MQ=255)
                                                         aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         <  1:2628769/68‑1 (MQ=255)
                                                         aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         <  1:3327770/68‑1 (MQ=255)
                                                         aTCGATAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGaaaa         <  1:1084050/68‑1 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCAcctgcc                            >  1:118063/3‑47 (MQ=39)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACg       >  1:284666/3‑68 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:1686953/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:1646717/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:657071/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:1469614/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:2277825/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:1823740/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:1073749/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:531306/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:1326178/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:1882218/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:2336842/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:2214028/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:3361660/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:3230613/3‑69 (MQ=255)
                                                           cgaTAGTCGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGCGCACGAAGTGGGAAAACGt      >  1:224411/3‑69 (MQ=255)
                                                                gTTGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGcc  >  1:2836692/1‑68 (MQ=255)
                                                                  |                                                                 
GTCATTAGTCAGCTCACACAGTACACCAGCCGGTTTAAAGCCTGCCAGCGTCATCAGATCAATAGTTGCTTCAGTATGACCGCCACGCGTCAGTACACCACCTGCCTGAGCGCGAAGTGGGAAAACGTGGCC  >  minE/2039381‑2039512

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: