Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	34.1%	noncoding (1219/2905 nt)	rrlA ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,359,876	0	C	T	34.1%	113.5 / 82.8	85	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (35/21);  new base T (21/8);  total (56/29)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCG  >  minE/2359815‑2359942                                                              |                                                                   cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                           <  1:2614026/71‑1 (MQ=31)     cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                                       >  1:982960/1‑71 (MQ=34)     cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                                       >  1:1724143/1‑71 (MQ=34)     cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                                       <  1:1572449/71‑1 (MQ=34)     cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                       <  1:3142201/71‑1 (MQ=255)     cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                       <  1:346414/71‑1 (MQ=255)      tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                      >  1:2479960/1‑71 (MQ=255)       ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                     >  1:33901/1‑71 (MQ=34)           cGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                   >  1:1429680/1‑69 (MQ=255)           cGTTACTTATGTCAGCATTCGCACTTCTGATACCTACAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                   >  1:1493397/1‑69 (MQ=25)            gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                <  1:1044328/71‑1 (MQ=34)            gTTACTTATGTCAGCATTAGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                                   >  1:712985/1‑68 (MQ=11)               aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAAc                                             >  1:3155421/1‑71 (MQ=34)               aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                              <  1:1831006/70‑1 (MQ=255)               aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                              <  1:3279495/70‑1 (MQ=255)               aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                              <  1:3275207/70‑1 (MQ=255)               aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAAc                                             >  1:682856/1‑71 (MQ=255)                     gTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                       >  1:3267138/1‑71 (MQ=255)                        aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                     <  1:949668/70‑1 (MQ=35)                        aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                     <  1:945720/70‑1 (MQ=35)                        aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                    >  1:2036231/1‑71 (MQ=34)                        aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                                    >  1:1046479/1‑71 (MQ=34)                        aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                         >  1:782118/1‑66 (MQ=11)                        aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                         >  1:2731854/1‑66 (MQ=11)                             tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                               >  1:57780/1‑55 (MQ=11)                             tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:1340370/1‑71 (MQ=255)                             tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:609476/1‑71 (MQ=255)                              cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:260759/1‑70 (MQ=35)                              cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                              >  1:1683142/1‑55 (MQ=11)                              cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                              >  1:1748464/1‑55 (MQ=11)                              cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                     >  1:3079029/1‑64 (MQ=255)                              cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                     >  1:284289/1‑64 (MQ=255)                              cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                     >  1:1110881/1‑64 (MQ=255)                              agCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                              >  1:961537/2‑55 (MQ=11)                               gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                              <  1:22728/70‑1 (MQ=35)                               gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                              >  1:1238462/1‑70 (MQ=35)                               gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                              >  1:477067/1‑70 (MQ=35)                                 aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:2558054/1‑67 (MQ=35)                                  cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                          >  1:2523233/1‑70 (MQ=255)                                  cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                           >  1:912591/1‑70 (MQ=255)                                     cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                          >  1:2946382/1‑68 (MQ=37)                                     cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                        >  1:1851500/1‑70 (MQ=40)                                     cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                        >  1:1189157/1‑70 (MQ=40)                                     cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAAAGAACGCTCCCCTACCCAACAacac                         >  1:1829806/1‑69 (MQ=255)                                        aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                     <  1:2547174/54‑1 (MQ=35)                                        aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                     <  1:1734987/54‑1 (MQ=35)                                        aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                     <  1:1948075/54‑1 (MQ=35)                                        aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                     <  1:2386666/54‑1 (MQ=35)                                        aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                     <  1:1430701/54‑1 (MQ=35)                                              ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACATCACATAGTGTCGc               >  1:2137463/1‑70 (MQ=255)                                              ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt              >  1:924656/1‑71 (MQ=255)                                              ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt              >  1:2090222/1‑71 (MQ=255)                                               cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                 >  1:676131/1‑67 (MQ=255)                                               cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc                 >  1:3188480/1‑67 (MQ=255)                                               cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg             <  1:2969483/71‑1 (MQ=255)                                               cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg             <  1:2607136/71‑1 (MQ=255)                                                 gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                         <  1:1440657/41‑1 (MQ=38)                                                 gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc           >  1:968100/1‑71 (MQ=255)                                                 gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc           >  1:2328205/1‑71 (MQ=255)                                                 gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc            >  1:1392526/1‑70 (MQ=255)                                                 gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc            >  1:966382/1‑70 (MQ=255)                                                 gCATACCTCACAGTACACCTTCACAGGCTTACAGa                                               <  1:1017131/35‑1 (MQ=14)                                                 gCATACCTCACAGTACACCTTCACAGGCTTACAGa                                               <  1:2688902/35‑1 (MQ=14)                                                 gCATACCTCACAGTACACCTTCACAGGCTTACAGa                                               <  1:1001653/35‑1 (MQ=14)                                                  cATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                                      >  1:2564914/1‑43 (MQ=14)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                              <  1:3447536/49‑1 (MQ=38)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                              <  1:3200085/49‑1 (MQ=38)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                              <  1:55083/49‑1 (MQ=38)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                              <  1:3023979/49‑1 (MQ=38)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                              <  1:2543887/49‑1 (MQ=38)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                              <  1:612461/49‑1 (MQ=38)                                                      ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                            >  1:3448545/1‑49 (MQ=38)                                                      ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                            >  1:2466695/1‑49 (MQ=38)                                                      ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                            >  1:2032944/1‑49 (MQ=38)                                                      ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                            >  1:1078912/1‑49 (MQ=38)                                                      ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                            >  1:1990233/1‑49 (MQ=38)                                                      ccTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                            >  1:741658/1‑49 (MQ=38)                                                       cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCt     <  1:3455054/71‑1 (MQ=255)                                                         cacaTCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCtt    >  1:1324367/1‑70 (MQ=255)                                                         cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                     >  1:604045/1‑37 (MQ=40)                                                         cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc   >  1:2929096/1‑71 (MQ=255)                                                          acaGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:2666068/1‑42 (MQ=255)                                                          acaGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:3352599/1‑42 (MQ=255)                                                          acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg  >  1:1215733/1‑71 (MQ=255)                                                          acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg  >  1:2411418/1‑71 (MQ=255)                                                              |                                                                   CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCG  >  minE/2359815‑2359942

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