Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,052	A→T	17.5%	I11I (ATT→ATA)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,052	0	A	T	17.5%	190.0 / 32.3	80	I11I (ATT→ATA)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base T (6/8);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGTT  >  minE/454983‑455122                                                                      |                                                                       gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       <  1:4334378/71‑1 (MQ=255) gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       >  1:328063/1‑71 (MQ=255) gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       <  1:3352644/71‑1 (MQ=255) gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       <  1:3796277/71‑1 (MQ=255) gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       <  1:4845549/71‑1 (MQ=255) gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       <  1:762769/71‑1 (MQ=255) gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       <  1:2087903/71‑1 (MQ=255) gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       <  1:2283969/71‑1 (MQ=255)  agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatat                                                                      <  1:386759/71‑1 (MQ=255)  agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       >  1:4696071/1‑70 (MQ=255)     aTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                                                                    <  1:1708427/70‑1 (MQ=255)       cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc                                                                 >  1:4318729/1‑71 (MQ=255)       cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc                                                                 >  1:2903474/1‑71 (MQ=255)        agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                                                >  1:626790/1‑71 (MQ=255)        agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                                                >  1:1768082/1‑71 (MQ=255)         gcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       >  1:494033/1‑63 (MQ=255)          taaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATcgcgcgcg                                                              >  1:4027102/4‑71 (MQ=255)          cagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                                                              >  1:925877/1‑71 (MQ=255)           agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                                             >  1:1929373/1‑71 (MQ=255)           agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                                             >  1:2968215/1‑71 (MQ=255)           aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                                                             >  1:1049272/3‑71 (MQ=255)           aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                                                             >  1:7948/3‑71 (MQ=255)              gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                                                          <  1:3150876/71‑1 (MQ=255)               cACGCATGCCGTAAAAGCTAAAATATTTCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                         <  1:1217570/71‑1 (MQ=255)               cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                                                                    >  1:957106/1‑60 (MQ=255)               cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                         <  1:4503315/71‑1 (MQ=255)               cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                         <  1:2786885/71‑1 (MQ=255)               cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                         <  1:4509645/71‑1 (MQ=255)               cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                                                          <  1:1933962/70‑1 (MQ=255)               cACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                                                         <  1:3778408/71‑1 (MQ=255)                 cGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                                       <  1:70855/71‑1 (MQ=255)                 cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa                                                       <  1:294826/71‑1 (MQ=255)                  gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                                      >  1:2119683/1‑70 (MQ=255)                   cATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGTTGTGATg                                                     >  1:3128724/1‑71 (MQ=255)                   cATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                                     <  1:121246/71‑3 (MQ=255)                     tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatat                                                                      <  1:4975880/52‑1 (MQ=255)                     tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                                     <  1:1998365/69‑3 (MQ=255)                       ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtg                                                 >  1:442478/1‑71 (MQ=255)                          tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                       <  1:4982778/46‑1 (MQ=255)                          tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                                                         <  1:5459242/60‑1 (MQ=255)                            aaaaCTGAAGTACTCCCAGATTTGTATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                                      >  1:3834926/1‑60 (MQ=255)                            aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc                                            <  1:1798406/70‑3 (MQ=255)                            aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc                                            <  1:4056179/70‑3 (MQ=255)                            aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc                                            <  1:3652356/70‑3 (MQ=255)                            aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                                                                    >  1:1839905/1‑47 (MQ=255)                            aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAt                                            >  1:3473302/1‑71 (MQ=255)                               gCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGc                                                            <  1:2532398/52‑1 (MQ=255)                                 tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                         <  1:515501/53‑1 (MQ=255)                                  aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                                             >  1:4483838/1‑64 (MQ=255)                                   aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa                                             >  1:5443618/1‑63 (MQ=255)                                   aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCataat                                     <  1:4807160/71‑1 (MQ=255)                                      tatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGtt                                               <  1:1532654/58‑1 (MQ=255)                                        ttCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc                                <  1:4033631/71‑1 (MQ=255)                                           ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTcc                               >  1:4394748/1‑69 (MQ=255)                                           ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcac                             <  1:5599699/71‑1 (MQ=255)                                           ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca                              >  1:3020669/1‑70 (MQ=255)                                           ccAGATTTGCAGAGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca                              >  1:1175136/1‑70 (MQ=255)                                            cAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcaca                            <  1:1510223/71‑1 (MQ=255)                                             aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                                       >  1:2925206/1‑43 (MQ=255)                                             aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                                       >  1:3796198/1‑43 (MQ=255)                                             aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACaa                           >  1:4977867/1‑71 (MQ=255)                                              gATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAt                          >  1:5553338/1‑71 (MQ=255)                                               aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         <  1:5153035/71‑1 (MQ=255)                                               aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         >  1:1924153/1‑71 (MQ=255)                                               aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         <  1:3458865/71‑1 (MQ=255)                                                 ttGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                       >  1:5031569/1‑71 (MQ=255)                                                 ttGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                       >  1:2460373/1‑71 (MQ=255)                                                 ttGCAGAGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                       >  1:737611/1‑71 (MQ=255)                                                   gCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                       <  1:5020296/69‑1 (MQ=255)                                                    cAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACgg                    >  1:3548632/1‑71 (MQ=255)                                                     aGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc                   <  1:1001850/71‑1 (MQ=255)                                                      gcgcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc                   <  1:4602369/70‑1 (MQ=255)                                                       cgcCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc                                            >  1:2435057/1‑41 (MQ=12)                                                       cgcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCaa                 >  1:4018329/1‑71 (MQ=255)                                                           aCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGtt                                               <  1:940477/37‑1 (MQ=255)                                                           aCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTcc                               <  1:1744881/53‑1 (MQ=255)                                                           aCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTcc                               <  1:2095013/53‑1 (MQ=255)                                                            cGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         <  1:3388961/58‑1 (MQ=255)                                                            cGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         <  1:5311807/58‑1 (MQ=255)                                                             gACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcac                             >  1:1080846/1‑53 (MQ=255)                                                             gACATAATAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCAc           <  1:1200868/71‑1 (MQ=255)                                                              aCATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTc             <  1:4935546/68‑1 (MQ=255)                                                               cATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGt         <  1:4355098/71‑1 (MQ=255)                                                                aTAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTa        <  1:5211152/71‑1 (MQ=255)                                                                aTAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTa        <  1:1986489/71‑1 (MQ=255)                                                                   gTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCataat                                     <  1:923398/39‑1 (MQ=255)                                                                   gTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       >  1:5243582/1‑69 (MQ=255)                                                                   gTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       >  1:3092661/1‑69 (MQ=255)                                                                    tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCAtaataa                                   >  1:389757/1‑40 (MQ=255)                                                                    tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCAtaataa                                   >  1:512561/1‑40 (MQ=255)                                                                    tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       <  1:4625531/68‑1 (MQ=255)                                                                    tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       <  1:3794278/68‑1 (MQ=255)                                                                    tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       <  1:1133868/68‑1 (MQ=255)                                                                     aaaTAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGt   >  1:2178223/1‑71 (MQ=255)                                                                      aaTAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGtt  >  1:457941/1‑71 (MQ=255)                                                                      |                                                                       GAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGTT  >  minE/454983‑455122

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