Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,055	A→C	16.1%	A10A (GCT→GCG)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,055	0	A	C	16.1%	158.1 / 31.9	81	A10A (GCT→GCG)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (34/34);  new base C (5/8);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

ATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGTT  >  minE/454987‑455122                                                                     |                                                                    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)      taaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATcgcgcgcg                                                              >  1:4027102/4‑71 (MQ=255)      cagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                                                              >  1:925877/1‑71 (MQ=255)       agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                                             >  1:2968215/1‑71 (MQ=255)       agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                                                             >  1:1929373/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)                 tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg                                                     <  1:1998365/69‑3 (MQ=255)                   ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtg                                                 >  1:442478/1‑71 (MQ=255)                      tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt                                                         <  1:5459242/60‑1 (MQ=255)                        aaaaCTGAAGTACTCCCAGATTTGTATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                                      >  1:3834926/1‑60 (MQ=255)                        aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc                                            <  1:4056179/70‑3 (MQ=255)                        aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc                                            <  1:1798406/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:2095013/53‑1 (MQ=255)                                                       aCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTcc                               <  1:1744881/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:512561/1‑40 (MQ=255)                                                                tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCAtaataa                                   >  1:389757/1‑40 (MQ=255)                                                                tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       <  1:3794278/68‑1 (MQ=255)                                                                tAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTaa       <  1:4625531/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)                                                                     aGCGCGCTGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                         <  1:3459574/45‑1 (MQ=255)                                                                     aGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATa     >  1:5512767/1‑65 (MQ=255)                                                                     aGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATa     >  1:3374281/1‑65 (MQ=255)                                                                     aGCGCGCGGCTGTGATGCGTGTTAATTCATAATAATTCCACAAtt                         >  1:3020230/1‑45 (MQ=255)                                                                     |                                                                    ATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGTT  >  minE/454987‑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: ‑