Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	40.2%	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	40.2%	71.2 / 94.7	92	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base A (21/16);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTG  >  minE/2359818‑2359931                                                  |                                                                 gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                             <  1:4314411/71‑1 (MQ=34) gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                             >  1:3216557/1‑71 (MQ=255) gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                             >  1:2000973/1‑71 (MQ=255) gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                             >  1:4602559/1‑71 (MQ=255)  cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                <  1:1836757/67‑1 (MQ=31)  cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                <  1:5153493/67‑1 (MQ=31)  cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                            >  1:331885/1‑71 (MQ=34)   tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                           >  1:5623932/1‑71 (MQ=34)   tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                           >  1:4465113/1‑71 (MQ=255)    ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                          >  1:1573809/1‑71 (MQ=255)     tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                          >  1:4457524/1‑70 (MQ=255)       tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                        >  1:4141641/1‑70 (MQ=35)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                      <  1:1137304/70‑1 (MQ=35)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     <  1:1456590/71‑1 (MQ=34)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     <  1:230144/71‑1 (MQ=34)         gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                     <  1:1034871/71‑1 (MQ=255)          ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     >  1:2651374/1‑70 (MQ=35)          ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     >  1:2527219/1‑70 (MQ=35)          ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     >  1:5405943/1‑70 (MQ=35)          ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCtt                                         <  1:4459936/66‑1 (MQ=255)             cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 <  1:468211/71‑1 (MQ=34)               tATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                               >  1:3215094/1‑71 (MQ=255)                 tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 <  1:1179536/67‑1 (MQ=35)                 tGTCAGCATTCGCACTTCTGATACCTCCAGCATACGTCACAGTACACCTTCACAGGCTTACAGAACGCTc                              <  1:3012951/70‑1 (MQ=11)                 tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                              <  1:3318125/70‑1 (MQ=255)                 tGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                              <  1:1450232/70‑1 (MQ=255)                  gTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                            >  1:474960/1‑71 (MQ=34)                   tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                        <  1:3607846/58‑1 (MQ=35)                   tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                           <  1:5434838/71‑1 (MQ=255)                    cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 >  1:4955679/1‑64 (MQ=35)                    cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 >  1:1059741/1‑64 (MQ=35)                    cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 >  1:5394983/1‑64 (MQ=35)                    cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                 >  1:786363/1‑64 (MQ=35)                     aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                          <  1:1678427/70‑1 (MQ=35)                     aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCt                          <  1:3078047/70‑1 (MQ=35)                     aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                          <  1:4822968/70‑1 (MQ=255)                     aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTAACAGGCTTACAGAACGCTCCCCTa                         <  1:1551669/71‑1 (MQ=11)                      gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                     <  1:3139555/58‑1 (MQ=11)                       cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                          >  1:1774057/1‑68 (MQ=255)                         ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                        >  1:2011287/1‑68 (MQ=35)                           cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                        >  1:2394865/1‑50 (MQ=38)                           cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                    >  1:3728057/1‑70 (MQ=35)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                   >  1:5267017/1‑55 (MQ=11)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                   >  1:1050604/1‑55 (MQ=11)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                   >  1:5664140/1‑55 (MQ=11)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                   >  1:2862491/1‑55 (MQ=11)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGACCGCTCCCCTACCcaa                    >  1:2506401/1‑70 (MQ=11)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                    >  1:3464800/1‑70 (MQ=255)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                    >  1:5597559/1‑70 (MQ=255)                           cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                    >  1:5040003/1‑70 (MQ=255)                            gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaaca                  <  1:2474110/71‑1 (MQ=255)                            gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                   <  1:1425669/70‑1 (MQ=255)                             cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                 <  1:5349487/71‑1 (MQ=34)                             cACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                     <  1:3668574/67‑1 (MQ=255)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     <  1:1491225/49‑1 (MQ=38)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     <  1:4719608/49‑1 (MQ=38)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                     <  1:2378503/49‑1 (MQ=38)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                >  1:2604141/1‑70 (MQ=35)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               <  1:5162549/71‑1 (MQ=35)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               <  1:1655642/71‑1 (MQ=35)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               >  1:720663/1‑71 (MQ=35)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               <  1:1923196/71‑1 (MQ=35)                               cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               <  1:271100/71‑1 (MQ=35)                               cTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg               <  1:3593738/71‑2 (MQ=255)                                ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                       <  1:3222575/46‑1 (MQ=38)                                ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcccc                           <  1:4416140/58‑1 (MQ=35)                                ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                              <  1:1190031/55‑1 (MQ=35)                                ttCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca               >  1:5214529/1‑70 (MQ=37)                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                     <  1:3713487/48‑1 (MQ=11)                                ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                <  1:1609437/69‑1 (MQ=255)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                              <  1:1627755/53‑1 (MQ=35)                                  cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca             >  1:361626/1‑70 (MQ=40)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt            <  1:1393065/68‑1 (MQ=40)                                     aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt            <  1:1131736/68‑1 (MQ=40)                                     aTACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                   <  1:1378131/61‑1 (MQ=255)                                       aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                 >  1:2301972/1‑61 (MQ=255)                                       aCCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                 >  1:827228/1‑61 (MQ=255)                                        ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc      <  1:50855/71‑1 (MQ=255)                                         cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg     <  1:4109123/71‑1 (MQ=255)                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                <  1:3853896/59‑1 (MQ=35)                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                <  1:165367/59‑1 (MQ=35)                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                <  1:1947825/59‑1 (MQ=35)                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaac                <  1:837182/59‑1 (MQ=35)                                          tCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc    >  1:170487/1‑71 (MQ=255)                                          tCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                        >  1:5323336/1‑35 (MQ=37)                                          tCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                        >  1:2651291/1‑35 (MQ=37)                                          tCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                        >  1:4084914/1‑35 (MQ=37)                                           ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca             >  1:4063043/1‑61 (MQ=40)                                           ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca             >  1:5630974/1‑61 (MQ=40)                                           ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt   >  1:4547525/1‑71 (MQ=255)                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg  >  1:3524338/1‑71 (MQ=255)                                            cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg  >  1:2333166/1‑71 (MQ=255)                                                  |                                                                 GCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTG  >  minE/2359818‑2359931

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