Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,755	T→G	7.3%	T220T (ACT→ACG)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,755	0	T	G	7.3%	208.0 / 12.4	82	T220T (ACT→ACG)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base G (3/3);  total (41/41)
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.53e-01

CGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGC  >  minE/2570686‑2570823                                                                      |                                                                     cGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTg                                                                     <  1:91890/71‑1 (MQ=255)   gCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGt                                                                    >  1:1712999/1‑70 (MQ=255)    cGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCa                                                                  <  1:2373341/71‑1 (MQ=255)     gAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAt                                                                 >  1:2562456/1‑71 (MQ=255)           tCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                           >  1:3058713/1‑71 (MQ=255)            cGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAtt                                                                >  1:1788102/1‑65 (MQ=255)             gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAt                                                                 >  1:3754501/1‑63 (MQ=255)             gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                         >  1:3660616/1‑71 (MQ=255)             gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                         >  1:2690411/1‑71 (MQ=255)             gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                         >  1:23045/1‑71 (MQ=255)                gATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGa                                                      <  1:2515169/71‑1 (MQ=255)                 aTTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCAt                                                                 <  1:372153/59‑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)                     tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGt                                                                    <  1:1062213/52‑1 (MQ=255)                     tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGt                                                                    <  1:1163982/52‑1 (MQ=255)                     tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                  <  1:3286176/70‑1 (MQ=255)                     tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAAt                                                 >  1:750207/1‑71 (MQ=255)                       aGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATgg                                               >  1:4085528/1‑71 (MQ=255)                       aGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATgg                                               >  1:3261773/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)                                      cTTACTGAAGCAAAATATAACCTTTAGCAGACTGTCa                                                                  <  1:3415530/37‑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:1015145/1‑71 (MQ=255)                                         aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACt                             >  1:2414018/1‑71 (MQ=255)                                          cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTg                            <  1:699814/71‑1 (MQ=255)                                          cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTg                            <  1:160589/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:1231582/1‑70 (MQ=255)                                                  aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                     >  1:3893099/1‑70 (MQ=255)                                                  aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                     >  1:3217210/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:145549/1‑71 (MQ=255)                                                       taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa               <  1:125495/71‑1 (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:364936/71‑1 (MQ=255)                                                         aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa             <  1:369893/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)                                                                      |                                                                     CGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGC  >  minE/2570686‑2570823

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