Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,886,204	G→A	intergenic (‑34/+635)	yobF ← / ← yebO	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,886,204	0	G	A	100.0%	114.9 / NA	31	intergenic (‑34/+635)	yobF/yebO	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/13);  total (18/13)

GTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATAC  >  REL606/1886132‑1886276                                                                         |                                                                         gtgtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTa                                                                       <  1:201846/76‑1 (MQ=255)  tgtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCtattat                                                                      <  1:194114/76‑1 (MQ=255)  tgtTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCtattat                                                                      <  2:315276/76‑1 (MQ=255)    tttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATca                                                                    <  2:860031/76‑1 (MQ=255)      tACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATcaca                                                                  >  1:1045170/1‑76 (MQ=255)      tACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATcaca                                                                  >  1:304766/1‑76 (MQ=255)        cTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATa                                                                <  2:71049/76‑1 (MQ=255)          cAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACa                                                              >  1:882990/1‑76 (MQ=255)           aGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACaa                                                             <  1:796497/76‑1 (MQ=255)           aGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACaa                                                             >  1:828802/1‑76 (MQ=255)           aGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACaa                                                             >  1:828838/1‑76 (MQ=255)                ttcttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTa                                                        <  2:1061868/76‑1 (MQ=255)                     tACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTAt                                                   >  1:80777/1‑76 (MQ=255)                      aCTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATg                                                  >  1:1010960/1‑76 (MQ=255)                          aaaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCaa                                              >  1:182639/1‑76 (MQ=255)                            aaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGc                                            >  2:411338/1‑76 (MQ=255)                            aaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGc                                            >  1:416980/1‑76 (MQ=255)                                 cacacaTAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCAttt                                       >  2:617266/1‑76 (MQ=255)                                     caTAAACAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTaa                                   <  2:258160/76‑1 (MQ=255)                                        aaaCAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACa                                >  1:213571/1‑76 (MQ=255)                                        aaaCAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACa                                >  1:213569/1‑76 (MQ=255)                                        aaaCAGAACTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACa                                >  1:213556/1‑76 (MQ=255)                                              aaCTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGtctc                          <  1:830718/76‑1 (MQ=255)                                               aCTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCg                         <  2:745685/76‑1 (MQ=255)                                               aCTGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCg                         <  2:838539/76‑1 (MQ=255)                                                 tGTACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGTCAAGCCATTTTTAAACATGTCTCGAt                       <  2:268076/76‑1 (MQ=255)                                                   tACCTCGTTTAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCa                     >  2:806801/1‑76 (MQ=255)                                                           ttAACCCGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGcaca             >  2:128584/1‑76 (MQ=255)                                                                 cGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAtctc       <  1:864892/76‑1 (MQ=255)                                                                 cGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAtctc       <  1:414401/76‑1 (MQ=255)                                                                 cGAAATCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAtctc       <  1:414389/76‑1 (MQ=255)                                                                     aTCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATa   >  2:246656/1‑76 (MQ=255)                                                                      tCTATTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATAc  >  2:883548/1‑76 (MQ=255)                                                                         |                                                                         GTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATAC  >  REL606/1886132‑1886276

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