Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,556	T→C	30.7%	noncoding (288/1542 nt)	rrsG ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,763,556	0	T	C	30.7%	118.0 / 78.7	88	noncoding (288/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (23/38);  new base C (11/16);  total (34/54)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  minE/1763502‑1763624                                                       |                                                                     gtaGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      <  1:249673/69‑1 (MQ=17) ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCtt                                                       <  1:947523/70‑1 (MQ=31) ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                      >  1:1070446/1‑71 (MQ=17) ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                      <  1:488761/71‑1 (MQ=17) ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                      <  1:1463072/71‑1 (MQ=17)   aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                    >  1:30238/1‑71 (MQ=21)     tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGa                                                  <  1:2968787/71‑1 (MQ=34)     tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGa                                                  <  1:131553/71‑1 (MQ=34)     tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGa                                                  <  1:24048/71‑1 (MQ=34)     tCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGt                                                    >  1:85136/1‑69 (MQ=25)       tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                >  1:971596/1‑71 (MQ=34)       tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                >  1:276426/1‑71 (MQ=21)       tGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGc                                                >  1:2277082/1‑71 (MQ=21)         gACCGTGTCTCAGTTCCAGTGTGGCTTGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAg                                                 <  1:1132258/68‑1 (MQ=25)         gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGc                                                <  1:832214/69‑1 (MQ=35)          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:1118318/71‑1 (MQ=34)          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:2263023/71‑1 (MQ=34)          aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:196189/71‑1 (MQ=34)            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      <  1:2994399/60‑1 (MQ=33)            cGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           <  1:1618049/71‑1 (MQ=21)             gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          <  1:1398828/71‑1 (MQ=34)             gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          <  1:2256201/71‑1 (MQ=34)             gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAc                                          <  1:2346561/71‑1 (MQ=21)              tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         >  1:2700916/1‑71 (MQ=34)               gtCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGa                                                  >  1:759790/1‑61 (MQ=11)               gtCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGcc                                               >  1:1348522/1‑64 (MQ=11)               gtCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGcc                                               >  1:1630127/1‑64 (MQ=11)                tctcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTg                                                      <  1:2471919/56‑1 (MQ=11)                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcccc                                       <  1:2553431/71‑1 (MQ=34)                tctcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcccc                                       <  1:459712/71‑1 (MQ=34)                  tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                     >  1:2147685/1‑71 (MQ=34)                  tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                     >  1:845958/1‑71 (MQ=34)                  tcAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCcca                                      >  1:2525004/1‑70 (MQ=35)                    aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                    >  1:2425376/1‑70 (MQ=35)                     gTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                      <  1:2626603/51‑1 (MQ=16)                      ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 <  1:1237903/71‑1 (MQ=39)                      ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                 <  1:1540399/71‑1 (MQ=39)                      ttCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:700277/59‑1 (MQ=255)                        ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaaca                                <  1:583570/70‑1 (MQ=255)                         cAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAg                              <  1:2577039/71‑1 (MQ=255)                           gtgtGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                            <  1:3095126/71‑1 (MQ=255)                            tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATGGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                           <  1:1259553/71‑1 (MQ=255)                            tgtgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                           <  1:1292392/71‑1 (MQ=255)                              tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                      <  1:960019/42‑1 (MQ=35)                                gCTGGTCATCCTCTCAGTCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATc                        <  1:3014807/70‑1 (MQ=255)                                gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    >  1:898966/1‑42 (MQ=38)                                gCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    >  1:2844678/1‑42 (MQ=38)                                  tGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATccc                      <  1:360200/70‑1 (MQ=255)                                   ggTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     >  1:3297030/1‑70 (MQ=255)                                    gTCTTCCTCTCAGACCAGCTAGGGATCGTCGCCTTTGTGAGCCGTTACCCCACCAACAAGCTa                           <  1:2590620/63‑1 (MQ=255)                                    gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     <  1:2713011/69‑1 (MQ=255)                                    gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     <  1:2659870/69‑1 (MQ=255)                                    gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     <  1:2485532/69‑1 (MQ=255)                                    gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                   <  1:3262853/71‑1 (MQ=255)                                    gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             >  1:1863440/1‑45 (MQ=255)                                    gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             >  1:1350009/1‑45 (MQ=255)                                    gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     <  1:1689289/69‑1 (MQ=255)                                    gTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     <  1:1979704/69‑1 (MQ=255)                                     tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt                  >  1:3095012/1‑71 (MQ=255)                                      cgtCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTGAGCCTTTa                                           >  1:1002807/3‑45 (MQ=11)                                       aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                    >  1:239725/1‑35 (MQ=40)                                       aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAACCTAATCCCATCTgg                >  1:3336392/1‑71 (MQ=255)                                       aTCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAccc                                        >  1:2367096/1‑47 (MQ=255)                                       aTCCTCTCAGACAAGCCAGGGATCGTCGCCTTGGTGAGCCGTTa                                           >  1:3115214/1‑44 (MQ=11)                                        tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:953631/41‑1 (MQ=255)                                        tCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                             <  1:2804409/41‑1 (MQ=255)                                         ccTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               >  1:1199714/1‑70 (MQ=255)                                            ctcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                         >  1:1419692/1‑41 (MQ=38)                                             tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcac            <  1:170357/69‑1 (MQ=255)                                             tcAGACCAGCCAAGGATCGTCCCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCa                     <  1:2451334/60‑1 (MQ=17)                                              cAGACCAGCTAGGGATCGTCGCCTTGGTGCGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc         >  1:1918379/1‑71 (MQ=255)                                              cAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc         >  1:503159/1‑71 (MQ=255)                                               aGACCAGCTAGGGATCGTCGCTTTGGTGAGCCGTTACCCCACCAACAAGc                             >  1:1597533/1‑50 (MQ=255)                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca           <  1:951459/68‑1 (MQ=255)                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:2581829/71‑1 (MQ=255)                                               aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:327474/71‑1 (MQ=255)                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:3335757/71‑1 (MQ=255)                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        >  1:2682884/1‑71 (MQ=255)                                               aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc        <  1:1848064/71‑1 (MQ=255)                                                gACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCg       <  1:2300424/71‑1 (MQ=255)                                                 aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa      >  1:3011801/1‑71 (MQ=255)                                                 aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                       <  1:1862634/54‑1 (MQ=255)                                                 aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATcc                       <  1:1698830/54‑1 (MQ=255)                                                   cAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa      >  1:2996018/1‑69 (MQ=255)                                                   cAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa      >  1:603063/1‑69 (MQ=255)                                                     gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc  <  1:950078/71‑1 (MQ=255)                                                     gCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGc  <  1:349756/71‑1 (MQ=255)                                                      cTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTggg               <  1:3107819/57‑1 (MQ=255)                                                       |                                                                     GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGC  >  minE/1763502‑1763624

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