Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	821,366	A→G	T4A (ACC→GCC)	unknown →	DNA‑binding domain of ModE / Molybdate‑binding domain of ModE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	821,366	0	A	G	100.0%	223.6 / NA	78	T4A (ACC→GCC)	unknown	DNA‑binding domain of ModE / Molybdate‑binding domain of ModE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (52/26);  total (52/26)

AAATTATATATCGCGCAGATTATACAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTACCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGT  >  1/821160‑821598                                                                                                                                                                                                               |                                                                                                                                                                                                                                         aaatTATATATCGCGCAGATTATACAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGaaa                                                                                                                                                                                                                                    >  1:406896/1‑213 (MQ=255)   atTATATATCGCGCAGATTATACAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTc                                                                                                                                                                                                                            <  1:803207/219‑1 (MQ=255)   atTATATATCGCGCAGATTATACAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAActc                                                                                                                                                                                                    <  1:98024/243‑1 (MQ=255)   atTATATATCGCGCAGATTATACAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGa                                                                                                                                                                                            >  1:75760/1‑251 (MQ=255)                 aGATTATACAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATc                                                                                                                                                                               >  1:530382/1‑250 (MQ=255)                  gATTATACAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAActc                                                                                                                                                                                                    >  1:98286/1‑228 (MQ=255)                    ttATACAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGtt                                                                                                                                                                           >  1:143890/1‑251 (MQ=255)                      ataCAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTc                                                                                                                                                                                                                            >  1:223610/1‑200 (MQ=255)                      ataCAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTtta                                                                                                                                                                         >  1:901144/1‑251 (MQ=255)                        aCAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTttatt                                                                                                                                                                       >  1:686345/1‑251 (MQ=255)                                  gcgcAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAAtt                                                                                                                                                             >  1:517917/1‑251 (MQ=255)                                     cAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGaa                                                                                                                                                          >  1:78648/1‑251 (MQ=255)                                                  tttttGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCa                                                                                                                                             >  1:772713/1‑251 (MQ=255)                                                        aaCAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAAcc                                                                                                                                       <  1:803657/251‑1 (MQ=255)                                                          cAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCaa                                                                                                                                     >  1:117293/1‑251 (MQ=255)                                                          cAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCaa                                                                                                                                     >  1:710121/1‑251 (MQ=255)                                                            aGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTAATTAAAAGAAATTGAAAAATGCGGCTCCAtt                                                                                                                                           <  1:229286/243‑1 (MQ=255)                                                                  gATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCaaa                                                                                                                             >  1:300327/1‑251 (MQ=255)                                                                  gATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCaaa                                                                                                                             >  1:616782/1‑251 (MQ=255)                                                                    tCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCaaaaa                                                                                                                           >  1:595979/1‑251 (MQ=255)                                                                    tCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCaaaaa                                                                                                                           >  1:124664/1‑251 (MQ=255)                                                                       caAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATc                                                                                                                                                                               >  1:740951/1‑196 (MQ=255)                                                                          gATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCa                                                                                                                     >  1:324397/1‑251 (MQ=255)                                                                                 tATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAg                                                                                                              <  1:522123/251‑1 (MQ=255)                                                                                  aTGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGt                                                                                                             <  1:495053/251‑1 (MQ=255)                                                                                          aaCGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAg                                                                                                     <  1:633813/251‑1 (MQ=255)                                                                                              aTTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAAcc                                                                                                                                       >  1:808897/1‑213 (MQ=255)                                                                                                         ttCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCTACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAg                                                                                                     <  1:449669/236‑1 (MQ=255)                                                                                                         ttCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAg                                                                                                     <  1:79439/236‑1 (MQ=255)                                                                                                           cTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCa                                                                                                                                                                                              >  1:398841/1‑145 (MQ=255)                                                                                                           cTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTAc                                                                                                          >  1:64018/1‑229 (MQ=255)                                                                                                              aTTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAg                                                                                 >  1:389199/1‑251 (MQ=255)                                                                                                                     aCAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAg                                                                                                     <  1:675990/224‑1 (MQ=255)                                                                                                                     aCAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGaa                                                                          >  1:675523/1‑251 (MQ=255)                                                                                                                          atatAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAAcc                                                                                                                                       >  1:467293/1‑185 (MQ=255)                                                                                                                          atatAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATc                                                                                          >  1:403316/1‑230 (MQ=255)                                                                                                                          atatAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATaaaa                                                                     >  1:569708/1‑251 (MQ=255)                                                                                                                          atatAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATaaaa                                                                     >  1:534900/1‑251 (MQ=255)                                                                                                                          atatAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATaaaa                                                                     >  1:419142/1‑251 (MQ=255)                                                                                                                           tataGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAAt                                                                    <  1:354361/251‑1 (MQ=255)                                                                                                                            ataGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAAcc                                                                                                                                       <  1:191325/183‑1 (MQ=255)                                                                                                                            ataGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAAc                                                                                                                                        <  1:198504/182‑1 (MQ=255)                                                                                                                            ataGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGa                                                                                            <  1:512122/226‑1 (MQ=255)                                                                                                                            ataGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAAt                                                                    >  1:122491/1‑250 (MQ=255)                                                                                                                            ataGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAAt                                                                    >  1:720434/1‑250 (MQ=255)                                                                                                                            ataGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAAt                                                                    >  1:468673/1‑250 (MQ=255)                                                                                                                            ataGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATc                                                                   <  1:527342/251‑1 (MQ=255)                                                                                                                            ataGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATc                                                                   <  1:365516/251‑1 (MQ=255)                                                                                                                               gTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCAAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAg                                                                 >  1:807315/1‑250 (MQ=255)                                                                                                                                 atgatgTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGcc                                                               >  1:874449/1‑250 (MQ=255)                                                                                                                                    atgTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATc                                                                                                                                                                                          >  1:618758/1‑124 (MQ=255)                                                                                                                                    atgTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGc                                                                                <  1:279872/230‑1 (MQ=255)                                                                                                                                    atgTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGaa                                                           >  1:393154/1‑251 (MQ=255)                                                                                                                                    atgTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGaa                                                           >  1:646735/1‑251 (MQ=255)                                                                                                                                        gCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCg                                                       >  1:819983/1‑251 (MQ=255)                                                                                                                                            ccGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGa                                                                                                                                                                 >  1:156411/1‑141 (MQ=255)                                                                                                                                            ccGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAg                                                                                                                  <  1:745731/188‑1 (MQ=255)                                                                                                                                            ccGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGcc                                                               <  1:842204/239‑1 (MQ=255)                                                                                                                                            ccGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTa                                                    >  1:659777/1‑250 (MQ=255)                                                                                                                                             cGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCa                                                                                                                                      <  1:851454/167‑1 (MQ=255)                                                                                                                                                        gTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATc                                                                   >  1:160796/1‑223 (MQ=255)                                                                                                                                                            tAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATcgg                                   <  1:709985/251‑1 (MQ=255)                                                                                                                                                              aaaCGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGc                                                                                <  1:293798/204‑1 (MQ=255)                                                                                                                                                                aCGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGaaa                                                                                                                                                         >  1:796430/1‑129 (MQ=255)                                                                                                                                                                aCGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATcggcgg                                >  1:338136/1‑250 (MQ=255)                                                                                                                                                                  gATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATc                                                                   <  1:285410/213‑1 (MQ=255)                                                                                                                                                                   atatatTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGaaa                                                                                                                                                         >  1:635633/1‑126 (MQ=255)                                                                                                                                                                   atatatTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAAt                                                                    <  1:866659/211‑1 (MQ=255)                                                                                                                                                                     atatTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAAcc                                                                                                                                       >  1:58919/1‑142 (MQ=255)                                                                                                                                                                     atatTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAAcc                                                                                                                                       >  1:19139/1‑142 (MQ=255)                                                                                                                                                                     atatTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAAt                                                                    <  1:497792/209‑1 (MQ=255)                                                                                                                                                                     atatTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATc                                                                   >  1:71246/1‑210 (MQ=255)                                                                                                                                                                     atatTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCaaaaa                          >  1:715644/1‑251 (MQ=255)                                                                                                                                                                     atatTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCaaaaa                          >  1:615051/1‑251 (MQ=255)                                                                                                                                                                       atTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATc                                                                   >  1:403271/1‑208 (MQ=255)                                                                                                                                                                       atTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAAcg                        >  1:201932/1‑251 (MQ=255)                                                                                                                                                                                  aaTTTTGAGGGGAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATc                                                                   <  1:418446/197‑1 (MQ=255)                                                                                                                                                                                             gAAAACCACATGCGAAGTGCCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCGGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGt  >  1:367901/1‑251 (MQ=255)                                                                                                                                                                                                               |                                                                                                                                                                                                                                         AAATTATATATCGCGCAGATTATACAACGATGAGCGCAATATAAAAATATTTTTGAACAAGAAAAGATCACAAGATTGCGTATGTGCCAAACGATTTATCGGCGTTCTTATTTAGAACAAAATATAGTATGATGTGCAACCGATAAACCATGTCATAAAACGATATATTCCACAGCGAATTTTGAGGGGAAAACCACATGCGAAGTACCGAAATTTTACTCACCATCAAATTGAAACAACAACTCTTCGCAGATCCGAAGCGGATCCGTTTATTAAAAGAAATTGAAAAATGCGGCTCCATTAACCAAGCGGCAAAAAATGCCAAAGTCAGTTACAAAAGCGCCTGGGATCACCTTGAAGCCATGAATAAAATCAGCCCGAAACCGTTACTGGAACGCAATATCGGCGGCAAAAACGGCGGCGGCACGGCGCTTACCGT  >  1/821160‑821598

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