Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,658,974	A→G	Q302Q (CAA→CAG)	ynfE →	putative selenate reductase YnfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,658,974	0	A	G	100.0%	194.7 / NA	52	Q302Q (CAA→CAG)	ynfE	putative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/29);  total (23/29)

GCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGC  >  NC_000913/1658781‑1659193                                                                                                                                                                                                  |                                                                                                                                                                                                                            gCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCg                                                                                                                                                                                                       <  1:499064/216‑1 (MQ=255)       gtggCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCa                                                                                                                                                                          >  1:513201/1‑239 (MQ=255)       gtggCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAg                                                                                                                                                                         >  1:722576/1‑240 (MQ=255)         ggCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCt                                                                                                                                                                                  <  2:251460/229‑1 (MQ=255)             tCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTAccc                                                                                                                                                                             >  2:364701/1‑230 (MQ=255)               aCTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTtc                                                                                                                                                                                                                      <  1:606964/187‑1 (MQ=255)               aCTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTtc                                                                                                                                                                                                                      >  2:606954/1‑187 (MQ=255)               aCTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTAccc                                                                                                                                                                             >  2:179185/1‑228 (MQ=255)               aCTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCAccc                                                                                                                                                                 >  1:520753/1‑240 (MQ=255)                       tcttGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATgg                                                                                                                                                         <  2:52664/240‑1 (MQ=255)                        cttGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGt                                                                                                                                                        <  1:598561/240‑1 (MQ=255)                         ttGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTc                                                                                                                                                       <  1:160861/240‑1 (MQ=255)                         ttGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTc                                                                                                                                                       <  1:133147/240‑1 (MQ=255)                           gAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCAc                                                                                                                                                     >  2:498027/1‑240 (MQ=255)                                aGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGt                                                                                                                                                        >  2:372112/1‑232 (MQ=255)                                       gagaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCta                                                                                                                                          <  1:644433/239‑1 (MQ=255)                                          aaaTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATAtt                                                                                                                                      <  1:234910/240‑1 (MQ=255)                                                  cGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATg                                                                                                                                                                      <  2:166149/200‑1 (MQ=255)                                                      aaaaTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGt                                                                                                                          <  2:238041/240‑1 (MQ=255)                                                              tGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATaa                                                                                                                  >  2:723874/1‑240 (MQ=255)                                                                ttATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTA‑CCGCAGATGCACCCAAAAATGGTCACTATaa                                                                                                                                                <  2:79679/207‑1 (MQ=255)                                                                   tCGACCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCTGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGt                                                                                                                          <  2:21477/227‑1 (MQ=255)                                                                          gcgATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGt                                                                                                                          <  2:355508/220‑1 (MQ=255)                                                                             atatACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAAcg                                                                                                     <  2:620439/238‑1 (MQ=255)                                                                                      tACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGc                                                                                          >  1:722590/1‑240 (MQ=255)                                                                                        cGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCtt                                                                                        >  1:6328/1‑240 (MQ=255)                                                                                          gCTGCCGGTCGTGATGACGAATGGCTCCCTATTCGCCCGGACACCGATGCCGCGCTGGTTGCGGGTAGTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGCTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAACCAGCGAAAACGCCGCAGTGGGCTTc                                                                                       >  1:301238/1‑239 (MQ=255)                                                                                             gCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCaa                                                                                   >  2:527290/1‑240 (MQ=255)                                                                                              ccGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGt                                                                                                                                                        >  1:277861/1‑170 (MQ=255)                                                                                              ccGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGt                                                                                                                                                        <  2:277858/170‑1 (MQ=255)                                                                                                  tCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAAcg                                                                                                     <  1:305116/217‑1 (MQ=255)                                                                                                   cGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTAc                                                                              <  2:334369/239‑1 (MQ=255)                                                                                                             aaTGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCgg                                                                   <  1:271590/240‑1 (MQ=255)                                                                                                             aaTGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCg                                                                    <  1:675623/239‑1 (MQ=255)                                                                                                              aNGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTc                                                                                                                                                       <  2:35328/155‑1 (MQ=255)                                                                                                                  cTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCATCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCggagga                                                               <  2:124177/239‑1 (MQ=255)                                                                                                                    cccTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCg                                                            <  2:174194/240‑1 (MQ=255)                                                                                                                            tcccGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGCTGCACCCAAAAAt                                                                                                                                                           >  2:507484/2‑137 (MQ=255)                                                                                                                             cccGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGt                                                                                              >  2:449472/1‑197 (MQ=255)                                                                                                                                 ggCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTg                                               >  1:614/1‑240 (MQ=255)                                                                                                                                                 tGGTTGCGGGTATTGCCTGTGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGc                                <  2:7238/239‑1 (MQ=255)                                                                                                                                                  ggTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGcac                              <  2:570119/240‑1 (MQ=255)                                                                                                                                                   gTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGcac                              >  1:712729/1‑239 (MQ=255)                                                                                                                                                         ggTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATaa                                                                                                                                                <  1:20850/119‑1 (MQ=255)                                                                                                                                                         ggTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATaa                                                                                                                                                >  2:20850/1‑119 (MQ=255)                                                                                                                                                          gTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAAc                      <  1:591985/240‑1 (MQ=255)                                                                                                                                                                          gATTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAAcc                     >  1:284789/1‑225 (MQ=255)                                                                                                                                                                           aTTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAg     <  1:403576/240‑1 (MQ=255)                                                                                                                                                                           aTTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAg     >  1:606912/1‑240 (MQ=255)                                                                                                                                                                           aTTAACGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAg     >  1:238691/1‑240 (MQ=255)                                                                                                                                                                              aaCGAAAATCTCGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGc  >  1:695828/1‑240 (MQ=255)                                                                                                                                                                                       ctcGTTGATCAGCCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGc                  <  2:722710/215‑1 (MQ=255)                                                                                                                                                                                                  |                                                                                                                                                                                                                            GCGGTGGTGGCATCACTTATCTTCTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGC  >  NC_000913/1658781‑1659193

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