Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,763,636	A→G	28.6%	noncoding (208/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,636	0	A	G	28.6%	97.0 / 50.9	70	noncoding (208/1542 nt)	rrsG	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base G (11/9);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTC  >  minE/1763569‑1763702                                                                    |                                                                   cTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGt                                                                 <  1:505915/71‑1 (MQ=255)  ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTc                                                                >  1:1695783/1‑71 (MQ=255)  ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTc                                                                >  1:1723624/1‑71 (MQ=255)  ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                                >  1:1220226/1‑71 (MQ=255)  ttGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTc                                                                >  1:1214099/1‑71 (MQ=255)   tGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTcc                                                               >  1:3300767/1‑71 (MQ=255)    ggTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTccc                                                              >  1:1577598/1‑71 (MQ=39)     gTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTc                                                                >  1:1547001/1‑68 (MQ=255)      tGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTccccc                                                            <  1:3177701/71‑1 (MQ=255)         gCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCctct                                                         >  1:2149591/1‑71 (MQ=39)          ccGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCtt                                                        <  1:3358638/71‑1 (MQ=255)           cGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCttt                                                       <  1:1525098/71‑1 (MQ=39)            gTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTg                                                      <  1:3031536/71‑1 (MQ=255)             ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                     <  1:1087258/71‑1 (MQ=39)             ttACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTgg                                                     <  1:3787567/71‑1 (MQ=39)              tACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGt                                                    >  1:4285383/1‑71 (MQ=255)               aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCttt                                                       >  1:4517999/1‑67 (MQ=255)               aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGt                                                    <  1:1110550/70‑1 (MQ=255)               aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGt                                                    >  1:1621961/1‑70 (MQ=255)               aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc                                                   >  1:4663740/1‑71 (MQ=255)               aCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTc                                                   >  1:832991/1‑71 (MQ=255)                ccccACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCt                                                  <  1:2425163/71‑1 (MQ=255)                ccccACCAACAAGCTAATCCCATCT‑GGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCt                                                  <  1:4291857/70‑1 (MQ=14)                   caccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                               >  1:1576669/1‑71 (MQ=39)                    accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGc                                               >  1:2639750/1‑70 (MQ=255)                    accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCg                                              <  1:966697/71‑1 (MQ=255)                    accaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGc                                               >  1:2150930/1‑70 (MQ=40)                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                             >  1:1872290/1‑71 (MQ=39)                     ccaACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGa                                             >  1:1381632/1‑71 (MQ=39)                     ccaACAAGCTAATCCCATCTGGGCACATCCGATCGCAAGAGGCCCGAAGGTCCCTCTCTTTGTTCTTGCGa                                             <  1:4049206/71‑1 (MQ=14)                       aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa                                                                    >  1:546695/1‑46 (MQ=40)                       aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCttt                                                       >  1:4387805/1‑59 (MQ=255)                       aacaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACa                                           >  1:2884910/1‑71 (MQ=38)                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAcgt                                          <  1:1039002/71‑3 (MQ=9)                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAGGGTCCCCCTCTTTGGTCTTGCGAcgt                                          <  1:1650901/71‑3 (MQ=9)                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCttt                                                       <  1:1633942/58‑1 (MQ=35)                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCtt                                                        <  1:1332741/57‑1 (MQ=35)                        acaaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACAt                                          <  1:2409129/71‑1 (MQ=38)                          aaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa                                        <  1:4137721/71‑1 (MQ=35)                          aaGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa                                        >  1:3614807/1‑71 (MQ=35)                                aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                   <  1:1063461/70‑1 (MQ=35)                                aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                   <  1:1623423/70‑1 (MQ=35)                                aTCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGt                                  <  1:3948916/71‑1 (MQ=34)                                   ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAtt                               <  1:1746430/71‑1 (MQ=34)                                   ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                <  1:980169/70‑1 (MQ=35)                                   ccATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAt                                >  1:3514238/1‑70 (MQ=35)                                   ccATCTCGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa                                        >  1:600767/1‑62 (MQ=25)                                    cATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCACTCTTTGGTCTTGCGACATTATGCGGTATTa                              <  1:2307381/71‑1 (MQ=21)                                        tGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTa                          >  1:2549556/1‑71 (MQ=34)                                            cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCgg                                   <  1:4411827/58‑1 (MQ=35)                                            cacaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGt                      <  1:4590587/71‑1 (MQ=34)                                             acaTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGtt                     >  1:2973031/1‑71 (MQ=34)                                               aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTa                                        <  1:646234/50‑1 (MQ=38)                                               aTCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                   <  1:1820788/71‑1 (MQ=34)                                                tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGc                            >  1:3193840/1‑61 (MQ=35)                                                tCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTc                   <  1:4464804/70‑1 (MQ=35)                                                 ccGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTa                              >  1:4059729/1‑58 (MQ=35)                                                  cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                <  1:3507602/71‑1 (MQ=34)                                                  cGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCag                >  1:1515507/1‑71 (MQ=34)                                                   gATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCTGTATTAGCTACCGTTTCCagt               <  1:1812767/71‑1 (MQ=21)                                                      ggCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtag             <  1:1607327/70‑1 (MQ=35)                                                        cAAGATGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCagtagt            <  1:4101261/69‑1 (MQ=25)                                                         aaGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTAt         >  1:1009649/1‑71 (MQ=34)                                                           gagGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcc       >  1:4149755/1‑71 (MQ=34)                                                              gCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTa          >  1:1653294/1‑65 (MQ=35)                                                              gCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTa          >  1:970810/1‑65 (MQ=35)                                                                ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGc                                     >  1:3617993/1‑36 (MQ=40)                                                                ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTa                                 >  1:2207488/1‑40 (MQ=40)                                                                ccGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTc  <  1:123297/71‑1 (MQ=34)                                                                 cGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTAtt                               <  1:2207065/41‑1 (MQ=38)                                                                  gAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATcccc     >  1:1780906/1‑66 (MQ=35)                                                                   aaGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTa          >  1:3429836/1‑60 (MQ=35)                                                                    |                                                                   CTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGGTCCCCCTCTTTGGTCTTGCGACATTATGCGGTATTAGCTACCGTTTCCAGTAGTTATCCCCCTC  >  minE/1763569‑1763702

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑