Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,235,443	A→G	intergenic (+147/‑145)	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,235,443	0	A	G	100.0%	76.0 / NA	29	intergenic (+147/‑145)	unknown/unknown	UPF0053 inner membrane protein YgdQ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/12);  total (17/12)

TCCATGTGCCGAAAGCGTATATTTACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGAATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTG  >  1/1235196‑1235578                                                                                                                                                                                                                                                        |                                                                                                                                        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)                         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=18)                                                                  aCCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAg                                                                      >  1:518045/1‑250 (MQ=255)                                                                           cGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGAGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGa                                                            >  1:471639/1‑251 (MQ=18)                                                                            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:249735/1‑251 (MQ=255)                                                                                                        aaaCAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCtt                               >  1:566179/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)                                                                                                                                                                                                                                                        |                                                                                                                                        TCCATGTGCCGAAAGCGTATATTTACTCCGCTATGGGCTTTTCCGTATTTGTGGAAATTCTCAATACCCAAATGCGTAAACGCTTGGGCGCTAAAATTTAGAGAAACAAAAGTGCGGTCGAAAAACACATTAAATTTCAACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGAATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTG  >  1/1235196‑1235578

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