Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,415,812	G→T	23.8%	noncoding (67/77 nt)	proM ←	tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,415,812	0	G	T	23.8%	99.2 / 16.8	42	noncoding (67/77 nt)	proM	tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base G (18/14);  new base T (5/5);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GGAGGGTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGA  >  minE/2415746‑2415879                                                                   |                                                                    ggagggTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGagagag                                                                 <  1:694645/71‑1 (MQ=255)      gTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGa                                                               <  1:3141884/68‑1 (MQ=255)                gCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTCGTCGGCGAGAGAGGATTCGAAc                                                        >  1:657077/1‑65 (MQ=255)                 cacaGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAc                                                 >  1:5593864/1‑71 (MQ=255)                 cacaGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAc                                                 >  1:1024987/1‑71 (MQ=255)                  acaGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAcc                                                <  1:4330351/71‑1 (MQ=255)                  acaGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAcc                                                <  1:3691230/71‑1 (MQ=255)                   caGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAccc                                               <  1:4726192/71‑1 (MQ=255)                   caGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAccc                                               <  1:4732647/71‑1 (MQ=255)                   caGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAcc                                                >  1:3209176/1‑70 (MQ=255)                     gAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCAc                                             >  1:2324492/1‑71 (MQ=255)                         aCAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCa                                              >  1:4755699/1‑66 (MQ=255)                         aCAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTgg                                          >  1:548331/1‑70 (MQ=255)                             aaaaCCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTccc                                      >  1:3533273/1‑70 (MQ=255)                             aaaaCCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCa                                     <  1:4770576/71‑1 (MQ=255)                                aCCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGAcc                                                <  1:5360046/57‑1 (MQ=255)                                   ccgccgAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGa                                                          <  1:4475856/44‑1 (MQ=255)                                        gAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                          >  1:3122587/1‑71 (MQ=255)                                        gAAGCGGGGTTAAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                          >  1:3793516/1‑71 (MQ=255)                                           gCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                          >  1:577056/1‑68 (MQ=255)                                              gggTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                 <  1:5097850/58‑1 (MQ=255)                                               ttttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCag                   >  1:3108423/3‑70 (MQ=12)                                                 ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCag                   <  1:2367851/69‑2 (MQ=12)                                                 ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                 >  1:5014009/1‑71 (MQ=12)                                                 ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                 >  1:468949/1‑71 (MQ=12)                                                 ttCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgc                           <  1:4029939/61‑1 (MQ=255)                                                   cAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                       >  1:3907592/1‑63 (MQ=255)                                                   cAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTg               >  1:4132185/1‑71 (MQ=255)                                                    aaaaTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCa                                              >  1:1731257/1‑39 (MQ=255)                                                    aaaaTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCa                                <  1:1771869/53‑1 (MQ=255)                                                       tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc           <  1:1308051/70‑2 (MQ=16)                                                       tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc           <  1:4180925/70‑2 (MQ=16)                                                       tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg             >  1:2772012/2‑69 (MQ=16)                                                       tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg             >  1:1922640/2‑69 (MQ=16)                                                       aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                       >  1:4671372/1‑59 (MQ=255)                                                       aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                       >  1:1358257/1‑59 (MQ=255)                                                       aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                       >  1:2281846/1‑59 (MQ=255)                                                        ttGGTCGGTGATAGAGGATTCGAACCTCCGAcccc                                              <  1:1540883/35‑2 (MQ=18)                                                        ttGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                 <  1:1723109/48‑1 (MQ=255)                                                        ttGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAAcc                                 <  1:4283347/48‑1 (MQ=255)                                                          ggTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTAc                      <  1:3823222/57‑1 (MQ=12)                                                              ggCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGcc    <  1:1147009/71‑1 (MQ=255)                                                               gCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGcc    >  1:1599202/1‑70 (MQ=255)                                                                 gagagaGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGa  >  1:3291001/1‑70 (MQ=255)                                                                   |                                                                    GGAGGGTGACGAAATGCACAGAAAACAAAAAACCCCGCCGAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGA  >  minE/2415746‑2415879

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