Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,604,050	A→T	intergenic (+118/‑174)	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,050	0	A	T	100.0%	91.4 / NA	33	intergenic (+118/‑174)	unknown/unknown	UPF0053 inner membrane protein YgdQ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (17/16);  total (17/16)

TTATCGAGAGCTTCAGCATCCATGTGCCGAAAGCGTATATTTACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTG  >  1/1603814‑1604214                                                                                                                                                                                                                                             |                                                                                                                                                                     ttATCGAGAGCTTCAGCATCCATGTGCCGAAAGCGTATATTTACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAAc                                                                                                                                                        <  1:14030/251‑1 (MQ=255)    tCGAGAGCTTCAGCATCCATGTGCCGAAAGCGTATATTTACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAgg                                                                                                                                                     <  1:54277/251‑1 (MQ=255)                   tCCATGTGCCGAAAGCGTATATTTACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGAtt                                                                                                                                      >  1:585994/1‑251 (MQ=255)                       tgtgCCGAAAGCGTATATTTACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTgg                                                                                                                                  <  1:895989/251‑1 (MQ=255)                                  cGTATATTTACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCG‑AAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGt                                                                                                                      <  1:245913/251‑1 (MQ=255)                                     atatTTACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTgg                                                                                                                    >  1:519324/1‑251 (MQ=255)                                         ttACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAAc                                                                                                                                                        <  1:581617/211‑1 (MQ=255)                                           aCTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCtgaac                                                                                                              <  1:51218/251‑5 (MQ=255)                                                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)                                                                                                                                       gTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAAc                                                                                                                                                        <  1:810272/117‑1 (MQ=17)                                                                                                                                                cacaTTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTaa         >  1:150356/1‑251 (MQ=255)                                                                                                                                                  caTTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTaaaa       >  1:46503/1‑251 (MQ=255)                                                                                                                                                       aaTTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg  <  1:336185/251‑1 (MQ=255)                                                                                                                                                       aaTTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg  <  1:688941/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                     TTATCGAGAGCTTCAGCATCCATGTGCCGAAAGCGTATATTTACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTG  >  1/1603814‑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: ‑