Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,761	T→G	7.5%	I222M (ATT→ATG) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,761	0	T	G	7.5%	204.1 / 12.9	80	I222M (ATT→ATG) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base G (3/3);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

TCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGC  >  minE/2570696‑2570828                                                                  |                                                                    tCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                                >  1:3058713/1‑71 (MQ=255)  cGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAtt                                                                     >  1:1788102/1‑65 (MQ=255)   gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                              >  1:2690411/1‑71 (MQ=255)   gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                              >  1:23045/1‑71 (MQ=255)   gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                              >  1:3660616/1‑71 (MQ=255)      gATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGa                                                           <  1:2515169/71‑1 (MQ=255)        ttGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                         <  1:1230583/71‑1 (MQ=255)         tgtgAGCTTAAGAGCGCAACTTACTGAAGCATAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                         >  1:1780225/1‑70 (MQ=255)         tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATACCCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                         >  1:3345346/1‑70 (MQ=255)         tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                         <  1:3676153/70‑1 (MQ=255)         tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                         <  1:1837824/70‑1 (MQ=255)         tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                         <  1:4065684/70‑1 (MQ=255)           tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                       <  1:3286176/70‑1 (MQ=255)           tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAAt                                                      >  1:750207/1‑71 (MQ=255)             aGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATgg                                                    >  1:3261773/1‑71 (MQ=255)             aGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATgg                                                    >  1:4085528/1‑71 (MQ=255)              gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATg                                                     <  1:3149440/69‑1 (MQ=255)              gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                                   <  1:455712/71‑1 (MQ=255)                 tAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAtt                                                                     <  1:2840234/50‑1 (MQ=255)                 tAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                         <  1:3671411/62‑1 (MQ=255)                  aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                       >  1:3621247/1‑63 (MQ=255)                  aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTa                                                 <  1:265145/69‑1 (MQ=255)                  aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTa                                                 <  1:1081723/69‑1 (MQ=255)                   agagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGt                                              <  1:2388650/71‑1 (MQ=255)                    gagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTAc                                                >  1:1181721/1‑68 (MQ=255)                     agCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTa                                                 >  1:2069328/1‑66 (MQ=255)                     agCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTa                                            <  1:823000/71‑1 (MQ=255)                       cgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACt                                          <  1:1422282/71‑1 (MQ=255)                       cgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACt                                          <  1:2094511/71‑1 (MQ=255)                            cTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCa                                                        <  1:2606398/52‑1 (MQ=255)                             ttACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTa                                    >  1:2128142/1‑71 (MQ=255)                              tACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAAt                                                      >  1:3039580/1‑52 (MQ=255)                              tACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTa                                    <  1:1388958/70‑1 (MQ=255)                              tACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTAc                                   >  1:2474137/1‑71 (MQ=255)                               aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTc                                                                    >  1:2420654/1‑37 (MQ=255)                               aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTc                                                                    >  1:3689858/1‑37 (MQ=255)                               aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                       >  1:3830835/1‑50 (MQ=255)                               aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                       >  1:1842152/1‑50 (MQ=255)                               aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCa                                                        >  1:1364942/1‑49 (MQ=255)                               aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                                  >  1:2414018/1‑71 (MQ=255)                               aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                                  >  1:1015145/1‑71 (MQ=255)                                cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTg                                 <  1:160589/71‑1 (MQ=255)                                cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTg                                 <  1:699814/71‑1 (MQ=255)                                cTGAAGCAAAATATAACC‑TGAGCAGACTGTCATTCgcgcgc                                                              >  1:3298477/1‑41 (MQ=255)                                 tGAAGCAAAATATAACCTTGAGCAGACTGTCATTcg                                                                   <  1:1273279/36‑1 (MQ=255)                                   aaGCAAAATATAACCTTGAGCAGACTGTCATTCACGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                <  1:478671/69‑1 (MQ=255)                                     gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTc                                         <  1:3472115/58‑1 (MQ=255)                                     gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCg                            >  1:3770683/1‑71 (MQ=255)                                     gCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGcga                                                       >  1:2898008/1‑42 (MQ=16)                                      cAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGc                           <  1:607265/71‑1 (MQ=255)                                       aaaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                          <  1:1629912/71‑1 (MQ=255)                                        aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                          >  1:3217210/1‑70 (MQ=255)                                        aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                          >  1:3893099/1‑70 (MQ=255)                                        aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                          >  1:1231582/1‑70 (MQ=255)                                         aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGtga                                            <  1:2039935/51‑3 (MQ=255)                                          atatAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                       >  1:884419/1‑39 (MQ=255)                                          atatAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                       >  1:3573855/1‑71 (MQ=255)                                           tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                       >  1:746901/1‑70 (MQ=255)                                           tataACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCggg                      >  1:90012/1‑71 (MQ=255)                                            ataACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCggg                      >  1:4124072/1‑70 (MQ=255)                                             taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATgg                                                    >  1:96724/1‑39 (MQ=255)                                             taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGccc                         <  1:24079/66‑1 (MQ=255)                                             taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                    <  1:125495/71‑1 (MQ=255)                                             taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                    >  1:145549/1‑71 (MQ=255)                                             taACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGt                     <  1:647673/70‑1 (MQ=255)                                             taACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGt                     <  1:173734/70‑1 (MQ=255)                                               aCCTTGAGCAGACTGTCATTCGCGCTCCGAGCAATGGCTACGTTACTCAGGTACTGATCC‑CCCAGTTACa                  <  1:1564847/70‑1 (MQ=255)                                               aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                  <  1:369893/71‑1 (MQ=255)                                               aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                  <  1:364936/71‑1 (MQ=255)                                                  ttGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTa                                            >  1:52763/1‑42 (MQ=255)                                                  ttGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                   >  1:91420/1‑67 (MQ=255)                                                          aCTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCa                        <  1:352150/54‑1 (MQ=255)                                                          aCTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGc       <  1:3178123/71‑1 (MQ=255)                                                             gTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                                <  1:2200964/43‑1 (MQ=255)                                                              tCATTCGCTCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCtt    <  1:2425951/70‑1 (MQ=255)                                                              tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTg                                 <  1:50974/41‑1 (MQ=255)                                                              tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTg                                 <  1:1440604/41‑1 (MQ=255)                                                              tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATc                              >  1:1157155/1‑44 (MQ=255)                                                               cATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATc                              >  1:781991/1‑43 (MQ=255)                                                               cATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGc  >  1:1136439/1‑71 (MQ=255)                                                                aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATcccc                           <  1:1082063/45‑3 (MQ=255)                                                                aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGccc                         <  1:2686883/47‑1 (MQ=255)                                                                 ttCGCGCGCCGAGCAATGGCTACGTTACTCAGTTACTGATCCGCCCAGGTa                    <  1:2577331/51‑1 (MQ=255)                                                                 ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                          <  1:3357757/45‑1 (MQ=255)                                                                 ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg        >  1:893123/1‑63 (MQ=255)                                                                  |                                                                    TCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGC  >  minE/2570696‑2570828

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