Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,604,103	A→G	intergenic (+171/‑121)	unknown → / → unknown	UPF0053 inner membrane protein YgdQ/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,604,103	0	A	G	100.0%	64.9 / NA	24	intergenic (+171/‑121)	unknown/unknown	UPF0053 inner membrane protein YgdQ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/9);  total (15/9)

GCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTG  >  1/1603861‑1604214                                                                                                                                                                                                                                                   |                                                                                                                gCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTtcat                                                                                                         >  1:197995/1‑251 (MQ=255) gCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTtcat                                                                                                         >  1:470504/1‑251 (MQ=255)         cTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGAttt                                                                                                  >  1:75044/1‑250 (MQ=255)                              tctcAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTg                                                                            >  1:247082/1‑251 (MQ=255)                               ctcAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTgg                                                                           >  1:157598/1‑251 (MQ=255)                                   aTACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGaaaa                                                                       >  1:248147/1‑251 (MQ=255)                                   aTACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGCGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGaaaa                                                                       >  1:429464/1‑251 (MQ=255)                                     aCCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAg                                                                      >  1:518045/1‑250 (MQ=255)                                              cGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGAGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGa                                                            >  1:471639/1‑251 (MQ=255)                                               gTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGaa                                                           >  1:286532/1‑251 (MQ=255)                                                   aCGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAagag                                                       <  1:121836/251‑1 (MQ=255)                                                    cGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAag                                                         >  1:212050/1‑248 (MQ=255)                                                              cTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg                                             <  1:292611/250‑1 (MQ=255)                                                                           aaaCAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCtt                               >  1:566179/1‑251 (MQ=255)                                                                           aaaCAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCtt                               >  1:249735/1‑251 (MQ=255)                                                                           aaaCAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCt                                <  1:465716/250‑1 (MQ=255)                                                                             aCAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTtctc                                                                                               <  1:764170/185‑1 (MQ=17)                                                                              cAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtata                            <  1:535620/251‑1 (MQ=255)                                                                               aaaaGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                           <  1:657698/251‑1 (MQ=255)                                                                                      cGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAccc                     <  1:342099/250‑1 (MQ=255)                                                                                                 cacaTTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTaa         >  1:150356/1‑251 (MQ=255)                                                                                                   caTTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTaaaa       >  1:46503/1‑251 (MQ=255)                                                                                                        aaTTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg  <  1:688941/251‑1 (MQ=255)                                                                                                        aaTTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg  <  1:336185/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                GCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTG  >  1/1603861‑1604214

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