Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,758	C→A	7.2%	V221V (GTC→GTA) ‡	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,758	0	C	A	7.2%	278.1 / 4.6	83	V221V (GTC→GTA) ‡	yibH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/39);  new base A (3/3);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

CGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGC  >  minE/2570689‑2570828                                                                      |                                                                       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: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)              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:4065684/70‑1 (MQ=255)                tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                         <  1:1837824/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:1081723/69‑1 (MQ=255)                         aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTa                                                 <  1:265145/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:2094511/71‑1 (MQ=255)                              cgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACt                                          <  1:1422282/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: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: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:173734/70‑1 (MQ=255)                                                    taACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGt                     <  1:647673/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)                                                                      |                                                                       CGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGC  >  minE/2570689‑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: ‑