Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	484,420	C→T	12.9%	L76L (CTG→TTG)	tolB →	periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	484,420	0	C	T	12.9%	247.1 / 22.4	93	L76L (CTG→TTG)	tolB	periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (48/33);  new base T (9/3);  total (57/36)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAT  >  minE/484353‑484481                                                                    |                                                              tatTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGc                                                            >  1:5032444/1‑71 (MQ=255)   ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:3968939/71‑1 (MQ=255)   ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:5024136/71‑1 (MQ=255)   ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:3934714/71‑1 (MQ=255)   ttGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:1911562/71‑1 (MQ=255)    tGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCca                                                          <  1:415750/70‑1 (MQ=255)         gcATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCAcagcag                                                    <  1:2230236/71‑1 (MQ=255)          cATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACagcagc                                                   <  1:3470395/71‑1 (MQ=255)           aTCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACagcagc                                                   <  1:4497911/70‑1 (MQ=255)           aTCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGcc                                                  <  1:304255/71‑1 (MQ=255)           aTCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGcc                                                  <  1:3116447/71‑1 (MQ=255)               ttGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGt                                              <  1:134081/71‑1 (MQ=255)                     ctgACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTgc                                         >  1:3200793/1‑70 (MQ=255)                     ctgACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTgc                                         <  1:1512489/70‑1 (MQ=255)                           tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:3322310/1‑69 (MQ=255)                           tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:5124175/1‑69 (MQ=255)                           tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:2568564/1‑69 (MQ=255)                           tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:417388/1‑69 (MQ=255)                           tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:1405030/1‑69 (MQ=255)                           tGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    >  1:1248139/1‑69 (MQ=255)                            gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGa                                   <  1:4171059/69‑1 (MQ=255)                            gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGa                                   <  1:3794876/69‑1 (MQ=255)                            gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGa                                   <  1:3697408/69‑1 (MQ=255)                            gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAg                                 >  1:1386852/1‑71 (MQ=255)                            gCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAg                                 >  1:3111809/1‑71 (MQ=255)                             cGTAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGt                                              <  1:4077991/57‑1 (MQ=255)                             cGTAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGt                                              <  1:4174140/57‑1 (MQ=255)                             cGTAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGt                                              <  1:4765968/57‑1 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCAcagca                                                     >  1:2256892/1‑48 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTa                                             >  1:986446/1‑56 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:188709/1‑57 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:2316276/1‑57 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:3124706/1‑57 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:4326772/1‑57 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:2189797/1‑57 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:1793924/1‑57 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGTGCTCGTTTGCCACAGCAGCCGGGTAg                                            >  1:3605645/1‑57 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTAc                              <  1:3891415/71‑1 (MQ=255)                               tAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTAc                              <  1:3475548/71‑1 (MQ=255)                                aaCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTg                                                             >  1:2156493/1‑39 (MQ=255)                                aaCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTAc                              >  1:3781454/1‑70 (MQ=255)                                 aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:248188/1‑43 (MQ=255)                                 aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:5114070/1‑43 (MQ=255)                                 aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:1080455/1‑43 (MQ=255)                                 aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:3999399/1‑43 (MQ=255)                                 aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:1439391/1‑43 (MQ=255)                                 aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:474664/1‑43 (MQ=255)                                 aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:3312105/1‑43 (MQ=255)                                 aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:2354194/1‑43 (MQ=255)                                 aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:767597/1‑43 (MQ=255)                                 aCAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCcaca                                                        >  1:3026030/1‑43 (MQ=255)                                   aGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:4912169/71‑1 (MQ=255)                                       gTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:3290992/67‑1 (MQ=255)                                       gTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:5161206/67‑1 (MQ=255)                                       gTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAcc                          <  1:4653797/67‑1 (MQ=255)                                          aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    <  1:4182551/54‑1 (MQ=255)                                          aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    <  1:531551/54‑1 (MQ=255)                                          aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAgg                                    <  1:5189153/54‑1 (MQ=255)                                          aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACaa                            <  1:172070/62‑1 (MQ=255)                                          aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACaa                            <  1:4624919/62‑1 (MQ=255)                                          aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCa                   >  1:572602/1‑71 (MQ=255)                                          aaTTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCa                   >  1:217982/1‑71 (MQ=255)                                            tttAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAAc                           >  1:3086278/1‑61 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGcc                                                  >  1:3363271/1‑35 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1408754/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:537950/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:627231/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:706034/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:744349/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:76171/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:105250/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1582753/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:2102796/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:3733995/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:398656/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1571011/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1794182/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:378018/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:3799750/1‑71 (MQ=255)                                               aaTCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTc              >  1:1921216/1‑71 (MQ=255)                                                 tCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCg            >  1:536670/1‑71 (MQ=255)                                                   cGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGc           <  1:1703316/70‑1 (MQ=255)                                                     ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAg                                 <  1:1191845/46‑1 (MQ=255)                                                     ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACt        <  1:58573/71‑1 (MQ=255)                                                     ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACt        <  1:3689644/71‑1 (MQ=255)                                                     ttAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACt        <  1:3202075/71‑1 (MQ=255)                                                       aGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTgg      >  1:3944476/1‑71 (MQ=255)                                                       aGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTgg      >  1:3463048/1‑71 (MQ=255)                                                           cgcgCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAt  >  1:2375499/1‑71 (MQ=255)                                                           cgcgCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAt  >  1:802516/1‑71 (MQ=255)                                                           cgcgCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGAATGGTCCGCACTGGGCAt  >  1:1718268/1‑71 (MQ=255)                                                                   tCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCa                         <  1:3198920/40‑1 (MQ=255)                                                                   tCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCa                         <  1:2400304/40‑1 (MQ=255)                                                                    |                                                              TATTGGCGGCATCGTTGCTGCTGACTTGCGTAACAGCGGTAAATTTAATCCGTTAGATCGCGCTCGTCTGCCACAGCAGCCGGGTAGTGCGCAGGAAGTACAACCAGCTGCATGGTCCGCACTGGGCAT  >  minE/484353‑484481

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