Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,571,225	A→T	16.4%	D377V (GAT→GTT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,571,225	0	A	T	16.4%	132.2 / 19.2	55	D377V (GAT→GTT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/12);  new base T (3/6);  total (37/18)
Fisher's exact test for biased strand distribution p-value = 4.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTAATTATTTTTGCATAT  >  minE/2571155‑2571291                                                                       |                                                                   ttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGa                                                                    >  1:2111540/1‑71 (MQ=255)  tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                                                   >  1:2770059/1‑71 (MQ=255)    aGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                                                                 <  1:999983/71‑1 (MQ=255)    aGCCATGTTT‑GGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                                                                <  1:326249/71‑1 (MQ=255)      ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                               <  1:3915486/71‑1 (MQ=255)      ccATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCa                                                                 <  1:1846212/69‑1 (MQ=255)       cATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtt                                                               >  1:77267/1‑70 (MQ=255)       cATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtt                                                               >  1:3698512/1‑70 (MQ=255)        aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGTTCATTga                                                             <  1:1162533/71‑1 (MQ=255)        aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                               >  1:1513791/1‑69 (MQ=255)        aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                                               >  1:3023694/1‑69 (MQ=255)         tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGTTCATTgag                                                            <  1:1721879/71‑1 (MQ=255)         tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgag                                                            <  1:2476756/71‑1 (MQ=255)          gtctCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAt                                                                >  1:2826972/4‑66 (MQ=255)          gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagg                                                           >  1:1719969/1‑71 (MQ=255)            ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                                         >  1:3113909/1‑71 (MQ=255)             tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                                                                >  1:2674656/1‑63 (MQ=255)             tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                                         >  1:1621007/1‑70 (MQ=255)             tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                                         >  1:1591226/1‑70 (MQ=255)             tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                                         <  1:2331410/70‑1 (MQ=255)               ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTTGATCATTgagggggag                                                      <  1:1659192/71‑1 (MQ=255)               ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggggag                                                      <  1:2609564/71‑1 (MQ=255)               ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                                                       <  1:3515525/70‑1 (MQ=255)                    tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGagggggaggggg                                                  >  1:3676318/1‑70 (MQ=255)                    tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGagggggaggggg                                                  >  1:2541047/1‑70 (MQ=255)                        gAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGTTCATTGAGGGGGAGGGGGGAGAt                                             <  1:965757/71‑1 (MQ=255)                        gAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagg                                                           <  1:3503709/57‑1 (MQ=255)                         aaaGTGCTGCTAAGAATGACCAGCTTGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt                                             >  1:3500077/1‑70 (MQ=255)                          aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                                         >  1:737039/1‑57 (MQ=255)                               ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggga                                                       >  1:71673/1‑54 (MQ=255)                               ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgag                                                            >  1:2333096/1‑49 (MQ=255)                                  ctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGTTCATTGAGGGGGAGGGGGGAGATTATGGAAc                                     <  1:1342200/69‑1 (MQ=255)                                  ctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATgg                                        >  1:2710222/1‑66 (MQ=255)                                  ctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACtt                                   <  1:3279660/71‑1 (MQ=255)                                    aaGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggggaggg                                                    >  1:1532508/1‑52 (MQ=255)                                    aaGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTg                                 >  1:1522160/1‑71 (MQ=255)                                      gAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCt                               >  1:2036300/1‑71 (MQ=255)                                        aTGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTcc                             >  1:2325304/1‑71 (MQ=255)                                         tGACCAGCTGGATGCATTATCTTTATTTGGTTCATTga                                                             <  1:994093/38‑1 (MQ=255)                                         tGACCAGCTGGATGCATTATCTTTATTTGGTTCATTga                                                             <  1:98407/38‑1 (MQ=255)                                         tGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAgggggaggggggag                                               >  1:232550/1‑52 (MQ=255)                                         tGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCt                            >  1:3083405/1‑71 (MQ=255)                                           aCCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTg                          >  1:1715649/1‑71 (MQ=255)                                                  ggATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTg                                 >  1:2817668/1‑57 (MQ=255)                                                   gATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGAc                   >  1:1655067/1‑70 (MQ=255)                                                    aTGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGAc                   >  1:1767080/1‑69 (MQ=255)                                                     tGCATTATCTTTATTTGGTTCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGc                                >  1:1528940/1‑55 (MQ=255)                                                     tGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGgaga                                              >  1:340802/1‑41 (MQ=255)                                                          tATCTTTATTTGGTTCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTAattatt           >  1:402068/1‑71 (MQ=255)                                                           aTCTTTATTTGGTTCATTGAGGGGGAGGGGGGAGATTa                                           >  1:1101065/1‑38 (MQ=255)                                                           aTCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGaa                                      <  1:502394/43‑1 (MQ=255)                                                           aTCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTAATTAttt          >  1:2972643/1‑71 (MQ=255)                                                            tCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTAATTAtttt         >  1:111600/1‑71 (MQ=255)                                                             cTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTAATTAttttt        >  1:3908975/1‑71 (MQ=255)                                                             cTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTAATTAttttt        >  1:2725510/1‑71 (MQ=255)                                                                    tGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTAATTATTTTTGCatat  >  1:2946507/1‑70 (MQ=255)                                                                       |                                                                   TTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAGGACTAATTATTTTTGCATAT  >  minE/2571155‑2571291

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