Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	688,340	Δ1 bp	7.8%	noncoding (72/88 nt)	serT ←	tRNA‑Ser

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	688,340	0	A	.	7.8%	141.7 / 8.1	51	noncoding (72/88 nt)	serT	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base . (2/2);  total (27/24)
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.49e-01

TCGACAGGCAAAAGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACA  >  minE/688272‑688406                                                                     |                                                                    tCGACAGGCAAAAGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGAt                                                                    <  1:460712/70‑1 (MQ=255)  cGACAGGCAAAAGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTc                                                                  >  1:2454627/1‑71 (MQ=255)   gACAGGCAAAAGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCg                                                                 <  1:2591961/71‑1 (MQ=255)     cAGGCAAAAGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGaa                                                               >  1:531891/1‑71 (MQ=255)     cAGGCAAAAGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGa                                                                >  1:1313780/1‑70 (MQ=255)      aGGCAAAAGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑c                                                             >  1:2454245/1‑71 (MQ=255)           aaaGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTgg                                                        <  1:484431/71‑1 (MQ=255)           aaaGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTg                                                         >  1:946803/1‑70 (MQ=255)           aaaGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTg                                                         >  1:2150932/1‑70 (MQ=255)             aGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTgg                                                        >  1:1102641/1‑69 (MQ=255)                    ggggTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTc                                               >  1:1936522/1‑71 (MQ=255)                      ggTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAG‑TTCGAATCTCTgc                                                        <  1:3539728/60‑2 (MQ=38)                      ggTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAG‑TTCGAATCTCTgc                                                        <  1:1994970/60‑2 (MQ=38)                        ttAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGAt                                                                    <  1:2336811/47‑1 (MQ=255)                        ttAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGAt                                                                    <  1:3312777/47‑1 (MQ=255)                        ttAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGt                                           >  1:2961282/1‑71 (MQ=255)                         tAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTc                                          >  1:27671/1‑71 (MQ=255)                          aGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCg                                         >  1:2652150/1‑71 (MQ=255)                             aTTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCttt                                                <  1:1490700/61‑1 (MQ=255)                              tttAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCgg                                     <  1:2204134/71‑1 (MQ=255)                              tttAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCgg                                     <  1:954192/71‑1 (MQ=255)                               ttAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCggg                                            >  1:2514759/1‑63 (MQ=255)                                tAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGtt                                   >  1:434044/1‑71 (MQ=255)                                tAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGt                                    <  1:743396/70‑1 (MQ=255)                                 aGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAAccc                                                   <  1:1452983/54‑1 (MQ=255)                                  gCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGtttt                                 <  1:1472074/71‑1 (MQ=255)                                   cTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCtt                                                 >  1:28239/1‑54 (MQ=255)                                   cTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTc                                <  1:2978651/71‑1 (MQ=255)                                    tAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGtt                                   <  1:3006430/67‑1 (MQ=255)                                         ccTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTc                                          <  1:2911987/55‑1 (MQ=255)                                         ccTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTc                                          <  1:1608549/55‑1 (MQ=255)                                          cttatcttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCg                         >  1:1973553/1‑71 (MQ=255)                                          cttatcttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCg                         >  1:2851936/1‑71 (MQ=255)                                          cttatcttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCg                         <  1:479776/71‑1 (MQ=255)                                          cttatcttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCg                         >  1:427594/1‑71 (MQ=255)                                           ttatcttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCg                         >  1:283473/1‑70 (MQ=255)                                           ttatcttatTTGGCGGAAGCGCAGAG‑TTCGAATCTCTgc                                                        >  1:746243/1‑38 (MQ=18)                                           ttatcttatTTGGCGGAAGCGCAGAG‑TTCGAATCTCTgc                                                        >  1:922208/1‑38 (MQ=18)                                             atcttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTc                                          <  1:1722859/51‑1 (MQ=255)                                             atcttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCaa                              >  1:2431333/1‑63 (MQ=255)                                             atcttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCaa                              >  1:1112578/1‑63 (MQ=255)                                             atcttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCa                               >  1:1862950/1‑62 (MQ=255)                                               cttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCttt                                                <  1:2881375/43‑1 (MQ=255)                                               cttatTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCttt                                                <  1:1001678/43‑1 (MQ=255)                                                   tttGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCGGTGCCTTc                 >  1:901062/1‑70 (MQ=255)                                                   tttGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCGGTGCCTTCa                <  1:487178/71‑1 (MQ=255)                                                         ggAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCGGTGCCTTCAACCGCt          <  1:3127910/71‑1 (MQ=255)                                                                cAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCcac   <  1:247773/71‑1 (MQ=255)                                                                 agagATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGttt                                  >  1:386749/1‑39 (MQ=255)                                                                 agagATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCcaca  >  1:1963955/1‑71 (MQ=255)                                                                 agagATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCcac   >  1:1700884/1‑70 (MQ=255)                                                                     |                                                                    TCGACAGGCAAAAGAAAAAGGGGTTAGCATTTAGCTAACCCCTTATCTTATTTGGCGGAAGCGCAGAGATTCGAA‑CTCTGGAACCCTTTCGGGTCGCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACA  >  minE/688272‑688406

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