Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	38.0%	noncoding (1228/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,867	0	G	A	38.0%	73.1 / 78.8	71	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (24/20);  new base A (19/8);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCC  >  minE/2359815‑2359933                                                     |                                                                   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:3279495/70‑1 (MQ=255)               aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                     <  1:1831006/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:945720/70‑1 (MQ=35)                        aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:949668/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:609476/1‑71 (MQ=255)                             tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                      >  1:1340370/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:1110881/1‑64 (MQ=255)                              cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                            >  1:284289/1‑64 (MQ=255)                              cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                            >  1:3079029/1‑64 (MQ=255)                              agCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                     >  1:961537/2‑55 (MQ=11)                               gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     >  1:1238462/1‑70 (MQ=35)                               gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     >  1:477067/1‑70 (MQ=35)                               gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:22728/70‑1 (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:1948075/54‑1 (MQ=35)                                        aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:2386666/54‑1 (MQ=35)                                        aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:1734987/54‑1 (MQ=35)                                        aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                            <  1:1430701/54‑1 (MQ=35)                                              ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACATCACATAGTGTCGc      >  1:2137463/1‑70 (MQ=255)                                              ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt     >  1:2090222/1‑71 (MQ=255)                                              ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt     >  1:924656/1‑71 (MQ=255)                                               cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc        >  1:3188480/1‑67 (MQ=255)                                               cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc        >  1:676131/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:2328205/1‑71 (MQ=255)                                                 gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc  >  1:968100/1‑71 (MQ=255)                                                 gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc   >  1:1392526/1‑70 (MQ=255)                                                 gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc   >  1:966382/1‑70 (MQ=255)                                                 gCATACCTCACAGTACACCTTCACAGGCTTACAGa                                      <  1:1001653/35‑1 (MQ=14)                                                 gCATACCTCACAGTACACCTTCACAGGCTTACAGa                                      <  1:2688902/35‑1 (MQ=14)                                                 gCATACCTCACAGTACACCTTCACAGGCTTACAGa                                      <  1:1017131/35‑1 (MQ=14)                                                  cATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                             >  1:2564914/1‑43 (MQ=14)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:3023979/49‑1 (MQ=38)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:612461/49‑1 (MQ=38)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:55083/49‑1 (MQ=38)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:3447536/49‑1 (MQ=38)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:3200085/49‑1 (MQ=38)                                                    tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                     <  1:2543887/49‑1 (MQ=38)                                                     |                                                                   CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCC  >  minE/2359815‑2359933

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