Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	46.1%	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	46.1%	52.6 / 109.3	78	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (24/18);  new base A (18/18);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCG  >  minE/2359818‑2359934                                                  |                                                                    gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                >  1:2877986/1‑71 (MQ=34) gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                >  1:2068957/1‑71 (MQ=34) gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                <  1:139878/71‑1 (MQ=34)   tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                               >  1:306494/1‑70 (MQ=255)     tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             >  1:757652/1‑70 (MQ=35)      aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                             <  1:2392173/69‑1 (MQ=255)       tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                           <  1:2522956/70‑1 (MQ=35)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                         <  1:2011619/70‑1 (MQ=35)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                        <  1:1446342/71‑1 (MQ=34)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                        <  1:244078/71‑1 (MQ=34)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                         <  1:1234611/70‑1 (MQ=255)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                        <  1:2578033/71‑1 (MQ=255)          ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                        >  1:1860133/1‑70 (MQ=255)          ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                       <  1:1320649/71‑1 (MQ=255)           tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                           >  1:1908192/1‑66 (MQ=255)           tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                           >  1:1770748/1‑66 (MQ=255)            aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                      <  1:1645731/70‑1 (MQ=35)            aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      >  1:1304051/1‑70 (MQ=255)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   >  1:927730/1‑53 (MQ=31)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTTACAGTACACCTTCACAGGCTTACAGAACGc                                   <  1:309652/69‑1 (MQ=11)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                >  1:1086804/1‑56 (MQ=11)                aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                             >  1:1642377/1‑59 (MQ=11)                aTGTCAGCATTCCCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                 >  1:305887/1‑71 (MQ=11)                 tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                >  1:2438209/1‑71 (MQ=255)                   tCAGCATTCGCACTTCTGATACCTCCAGCATGCGTCACAGCACACCTTCACAGGCTTACAGa                                       <  1:2934585/62‑1 (MQ=25)                   tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                              <  1:1281251/71‑1 (MQ=255)                    cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                      >  1:1102336/1‑62 (MQ=35)                    cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                              >  1:720535/1‑70 (MQ=35)                     aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                             <  1:1667240/70‑1 (MQ=255)                     aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                             <  1:2878774/70‑1 (MQ=255)                     aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                             <  1:1380882/70‑1 (MQ=255)                      gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                           <  1:2199996/71‑1 (MQ=255)                      gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                           >  1:1961986/1‑71 (MQ=255)                       cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcaca                                                 >  1:586703/1‑48 (MQ=38)                       cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                          >  1:2725492/1‑55 (MQ=35)                       cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:211403/1‑71 (MQ=34)                       cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                   >  1:497636/1‑46 (MQ=11)                       cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:308094/1‑71 (MQ=255)                       cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:2931378/1‑71 (MQ=255)                        aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                         <  1:1653236/71‑1 (MQ=34)                         ttCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                <  1:1139034/47‑1 (MQ=11)                           cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                           <  1:1548328/50‑1 (MQ=38)                           cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                       >  1:724396/1‑70 (MQ=35)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                           >  1:528784/1‑50 (MQ=11)                            gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCtt                                            <  1:1225384/48‑1 (MQ=11)                              aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                           <  1:1162471/47‑1 (MQ=38)                               cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      <  1:33267/51‑1 (MQ=11)                               cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                  >  1:1881897/1‑70 (MQ=255)                                ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                      <  1:1713944/50‑1 (MQ=38)                                ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                  >  1:1526864/1‑70 (MQ=37)                                ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                  >  1:1903302/1‑70 (MQ=37)                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                         <  1:426962/47‑1 (MQ=11)                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                 <  1:1432454/55‑1 (MQ=11)                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                 <  1:108771/55‑1 (MQ=11)                                 tCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacgc                 >  1:1826267/1‑68 (MQ=11)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                >  1:799297/1‑70 (MQ=40)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                >  1:920005/1‑70 (MQ=40)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                >  1:1374879/1‑70 (MQ=40)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                >  1:2169013/1‑70 (MQ=40)                                  cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacgc                 <  1:826395/69‑3 (MQ=11)                                   tGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACg                                    <  1:1927841/49‑1 (MQ=11)                                    gATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                             <  1:2790287/39‑1 (MQ=40)                                    gATACCTCCAGCATGCCTCACAGCAAACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAg             >  1:313576/1‑71 (MQ=255)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt               <  1:2048371/68‑1 (MQ=40)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt               <  1:2045103/68‑1 (MQ=40)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt            <  1:1723937/71‑1 (MQ=255)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt            <  1:398862/71‑1 (MQ=255)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt            <  1:1844809/71‑1 (MQ=255)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt            >  1:1085545/1‑71 (MQ=255)                                       aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt          >  1:1194434/1‑71 (MQ=255)                                       aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      >  1:630096/1‑43 (MQ=14)                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc       >  1:163174/1‑71 (MQ=255)                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc       >  1:1619741/1‑71 (MQ=255)                                          tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      >  1:2862434/1‑40 (MQ=14)                                             aGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                         >  1:835440/1‑50 (MQ=38)                                             aGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGc    <  1:2711943/71‑1 (MQ=255)                                               cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg  >  1:411978/1‑71 (MQ=255)                                               cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg  >  1:1735442/1‑71 (MQ=255)                                                  |                                                                    GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCG  >  minE/2359818‑2359934

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