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%	176.5 / NA	56	intergenic (+156/‑28)	alaV/rrlH	tRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/19);  total (37/19)

AAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCA  >  minE/204492‑204621                                                                 |                                                                  aaaTTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGc                                                              >  1:1186791/1‑70 (MQ=255) aaaTTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGc                                                              >  1:2469101/1‑70 (MQ=255)   aTTGAAACACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGc                                                              >  1:949870/1‑68 (MQ=255)       aaaCACTGAACAATGAAAGTTGTTGGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                                       <  1:1945073/71‑1 (MQ=11)       aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaaa                                                        >  1:2142937/1‑70 (MQ=255)       aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGaaa                                                        >  1:25514/1‑70 (MQ=255)       aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                                       <  1:2085822/71‑1 (MQ=255)       aaaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                                       <  1:1163429/71‑1 (MQ=255)        aaCACTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACa                                                      >  1:1606746/1‑71 (MQ=255)         acacTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACa                                                      >  1:831749/1‑70 (MQ=255)         acacTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACAt                                                     >  1:1801252/1‑71 (MQ=255)         acacTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACAt                                                     >  1:1711261/1‑71 (MQ=255)          cacTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATc                                                    >  1:2329849/1‑71 (MQ=255)          cacTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATc                                                    >  1:2866282/1‑71 (MQ=255)           acTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCt                                                   <  1:768425/71‑1 (MQ=255)               aacaatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCgg                                               >  1:550037/1‑71 (MQ=255)                acaatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCggg                                              >  1:1225130/1‑71 (MQ=255)                acaatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCggg                                              >  1:812463/1‑71 (MQ=255)                 caatGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGt                                             >  1:2741695/1‑71 (MQ=255)                  aatgAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt                                            <  1:905973/68‑1 (MQ=255)                      aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgtg                                         <  1:832148/70‑1 (MQ=25)                            gttCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtg                                           <  1:1256520/62‑1 (MQ=25)                            gttCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGtt                                    >  1:347715/1‑69 (MQ=25)                               cGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc                                <  1:986707/70‑1 (MQ=25)                                gTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGa                              <  1:2263840/71‑1 (MQ=21)                                   aGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTa                                   >  1:969024/1‑63 (MQ=255)                                    gTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGtt                                    >  1:136573/1‑61 (MQ=255)                                    gTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGtt                                    >  1:525151/1‑61 (MQ=255)                                    gTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTaa                          <  1:1643008/71‑1 (MQ=21)                                       tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTc                                                 >  1:1755983/1‑45 (MQ=255)                                       tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTc                                                 >  1:2950967/1‑45 (MQ=255)                                       tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                       >  1:102014/1‑71 (MQ=21)                                        ctcAAATTTTCGCAACACGATGATGGATCGCAAGAAAc                                                       <  1:1788831/38‑1 (MQ=255)                                         tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATc                                                    <  1:2333897/40‑1 (MQ=255)                                         tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCg                               <  1:2306813/61‑1 (MQ=25)                                         tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                       <  1:2910091/69‑1 (MQ=25)                                          cAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtg                                           >  1:113660/1‑48 (MQ=255)                                          cAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa                     >  1:682149/1‑70 (MQ=25)                                           aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGt                      >  1:1355468/1‑68 (MQ=25)                                           aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGt                      <  1:2701579/68‑1 (MQ=25)                                           aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTaca                   >  1:713227/1‑71 (MQ=21)                                           aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTaca                   >  1:841607/1‑71 (MQ=21)                                           aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa                     <  1:1308449/69‑1 (MQ=25)                                             aTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                       >  1:1844736/1‑65 (MQ=25)                                             aTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                       >  1:642065/1‑65 (MQ=25)                                             aTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                       >  1:643689/1‑65 (MQ=25)                                             aTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTacac                  <  1:2925776/70‑1 (MQ=25)                                             aTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACg                 >  1:753676/1‑71 (MQ=21)                                                 tCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg              <  1:470335/70‑1 (MQ=25)                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc       >  1:2452985/1‑71 (MQ=21)                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc       >  1:279436/1‑71 (MQ=21)                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGccc       >  1:1060693/1‑71 (MQ=21)                                                       cacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGcc        >  1:972401/1‑70 (MQ=25)                                                          gatgatGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTa                                   <  1:570898/40‑1 (MQ=255)                                                          gatgatGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTgg    >  1:1693856/1‑71 (MQ=21)                                                          gatgatGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTgg    >  1:1517525/1‑71 (MQ=21)                                                            tgatggATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCa  <  1:1271710/71‑1 (MQ=21)                                                            tgatggATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCa  <  1:1268059/71‑1 (MQ=21)                                                              atggATCGCAAGAAACATCTTCGGGTTGTGAGGtt                                    <  1:2218773/35‑1 (MQ=255)                                                                 |                                                                  AAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCA  >  minE/204492‑204621

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