Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	571,482	C→T	35.1%	A25A (GCG→GCA)	ybjE ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	571,482	0	C	T	35.1%	66.0 / 50.8	91	A25A (GCG→GCA)	ybjE	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (35/24);  new base T (17/15);  total (52/39)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

ACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACA  >  minE/571422‑571543                                                             |                                                              aCCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCtt                                                        <  1:2271470/68‑1 (MQ=255) aCCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCtt                                                        <  1:262683/68‑1 (MQ=255)    cATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCttgttg                                                    >  1:192473/1‑69 (MQ=255)    cATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCttgttg                                                    >  1:1197802/1‑69 (MQ=255)    cATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCttgttg                                                    >  1:765305/1‑69 (MQ=255)        aaaaaGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACg                                                              >  1:3213368/1‑55 (MQ=255)        aaaaaGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGa                                                >  1:2413375/1‑69 (MQ=255)        aaaaaGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGa                                                >  1:251363/1‑69 (MQ=255)        aaaaaGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGa                                                >  1:1948337/1‑69 (MQ=255)        aaaaaGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGa                                                >  1:1959839/1‑69 (MQ=255)        aaaaaGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGa                                                >  1:1201667/1‑69 (MQ=255)        aaaaaGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGa                                                >  1:780154/1‑69 (MQ=255)        aaaaaGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGa                                                >  1:833322/1‑69 (MQ=255)        aaaaaGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGa                                                >  1:2803601/1‑69 (MQ=255)                   aaGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGa                                      >  1:3001226/1‑68 (MQ=255)                   aaGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGa                                      >  1:3049123/1‑68 (MQ=255)                   aaGGTAAACCATACAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGa                                      >  1:3133401/1‑68 (MQ=255)                    aGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAg                                     >  1:2241009/1‑68 (MQ=255)                    aGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAg                                     >  1:2163737/1‑68 (MQ=255)                    aGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAg                                     >  1:3333647/1‑68 (MQ=255)                    aGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAg                                     >  1:1907182/1‑68 (MQ=255)                    aGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAg                                     >  1:336545/1‑68 (MQ=255)                    aGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAg                                     >  1:2427198/1‑68 (MQ=255)                    aGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAg                                     >  1:1875555/1‑68 (MQ=255)                    aGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAg                                     >  1:1705916/1‑68 (MQ=255)                      gTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTa                                  >  1:1151650/1‑69 (MQ=255)                      gTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTa                                  >  1:2666156/1‑69 (MQ=255)                      gTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTa                                  >  1:150395/1‑69 (MQ=255)                       tAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTaa                                 >  1:482639/1‑69 (MQ=255)                       tAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTaa                                 >  1:2566483/1‑69 (MQ=255)                       tAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTaa                                 >  1:1495842/1‑69 (MQ=255)                          aCCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCgg                                           <  1:1675426/56‑1 (MQ=255)                           ccATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAAccc                              >  1:304462/1‑68 (MQ=255)                           ccATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAAccc                              >  1:2568203/1‑68 (MQ=255)                           ccATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAAccc                              >  1:16289/1‑68 (MQ=255)                           ccATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAAccc                              >  1:1246572/1‑68 (MQ=255)                           ccATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAAccc                              >  1:278054/1‑68 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:2470185/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:2733309/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:2654738/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:2466331/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:2860802/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:2808787/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:2179065/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:1045206/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:1299334/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:139038/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:1543131/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:554885/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:470865/50‑1 (MQ=255)                                    ttAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:330081/50‑1 (MQ=255)                                     tAATAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGAAGCGGAATg                                       <  1:2264307/49‑1 (MQ=38)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:584592/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:434888/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:41419/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:654384/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:907635/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:981645/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:10238/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:991767/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:2674485/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:1037530/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:1075767/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:118974/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:1637579/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:1639822/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:1945876/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:2658674/1‑39 (MQ=37)                                      aaTAGCTGATTAATAACTCTTAATGCAGCTTTTTGGCGa                                                >  1:261663/1‑39 (MQ=37)                                            tGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAg             <  1:827633/68‑1 (MQ=255)                                            tGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAg             <  1:1365397/68‑1 (MQ=255)                                            tGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAg             <  1:694936/68‑1 (MQ=255)                                            tGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAg             <  1:2923059/68‑1 (MQ=255)                                            tGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAg             <  1:501417/68‑1 (MQ=255)                                            tGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAg             <  1:1177407/68‑1 (MQ=255)                                            tGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAg             <  1:1890759/68‑1 (MQ=255)                                            tGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAg             <  1:1919288/68‑1 (MQ=255)                                                 aataaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCGACAATCAGGGGAACCAGAATGa        >  1:2666003/1‑68 (MQ=255)                                                   taaCTTTTAACGCAGCTTGTTGGCGAAGCGTAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:1384364/68‑1 (MQ=255)                                                   taaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:1863633/68‑1 (MQ=255)                                                   taaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:1683770/68‑1 (MQ=255)                                                   taaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:3244374/68‑1 (MQ=255)                                                   taaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:1066977/68‑1 (MQ=255)                                                   taaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:2522146/68‑1 (MQ=255)                                                    aaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:205280/67‑1 (MQ=255)                                                    aaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:747250/67‑1 (MQ=255)                                                    aaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:1690334/67‑1 (MQ=255)                                                    aaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:2528458/67‑1 (MQ=255)                                                    aaCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:2030062/67‑1 (MQ=255)                                                     aCTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGAtt      <  1:1865230/66‑1 (MQ=255)                                                       ttttAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACa  <  1:2726675/68‑1 (MQ=255)                                                             |                                                              ACCCATAAAAAAGAGAATAAGGTAAACCATCCAGCTTAATAGCTGATTAATAACTTTTAACGCAGCTTGTTGGCGAAGCGGAATGAGGTAACCCACAATCAGGGGAACCAGAATGATTAACA  >  minE/571422‑571543

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: ‑