Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,359,227	G→C	35.9%	noncoding (1868/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,227	0	G	C	35.9%	104.4 / 74.0	89	noncoding (1868/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (22/35);  new base C (17/15);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CGCTACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGTATA  >  minE/2359159‑2359278                                                                     |                                                    cGCTACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCg                                                     >  1:2360182/1‑69 (MQ=40) cGCTACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCg                                                     >  1:3429729/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:371375/1‑69 (MQ=255)    tACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCt                                                  >  1:2790557/1‑69 (MQ=255)    tACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCt                                                  >  1:3156533/1‑69 (MQ=255)    tACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCt                                                  <  1:3074728/69‑1 (MQ=255)        ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:516950/68‑1 (MQ=255)        ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:2372874/68‑1 (MQ=255)        ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:1620996/68‑1 (MQ=255)        ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:1456873/68‑1 (MQ=255)        ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:2967308/68‑1 (MQ=255)        ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:1326090/68‑1 (MQ=255)        ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:3351112/68‑1 (MQ=255)        ttAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:1080134/68‑1 (MQ=255)             aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCTCCTTGCGGCTAACCCCATc                                         >  1:871807/1‑69 (MQ=255)             aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:1444361/1‑69 (MQ=255)             aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:876789/1‑69 (MQ=255)             aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:1582920/1‑69 (MQ=255)             aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:848908/1‑69 (MQ=255)             aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:1242018/1‑69 (MQ=255)             aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:951339/1‑69 (MQ=255)             aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:1743163/1‑69 (MQ=255)             aCCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATc                                         >  1:2320720/1‑69 (MQ=255)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAAcccc                                            >  1:2962442/1‑62 (MQ=14)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:1120658/1‑69 (MQ=14)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:536257/1‑69 (MQ=14)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:2338499/1‑69 (MQ=14)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:2350209/1‑69 (MQ=14)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:2416487/1‑69 (MQ=14)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:3207478/1‑69 (MQ=14)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:2825189/1‑69 (MQ=14)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:1187568/1‑69 (MQ=14)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAAtt                                     >  1:317776/1‑69 (MQ=14)                 ttATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCACATCAAtt                                     >  1:514847/1‑69 (MQ=14)                  tataGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTa                                    <  1:2297513/69‑1 (MQ=255)                  tataGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTa                                    <  1:2337875/69‑1 (MQ=255)                  tataGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTa                                    <  1:2496967/69‑1 (MQ=255)                     aGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAAc                                               <  1:1475938/55‑1 (MQ=255)                         aCGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTcc                             >  1:1934083/1‑69 (MQ=255)                         aCGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTcc                             >  1:854352/1‑69 (MQ=255)                         aCGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTcc                             >  1:1877348/1‑69 (MQ=255)                         aCGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTcc                             >  1:1753492/1‑69 (MQ=255)                            gccgccGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGc                          <  1:3256986/69‑1 (MQ=255)                                 cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:45192/69‑1 (MQ=14)                                 cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:2935830/69‑1 (MQ=14)                                 cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:3018823/69‑1 (MQ=14)                                 cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:2740235/69‑1 (MQ=14)                                 cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:763638/69‑1 (MQ=14)                                 cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:433049/69‑1 (MQ=14)                                 cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:1191215/69‑1 (MQ=14)                                 cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:1074535/69‑1 (MQ=14)                                 cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:1883889/69‑1 (MQ=14)                                 cgTTTACCGGGGCTTCGATCAAGAGCTTCGCGTTACCGCTAACCCCATCAATTAACCTTCCGGCACCgg                     <  1:207422/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:83708/69‑1 (MQ=255)                                    ttACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCa                  <  1:2107014/69‑1 (MQ=255)                                    ttACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCa                  <  1:2385606/69‑1 (MQ=255)                                       ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:1580331/68‑1 (MQ=255)                                       ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:1604566/68‑1 (MQ=255)                                       ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:2809240/68‑1 (MQ=255)                                       ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:1276402/68‑1 (MQ=255)                                       ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:1016355/68‑1 (MQ=255)                                       ccGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAgg                <  1:1003773/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:2097343/1‑69 (MQ=255)                                         ggggCTTCGATCAAGAGCTTCGCGTTGCCGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGt             >  1:2070560/1‑69 (MQ=255)                                         ggggCTTCGATCAAGAGCTTCGCGTTGCCGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGt             >  1:1639011/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:989871/68‑1 (MQ=255)                                                   tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:1038847/68‑1 (MQ=255)                                                   tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:1173970/68‑1 (MQ=255)                                                   tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:1195084/68‑1 (MQ=255)                                                   tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:1661492/68‑1 (MQ=255)                                                   tCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGta    <  1:2175533/68‑1 (MQ=255)                                                    cAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGtata  >  1:160379/1‑69 (MQ=255)                                                    cAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGtata  >  1:2911836/1‑69 (MQ=255)                                                           ttCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCg      <  1:265810/58‑1 (MQ=255)                                                                cTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCg      <  1:1858863/53‑1 (MQ=255)                                                                     |                                                    CGCTACCTTAGGACCGTTATAGTTACGGCCGCCGTTTACCGGGGCTTCGATCAAGAGCTTCGCCTTGCGGCTAACCCCATCAATTAACCTTCCGGCACCGGGCAGGCGTCACACCGTATA  >  minE/2359159‑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: ‑