Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,279,749	Δ1 bp	coding (633/1101 nt)	trmA →	tRNA (uracil‑5‑)‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,279,746	0	A	.	100.0%	168.4 / NA	41	coding (630/1101 nt)	trmA	tRNA (uracil‑5‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (22/19);  total (22/19)

CTTTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAATT  >  minE/2279677‑2279810                                                                      |                                                                 cTTTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaa                                                                 <  1:138024/71‑1 (MQ=255)  tttACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCa                                                                  >  1:902724/1‑69 (MQ=255)   ttACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaag                                                               >  1:2065599/1‑70 (MQ=255)      cccAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaa                                                                 <  1:1137642/66‑1 (MQ=255)       ccAGCCGAACTCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGa                                                          <  1:49140/71‑1 (MQ=255)       ccAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGa                                                          <  1:486080/71‑1 (MQ=255)          gCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGAGGTAACCAAAGGCTCaaa                                                                >  1:2079147/1‑63 (MQ=255)          gCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaa                                                                >  1:773499/1‑63 (MQ=255)          gCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaa                                                                >  1:2870589/1‑63 (MQ=255)          gCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGAtt                                                        >  1:526759/1‑70 (MQ=255)             gAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaa                                                                 >  1:627544/1‑59 (MQ=255)             gAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaa                                                                 >  1:623925/1‑59 (MQ=255)             gAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACt                                                    <  1:504121/71‑1 (MQ=255)              aaCGCGGCGATGAATATTCAGATGCTGTAATGGGCGCTGGACGTAACCAAAGGCTCaaagg                                                              <  1:1305348/61‑3 (MQ=255)              aaCGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaagg                                                              >  1:326786/1‑59 (MQ=255)              aaCGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTg                                                   >  1:1908356/1‑71 (MQ=255)                     cgATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTg                                             >  1:2782639/1‑70 (MQ=255)                       aTGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTa                                           >  1:2434595/1‑70 (MQ=255)                        tGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACt                                         >  1:2292245/1‑71 (MQ=255)                            tatTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCgg                                     >  1:1832283/1‑71 (MQ=255)                            tatTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCgg                                     <  1:605004/71‑1 (MQ=255)                             atTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCgg                                     >  1:2413812/1‑70 (MQ=255)                              ttCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGc                                    <  1:486195/70‑1 (MQ=255)                              ttCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGc                                    >  1:1614837/1‑70 (MQ=255)                               tCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGa                                                 >  1:1085051/1‑56 (MQ=255)                               tCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCaa                                  >  1:1426746/1‑71 (MQ=255)                                 aGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACg                                >  1:2810234/1‑71 (MQ=255)                                   aTGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGt                              <  1:1823034/71‑1 (MQ=255)                                    tGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaa                                                                >  1:1285997/1‑37 (MQ=255)                                    tGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaa                                                                >  1:1524275/1‑37 (MQ=255)                                    tGCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGt                              >  1:2846077/1‑70 (MQ=255)                                     gCTGGAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTAc                                                     <  1:663136/46‑1 (MQ=39)                                         gAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACttt                        <  1:2091749/71‑1 (MQ=255)                                         gAATGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACt                          <  1:2824682/69‑1 (MQ=255)                                           aTGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                       >  1:2288780/1‑70 (MQ=255)                                            tGGGCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCaa                                  <  1:2604592/58‑1 (MQ=255)                                             gggCGCTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCAt                    >  1:2696620/1‑71 (MQ=255)                                               gcgcTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTa                  <  1:805630/71‑1 (MQ=255)                                                cgcTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAg                 <  1:672176/71‑1 (MQ=255)                                                  cTGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAgcg               <  1:368073/71‑1 (MQ=255)                                                   tGGACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtt                         <  1:2137988/60‑1 (MQ=255)                                                     gACGTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGc                                    <  1:2319796/47‑1 (MQ=38)                                                        gTAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTg                                             <  1:1862145/35‑1 (MQ=38)                                                         tAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCaa                                  <  1:128745/45‑1 (MQ=38)                                                         tAACCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTa                             <  1:509650/50‑1 (MQ=39)                                                          aaCCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGcgcg       >  1:2839954/1‑71 (MQ=255)                                                          aaCCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGcgcg       >  1:1551248/1‑71 (MQ=255)                                                           aCCAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAAc                           >  1:2528020/1‑50 (MQ=39)                                                            ccAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt      >  1:361146/1‑70 (MQ=255)                                                            ccAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt      >  1:1313903/1‑70 (MQ=255)                                                             cAAAGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa    <  1:2849682/71‑1 (MQ=255)                                                               aaGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAAtt  >  1:1558746/1‑71 (MQ=255)                                                               aaGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAAtt  >  1:899044/1‑71 (MQ=255)                                                               aaGGCTC‑AAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAAtt  >  1:902787/1‑71 (MQ=255)                                                                      |                                                                 CTTTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAATT  >  minE/2279677‑2279810

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