Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	43,706	T→C	25.5%	G74G (GGT→GGC)	yaaU →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	43,706	0	T	C	25.5%	95.3 / 29.6	82	G74G (GGT→GGC)	yaaU	predicted transporter
Reads supporting (aligned to +/- strand):  ref base T (40/21);  new base C (9/12);  total (49/33)
Fisher's exact test for biased strand distribution p-value = 7.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

TGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGAT  >  minE/43640‑43771                                                                   |                                                                  tGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTta                                                                 <  1:2578000/69‑1 (MQ=255)                   gggCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCg                                                <  1:419611/68‑1 (MQ=255)                   gggCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCg                                                <  1:3340088/68‑1 (MQ=255)                   gggCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCg                                                <  1:248031/68‑1 (MQ=255)                   gggCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCg                                                <  1:837519/68‑1 (MQ=255)                   gggCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCg                                                <  1:1956391/68‑1 (MQ=255)                   gggCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCg                                                <  1:862948/68‑1 (MQ=255)                    ggCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGa                                              >  1:941102/1‑69 (MQ=255)                    ggCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGa                                              >  1:2649551/1‑69 (MQ=255)                           gAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGgcgc                                        >  1:521121/1‑68 (MQ=255)                           gAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGgcgc                                        >  1:2340525/1‑68 (MQ=255)                           gAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGgcgc                                        >  1:2152314/1‑68 (MQ=255)                           gAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGgcgc                                        >  1:2939592/1‑68 (MQ=255)                           gAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGgcgc                                        >  1:3037991/1‑68 (MQ=255)                           gAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGgcgc                                        >  1:1909849/1‑68 (MQ=255)                           gAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGgcgc                                        >  1:3142970/1‑68 (MQ=255)                           gAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGgcgc                                        >  1:3216126/1‑68 (MQ=255)                           gAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGgcg                                         >  1:2221137/1‑67 (MQ=255)                                cTCGCCGGGCTGTTCGTTGTCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATg                                  <  1:406901/69‑1 (MQ=255)                                cTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAAt                                   >  1:334212/1‑68 (MQ=255)                                cTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAAt                                   >  1:2373158/1‑68 (MQ=255)                                cTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATg                                  <  1:186004/69‑1 (MQ=255)                                cTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATg                                  <  1:3296941/69‑1 (MQ=255)                                cTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATg                                  <  1:1911336/69‑1 (MQ=255)                                cTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATg                                  >  1:699910/1‑69 (MQ=255)                                cTCGCCGGGCTGGTCGTTGGCACATCGCTGTTTGGTTATATTTCCGAt                                                       >  1:398391/1‑48 (MQ=255)                                     cGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTc                               <  1:3315357/67‑1 (MQ=255)                                     cGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTc                               <  1:1003306/67‑1 (MQ=255)                                     cGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTc                               <  1:424723/67‑1 (MQ=255)                                     cGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTc                               <  1:2721220/67‑1 (MQ=255)                                     cGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTc                               <  1:696335/67‑1 (MQ=255)                                     cGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTc                               <  1:2598187/67‑1 (MQ=255)                                      gggCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTcc                              >  1:1888890/1‑67 (MQ=255)                                      gggCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTc                            >  1:883206/1‑69 (MQ=255)                                      gggCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTc                            >  1:1747473/1‑69 (MQ=255)                                          tGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCAtt                         >  1:3166961/1‑68 (MQ=255)                                          tGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCAtt                         >  1:527871/1‑68 (MQ=255)                                          tGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCAtt                         >  1:881883/1‑68 (MQ=255)                                          tGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCAtt                         >  1:1682947/1‑68 (MQ=255)                                          tGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCAtt                         >  1:591546/1‑68 (MQ=255)                                          tGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCAtt                         >  1:645248/1‑68 (MQ=255)                                          tGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCAtt                         >  1:2767262/1‑68 (MQ=255)                                           gTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGa                       <  1:1256713/69‑1 (MQ=255)                                           gTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGa                       <  1:3320762/69‑1 (MQ=255)                                              cGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGa                                              >  1:537110/1‑43 (MQ=255)                                              cGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:183658/1‑60 (MQ=255)                                              cGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:18673/1‑60 (MQ=255)                                              cGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:1198172/1‑60 (MQ=255)                                              cGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:533164/1‑60 (MQ=255)                                              cGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:1551552/1‑60 (MQ=255)                                              cGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:3053791/1‑60 (MQ=255)                                              cGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:1370031/1‑60 (MQ=255)                                              cGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:856275/1‑60 (MQ=255)                                              cGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:1105903/1‑60 (MQ=255)                                              cGTTCGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:458697/1‑60 (MQ=255)                                                     acaTCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGcc             >  1:2624601/1‑69 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:2421739/59‑1 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:1167698/59‑1 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:1566720/59‑1 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:3427334/59‑1 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:3418541/59‑1 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:1695983/59‑1 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:2186188/59‑1 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:140598/59‑1 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:2962365/59‑1 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:3269807/59‑1 (MQ=255)                                                       aTCGCTGTTTGGCTATATCTACGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGata                     <  1:2342880/59‑1 (MQ=255)                                                          gCTGTTTGGTTATATTTCCGATAAAGTCGGACGgcgc                                        >  1:3204899/1‑37 (MQ=255)                                                          gCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:2892922/1‑48 (MQ=255)                                                          gCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:3434392/1‑48 (MQ=255)                                                          gCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:315131/1‑48 (MQ=255)                                                          gCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:3150911/1‑48 (MQ=255)                                                          gCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:2200501/1‑48 (MQ=255)                                                          gCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCgg        <  1:2741433/69‑1 (MQ=255)                                                          gCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCgg        <  1:1587041/69‑1 (MQ=255)                                                          gCTGTTAGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCt                             >  1:442401/1‑48 (MQ=255)                                                           gtGTTTGGCTATATCTCCGATAAAGTCGGACGGCGCATAATGTTCCTCATTGata                     <  1:1052513/54‑1 (MQ=255)                                                                tGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGat  >  1:324526/1‑69 (MQ=255)                                                                tGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGat  >  1:2135183/1‑69 (MQ=255)                                                                tGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGat  >  1:2082943/1‑69 (MQ=255)                                                                tGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGat  >  1:1405806/1‑69 (MQ=255)                                                                tGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGat  >  1:1047544/1‑69 (MQ=255)                                                                   |                                                                  TGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGTTATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGAT  >  minE/43640‑43771

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