Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	204,556	A→G	intergenic (+156/‑28)	alaV → / → rrlH	tRNA‑Ala/23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	204,556	0	A	G	100.0%	229.3 / NA	73	intergenic (+156/‑28)	alaV/rrlH	tRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/47);  total (26/47)

GAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTG  >  minE/204497‑204609                                                            |                                                      gAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaa                                             <  1:4358461/70‑1 (MQ=255) gAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaa                                             <  1:499243/70‑1 (MQ=255) gAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaa                                             <  1:4770654/70‑1 (MQ=255) gAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaa                                             <  1:839955/70‑1 (MQ=255) gAAACACTGAACAATGAAAGTTGTTCGTGAGCCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaa                                             <  1:4599594/70‑1 (MQ=11)  aaaGACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                           <  1:4651130/71‑1 (MQ=11)  aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaa                                             <  1:2020081/69‑1 (MQ=255)  aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                           <  1:3252627/71‑1 (MQ=255)  aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                           <  1:2345633/71‑1 (MQ=255)  aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                           <  1:2710962/71‑1 (MQ=255)  aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                           <  1:2852114/71‑1 (MQ=255)  aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                           <  1:1750298/71‑1 (MQ=255)  aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                           <  1:4781840/71‑1 (MQ=255)  aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                           <  1:3895398/71‑1 (MQ=255)  aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                           <  1:604297/71‑1 (MQ=255)       cTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCtt                                      >  1:2919056/1‑71 (MQ=255)        tGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaa                                             <  1:1235814/63‑1 (MQ=255)            caatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATc                                        >  1:793553/1‑64 (MQ=255)            caatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCt                                       >  1:1356349/1‑65 (MQ=255)            caatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGt                                 >  1:3509251/1‑71 (MQ=255)            caatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGt                                 >  1:3885044/1‑71 (MQ=255)            caatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGt                                 >  1:4140465/1‑71 (MQ=255)            caatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGt                                 >  1:4093812/1‑71 (MQ=255)            caatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGt                                 >  1:433303/1‑71 (MQ=255)            caatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGt                                 >  1:2091433/1‑71 (MQ=255)            caatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGt                                 >  1:1418037/1‑71 (MQ=255)                 aaaTTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGa                            <  1:4113060/71‑1 (MQ=255)                 aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGa                            <  1:4610927/71‑1 (MQ=21)                 aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGa                            <  1:2098954/71‑1 (MQ=21)                 aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGa                            <  1:4203595/71‑1 (MQ=21)                          cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:2441704/70‑1 (MQ=25)                          cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:4599483/70‑1 (MQ=25)                          cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:19506/70‑1 (MQ=25)                          cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:3297281/70‑1 (MQ=25)                          cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:4676415/70‑1 (MQ=25)                          cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:1571374/70‑1 (MQ=25)                          cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:746247/70‑1 (MQ=25)                                  tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg           >  1:4616321/1‑71 (MQ=21)                                  tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg           >  1:666625/1‑71 (MQ=21)                                  tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg           >  1:1272702/1‑71 (MQ=21)                                  tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg           >  1:2903728/1‑71 (MQ=21)                                  tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg           >  1:2655170/1‑71 (MQ=21)                                  tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg           >  1:890888/1‑71 (MQ=21)                                  tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg           >  1:2091218/1‑71 (MQ=21)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:1067690/60‑1 (MQ=255)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:1393355/60‑1 (MQ=255)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:4386470/60‑1 (MQ=255)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:4749356/60‑1 (MQ=255)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:2376266/60‑1 (MQ=255)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa         >  1:3954208/1‑71 (MQ=21)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa         >  1:494860/1‑71 (MQ=21)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa         >  1:4917213/1‑71 (MQ=21)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa         >  1:4670822/1‑71 (MQ=21)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa         >  1:227181/1‑71 (MQ=21)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa         >  1:3515980/1‑71 (MQ=21)                                    tcAAATTTTCGCAACACGATGATGGATCGCAAGAAAAATCTTCGGGTTGTGAGGTTAAGc                    <  1:3777296/60‑1 (MQ=11)                                    tcAAATTTTCGCAACACGATGATGGATCACAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:1227792/60‑1 (MQ=11)                                    tcAAAATTTCGCAACACGATGATGGATCGCAAGAAACATc                                        >  1:4572985/1‑40 (MQ=14)                                     cAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                    <  1:1297759/59‑1 (MQ=255)                                      aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGa                  >  1:550415/1‑60 (MQ=255)                                      aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGACGTTAAGCGa                  >  1:1854814/1‑60 (MQ=11)                                       aaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTacac      <  1:2227839/71‑1 (MQ=21)                                       aaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTacac      <  1:3619939/71‑1 (MQ=21)                                            tCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg  <  1:791804/70‑1 (MQ=25)                                            tCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg  <  1:751543/70‑1 (MQ=25)                                            tCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg  <  1:1073871/70‑1 (MQ=25)                                            tCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg  <  1:479393/70‑1 (MQ=25)                                            tCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg  <  1:3536802/70‑1 (MQ=25)                                            tCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg  <  1:2368575/70‑1 (MQ=25)                                            tCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg  <  1:2050846/70‑1 (MQ=25)                                             cGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg  <  1:3398223/69‑1 (MQ=25)                                             cGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg  <  1:2042894/69‑1 (MQ=25)                                             cGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg  <  1:967840/69‑1 (MQ=25)                                                            |                                                      GAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTG  >  minE/204497‑204609

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