Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,570,764	C→T	8.0%	R223R (CGC→CGT)	yibH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,570,764	0	C	T	8.0%	235.2 / 8.0	75	R223R (CGC→CGT)	yibH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base T (3/3);  total (34/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 = 8.62e-01

TCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCT  >  minE/2570696‑2570832                                                                     |                                                                     tCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                                    >  1:3058713/1‑71 (MQ=255)   gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                                  >  1:23045/1‑71 (MQ=255)   gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                                  >  1:2690411/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:1837824/70‑1 (MQ=255)         tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                             <  1:4065684/70‑1 (MQ=255)         tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                             <  1:3676153/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)                 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)                            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)                               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:699814/71‑1 (MQ=255)                                cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTg                                     <  1:160589/71‑1 (MQ=255)                                cTGAAGCAAAATATAACC‑TGAGCAGACTGTCATTCgcgcgc                                                                  >  1:3298477/1‑41 (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:3217210/1‑70 (MQ=255)                                        aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGcc                              >  1:3893099/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:173734/70‑1 (MQ=255)                                             taACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGt                         <  1:647673/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)                                                             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)                                                                   cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGa                                    <  1:1160454/37‑1 (MQ=255)                                                                   cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGa                                    <  1:1250970/37‑1 (MQ=255)                                                                   cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGa                                    <  1:4100317/37‑1 (MQ=255)                                                                   cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAg               >  1:1584515/1‑58 (MQ=255)                                                                   cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc              >  1:3002171/1‑59 (MQ=255)                                                                   cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCt  >  1:3876204/1‑71 (MQ=255)                                                                   cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCt  >  1:3700257/1‑71 (MQ=255)                                                                     |                                                                     TCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCT  >  minE/2570696‑2570832

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