Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,876	C→T	27.4%	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	27.4%	124.4 / 54.2	73	noncoding (1219/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (38/15);  new base T (10/10);  total (48/25)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCG  >  minE/2359818‑2359942                                                           |                                                                   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                        >  1:24964/1‑71 (MQ=34) gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                        >  1:1412868/1‑71 (MQ=255) gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                        <  1:66456/71‑1 (MQ=255)  cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                       >  1:469020/1‑71 (MQ=255)    ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                     >  1:1976749/1‑71 (MQ=255)     tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                     <  1:91623/70‑1 (MQ=35)     tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                      <  1:1054478/69‑1 (MQ=255)      aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                   >  1:1897381/1‑71 (MQ=34)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                   <  1:1976499/68‑1 (MQ=35)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                                <  1:458686/71‑1 (MQ=255)          ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                               <  1:450087/71‑1 (MQ=255)           tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                   >  1:502575/1‑66 (MQ=35)           tACATATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                           >  1:148121/1‑58 (MQ=9)            aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                              >  1:742612/1‑70 (MQ=255)              ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                              >  1:878157/1‑68 (MQ=35)              ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                              >  1:963883/1‑68 (MQ=35)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                                  >  1:995360/1‑62 (MQ=35)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                           <  1:95/69‑1 (MQ=255)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                          >  1:1343740/1‑70 (MQ=255)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                         >  1:617124/1‑71 (MQ=255)                      gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                          <  1:383399/64‑1 (MQ=11)                      gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                                   <  1:1631892/71‑1 (MQ=255)                       cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                                      >  1:1598157/1‑67 (MQ=35)                       cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                  <  1:117699/71‑1 (MQ=34)                         ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                 >  1:137019/1‑70 (MQ=35)                         ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGc                                           <  1:1477282/60‑1 (MQ=11)                         ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                                 >  1:300083/1‑70 (MQ=255)                           cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                           >  1:1403569/1‑42 (MQ=32)                           cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                       >  1:93021/1‑62 (MQ=35)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:844030/1‑70 (MQ=255)                            gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                     <  1:779038/47‑1 (MQ=38)                             cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                            <  1:1940499/71‑1 (MQ=34)                              aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                         >  1:1672312/1‑57 (MQ=35)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                          >  1:1596312/1‑71 (MQ=35)                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:1823441/1‑65 (MQ=255)                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                          >  1:161980/1‑37 (MQ=37)                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                          >  1:494932/1‑37 (MQ=37)                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTcac                                                          >  1:444005/1‑37 (MQ=37)                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                           >  1:807939/1‑36 (MQ=33)                                 tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                        >  1:553395/1‑71 (MQ=39)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                        >  1:1729478/1‑38 (MQ=40)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                        <  1:1574592/70‑1 (MQ=40)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                        >  1:1679362/1‑70 (MQ=40)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                        >  1:1799976/1‑70 (MQ=40)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                       <  1:1512889/68‑1 (MQ=255)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt                    >  1:366999/1‑71 (MQ=255)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt                    <  1:1340718/71‑1 (MQ=255)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg                     >  1:353240/1‑70 (MQ=255)                                       aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                  >  1:343933/1‑71 (MQ=255)                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc               >  1:1302178/1‑71 (MQ=255)                                           ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt              >  1:392370/1‑71 (MQ=255)                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:524157/1‑53 (MQ=35)                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                          >  1:520875/1‑58 (MQ=38)                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg             >  1:1198048/1‑71 (MQ=255)                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg             >  1:1461840/1‑71 (MQ=255)                                            aaGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                               >  1:514411/2‑37 (MQ=255)                                              gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                            <  1:2050315/38‑1 (MQ=40)                                              gCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                       <  1:80719/43‑1 (MQ=14)                                               cATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                              >  1:879528/1‑35 (MQ=40)                                               cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg          >  1:1261977/1‑71 (MQ=255)                                                 tGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCa        <  1:1081925/71‑1 (MQ=255)                                                    cTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                           <  1:2088661/49‑1 (MQ=255)                                                    cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                          <  1:1620731/50‑1 (MQ=40)                                                    cTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCt     <  1:1008048/71‑1 (MQ=255)                                                      cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc   >  1:668492/1‑71 (MQ=255)                                                       acaGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                               >  1:1038217/1‑42 (MQ=255)                                                       acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                        >  1:433677/1‑49 (MQ=255)                                                       acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc   <  1:987269/70‑1 (MQ=255)                                                       acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg  <  1:1783151/71‑1 (MQ=255)                                                       acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg  <  1:789733/71‑1 (MQ=255)                                                       acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg  >  1:1482366/1‑71 (MQ=255)                                                       acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg  >  1:6844/1‑71 (MQ=255)                                                       acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg  >  1:662137/1‑71 (MQ=255)                                                           |                                                                   GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCG  >  minE/2359818‑2359942

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