Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,867	G→A	28.9%	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	28.9%	166.2 / 72.5	90	noncoding (1228/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (36/28);  new base A (15/11);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGC  >  minE/2359816‑2359935                                                    |                                                                     ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                    <  1:1219821/70‑1 (MQ=255)   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                 <  1:1967300/71‑1 (MQ=34)         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            >  1:1920024/1‑70 (MQ=255)           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                         <  1:1971451/71‑1 (MQ=34)           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAg                                         <  1:470092/71‑1 (MQ=255)            ttACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGa                                        <  1:2912684/71‑1 (MQ=34)             tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                       >  1:2429440/1‑71 (MQ=34)              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                          >  1:1947817/1‑67 (MQ=35)              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAAc                                      >  1:3236727/1‑71 (MQ=255)              aCTTATGTCAGAATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                          >  1:1806242/1‑67 (MQ=25)                ttATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                              <  1:1864834/61‑1 (MQ=255)                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                 <  1:1015999/71‑1 (MQ=34)                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                  <  1:1264742/70‑1 (MQ=35)                     tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                                <  1:945879/70‑1 (MQ=35)                     tCAGCATTCGCACTTCTGATACCTCCAGCATGCCTAACAGCACACCTTCACAGGCTTACAGAACGCTcccc                               <  1:2477947/71‑1 (MQ=21)                     tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcccc                               >  1:95498/1‑71 (MQ=255)                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                          >  1:599781/1‑58 (MQ=35)                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                          >  1:1049703/1‑58 (MQ=35)                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                             >  1:969507/1‑71 (MQ=34)                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                             >  1:1181780/1‑71 (MQ=34)                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                              <  1:1166160/70‑1 (MQ=255)                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                             >  1:2996662/1‑71 (MQ=255)                        gCATTCGCACTTCTGATACCTCC‑‑CATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                          <  1:3304396/71‑1 (MQ=11)                         cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                    >  1:2696427/1‑46 (MQ=32)                         cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                    >  1:1399405/1‑46 (MQ=32)                         cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                    >  1:2930220/1‑46 (MQ=32)                         cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                            >  1:2549612/1‑70 (MQ=255)                          aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                           >  1:3256678/1‑70 (MQ=35)                          aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                           >  1:205052/1‑70 (MQ=35)                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                          >  1:2221000/1‑70 (MQ=35)                           ttCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAccc                          >  1:2420337/1‑70 (MQ=35)                            tCGCACTTCTGATACCTCCAGCATGCCTCCCAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                        >  1:1482346/1‑71 (MQ=255)                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTcc                                 >  1:3148184/1‑61 (MQ=11)                             cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                        >  1:2638632/1‑70 (MQ=255)                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                            <  1:2095182/65‑1 (MQ=35)                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            <  1:331304/49‑1 (MQ=11)                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            <  1:342486/49‑1 (MQ=11)                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            <  1:107829/49‑1 (MQ=11)                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            <  1:1308939/49‑1 (MQ=11)                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            <  1:2027009/49‑1 (MQ=11)                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                            <  1:2103543/49‑1 (MQ=11)                               cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                     <  1:2316181/71‑1 (MQ=34)                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                     >  1:382897/1‑70 (MQ=35)                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                     >  1:1490019/1‑70 (MQ=35)                                 cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACacaacaaca                   >  1:1840053/1‑71 (MQ=25)                                 cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                   <  1:2189016/71‑1 (MQ=35)                                 cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                   <  1:917082/71‑1 (MQ=35)                                 cTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAacaaca                   <  1:1234766/71‑1 (MQ=35)                                   tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                 >  1:745389/1‑71 (MQ=39)                                   tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                 >  1:50013/1‑71 (MQ=39)                                   tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                 >  1:962524/1‑71 (MQ=39)                                   tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                 >  1:2481265/1‑71 (MQ=39)                                    cTGATACCTCCGGCATGCCTCACAGCACACCTTCACAGGCt                                              <  1:2216411/41‑1 (MQ=255)                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                              <  1:350818/41‑1 (MQ=38)                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                              <  1:2068936/41‑1 (MQ=38)                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                              <  1:2568887/41‑1 (MQ=38)                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                     >  1:1292132/1‑50 (MQ=38)                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                     >  1:577036/1‑50 (MQ=38)                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                       <  1:667999/64‑1 (MQ=35)                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac                  <  1:885586/69‑1 (MQ=39)                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                 <  1:3173601/70‑1 (MQ=40)                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                 <  1:861749/70‑1 (MQ=40)                                    cTGATACCTCCAGCATACCTCACAGTACACCTTCACAgg                                                >  1:2408508/1‑39 (MQ=14)                                    cTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                               >  1:1197054/1‑40 (MQ=14)                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACAt                <  1:1149781/68‑1 (MQ=40)                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt             >  1:672591/1‑71 (MQ=255)                                         aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt           >  1:1872661/1‑71 (MQ=255)                                         aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt           >  1:3360355/1‑71 (MQ=255)                                         aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt           >  1:1868908/1‑71 (MQ=255)                                         aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt           >  1:213144/1‑71 (MQ=255)                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt           <  1:2052295/70‑1 (MQ=255)                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc          <  1:3525/71‑1 (MQ=255)                                          ccTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTc          <  1:2469519/71‑1 (MQ=255)                                           cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg         <  1:3101752/71‑1 (MQ=255)                                            tCCAGCATGCCTAACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGc        >  1:1405040/1‑71 (MQ=255)                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                >  1:1183944/1‑47 (MQ=11)                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                >  1:3303068/1‑47 (MQ=11)                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                >  1:3177353/1‑47 (MQ=11)                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                                >  1:76169/1‑47 (MQ=11)                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                     >  1:159023/1‑58 (MQ=11)                                            tCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                     >  1:709563/1‑58 (MQ=11)                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                       >  1:811106/1‑38 (MQ=40)                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                       >  1:2083814/1‑38 (MQ=40)                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                       >  1:2801469/1‑38 (MQ=40)                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                       >  1:241427/1‑38 (MQ=40)                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg      >  1:905027/1‑71 (MQ=255)                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg      <  1:2336701/71‑1 (MQ=255)                                                 cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg   >  1:776948/1‑71 (MQ=255)                                                 cATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacgc                  <  1:64735/56‑3 (MQ=11)                                                  aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGc  <  1:1265842/71‑1 (MQ=255)                                                    |                                                                     CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGC  >  minE/2359816‑2359935

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