Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,225	G→A	28.1%	noncoding (1870/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,225	0	G	A	28.1%	115.9 / 37.9	92	noncoding (1870/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (31/35);  new base A (11/15);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGTATA  >  minE/2359158‑2359278                                                                    |                                                      tCGCTACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGc                                                      >  1:910568/1‑69 (MQ=39)  cGCTACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCg                                                     >  1:3429729/1‑69 (MQ=40)  cGCTACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCg                                                     >  1:2360182/1‑69 (MQ=40)     tACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGc                                                   >  1:1327380/1‑68 (MQ=255)     tACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCt                                                  >  1:341329/1‑69 (MQ=255)     tACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCt                                                  >  1:1056616/1‑69 (MQ=255)     tACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCt                                                  >  1:3156533/1‑69 (MQ=255)     tACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCt                                                  >  1:371375/1‑69 (MQ=255)     tACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCt                                                  >  1:2790557/1‑69 (MQ=255)     tACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCt                                                  <  1:3074728/69‑1 (MQ=255)         ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:1620996/68‑1 (MQ=255)         ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:2372874/68‑1 (MQ=255)         ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:516950/68‑1 (MQ=255)         ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:1456873/68‑1 (MQ=255)         ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:1326090/68‑1 (MQ=255)         ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:1080134/68‑1 (MQ=255)         ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:3351112/68‑1 (MQ=255)         ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:2967308/68‑1 (MQ=255)              aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCTCCTTGCGGCTAACCCCATc                                         >  1:871807/1‑69 (MQ=255)              aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:1743163/1‑69 (MQ=255)              aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:2320720/1‑69 (MQ=255)              aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:876789/1‑69 (MQ=255)              aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:848908/1‑69 (MQ=255)              aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:951339/1‑69 (MQ=255)              aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:1582920/1‑69 (MQ=255)              aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:1242018/1‑69 (MQ=255)              aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:1444361/1‑69 (MQ=255)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAAcccc                                            >  1:2962442/1‑62 (MQ=14)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:317776/1‑69 (MQ=14)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:536257/1‑69 (MQ=14)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:2338499/1‑69 (MQ=14)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:1120658/1‑69 (MQ=14)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:2350209/1‑69 (MQ=14)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:1187568/1‑69 (MQ=14)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:2825189/1‑69 (MQ=14)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:3207478/1‑69 (MQ=14)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:2416487/1‑69 (MQ=14)                  ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCACATCAAtt                                     >  1:514847/1‑69 (MQ=14)                   tataGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTa                                    <  1:2337875/69‑1 (MQ=255)                   tataGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTa                                    <  1:2496967/69‑1 (MQ=255)                   tataGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTa                                    <  1:2297513/69‑1 (MQ=255)                      aGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:1475938/55‑1 (MQ=255)                          aCGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTcc                             >  1:1877348/1‑69 (MQ=255)                          aCGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTcc                             >  1:1934083/1‑69 (MQ=255)                          aCGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTcc                             >  1:1753492/1‑69 (MQ=255)                          aCGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTcc                             >  1:854352/1‑69 (MQ=255)                             gccgccGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGc                          <  1:3256986/69‑1 (MQ=255)                                  cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:433049/69‑1 (MQ=14)                                  cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:2740235/69‑1 (MQ=14)                                  cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:3018823/69‑1 (MQ=14)                                  cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:45192/69‑1 (MQ=14)                                  cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:2935830/69‑1 (MQ=14)                                  cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:207422/69‑1 (MQ=14)                                  cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:1883889/69‑1 (MQ=14)                                  cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:763638/69‑1 (MQ=14)                                  cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:1191215/69‑1 (MQ=14)                                  cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:1074535/69‑1 (MQ=14)                                   gTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:692941/68‑1 (MQ=14)                                   gTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:354386/68‑1 (MQ=14)                                   gTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:1252546/68‑1 (MQ=14)                                     ttACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCa                  <  1:2385606/69‑1 (MQ=255)                                     ttACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCa                  <  1:2107014/69‑1 (MQ=255)                                     ttACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCa                  <  1:83708/69‑1 (MQ=255)                                        ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:1003773/68‑1 (MQ=255)                                        ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:2809240/68‑1 (MQ=255)                                        ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:1604566/68‑1 (MQ=255)                                        ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:1580331/68‑1 (MQ=255)                                        ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:1276402/68‑1 (MQ=255)                                        ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:1016355/68‑1 (MQ=255)                                          ggggCTTCGATCAAGAGCTTCGCGTTGCCGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCg              >  1:1397176/1‑68 (MQ=255)                                          ggggCTTCGATCAAGAGCTTCGCGTTGCCGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGt             >  1:624051/1‑69 (MQ=255)                                          ggggCTTCGATCAAGAGCTTCGCGTTGCCGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGt             >  1:1013263/1‑69 (MQ=255)                                          ggggCTTCGATCAAGAGCTTCGCGTTGCCGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGt             >  1:1639011/1‑69 (MQ=255)                                          ggggCTTCGATCAAGAGCTTCGCGTTGCCGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGt             >  1:2070560/1‑69 (MQ=255)                                          ggggCTTCGATCAAGAGCTTCGCGTTGCCGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGt             >  1:2097343/1‑69 (MQ=255)                                            ggCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTca           <  1:2781477/69‑1 (MQ=255)                                            ggCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTca           <  1:2719039/69‑1 (MQ=255)                                             gCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTa                                    >  1:2506085/1‑43 (MQ=255)                                              cTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:570561/57‑1 (MQ=14)                                                   aTCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCg      <  1:702836/67‑1 (MQ=255)                                                   aTCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACAGGGCAGGCGTCACACCg      <  1:1870270/67‑1 (MQ=255)                                                    tCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:1390575/51‑1 (MQ=14)                                                    tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:2175533/68‑1 (MQ=255)                                                    tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:1661492/68‑1 (MQ=255)                                                    tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:1195084/68‑1 (MQ=255)                                                    tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:1173970/68‑1 (MQ=255)                                                    tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:1038847/68‑1 (MQ=255)                                                    tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:989871/68‑1 (MQ=255)                                                     cAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGtata  >  1:2911836/1‑69 (MQ=255)                                                     cAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGtata  >  1:160379/1‑69 (MQ=255)                                                            ttCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCg      <  1:265810/58‑1 (MQ=255)                                                                 cTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCg      <  1:1858863/53‑1 (MQ=255)                                                                    |                                                      TCGCTACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGTATA  >  minE/2359158‑2359278

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