Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,890,626	G→A	L425F (CTC→TTC)	cysJ ←	sulfite reductase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,890,626	0	G	A	100.0%	183.8 / NA	59	L425F (CTC→TTC)	cysJ	sulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/26);  total (33/26)

TCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAAT  >  NC_000913/2890337‑2890917                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                    tCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGa                                                                                                                                                                                                                                                                                            >  2:310671/1‑299 (MQ=255)  cgcgcgACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAGGCTGGa                                                                                                                                                                                                                                                                                            >  1:227684/1‑298 (MQ=255)      cgACCAAGCAAGATCGATACGTGTCAGTACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGGAGTCTTACGTAAAGTGCGGATTACCAAAGAACAGCAAGTTTTTACCTAGCGCTTCGTCGGCGGCGCGTTGCTGCATACAGGCGCGGACCGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGACGAACAAGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGcccc                                                                                                                                                                                                                                                                                      >  2:25150/1‑297 (MQ=255)             gCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCNGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCa                                                                                                                                                                                                                                                                               <  2:517793/300‑1 (MQ=255)             gCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCa                                                                                                                                                                                                                                                                               <  1:523116/300‑1 (MQ=255)                               aGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGAGGCCTTCCAc                                                                                                                                                                                                                                                              <  1:233106/299‑1 (MQ=255)                                gCACGCCATCTTTGACTTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcgt                                                                                                                                                                                                                                                            <  2:668182/300‑1 (MQ=255)                                       aTCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATGAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAAGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCAAGGGCGCGGCCTTCCAcgtcgt                                                                                                                                                                                                                                                         >  2:354720/1‑296 (MQ=255)                                       aTCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATCCCGGTGCCTGGGCCAATCATAATCAACGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCAACACCGGCACGGGCGCGGCCTTCCACGTAGTAAc                                                                                                                                                                                                                                                      >  2:166961/1‑299 (MQ=255)                                        tCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGTAAGTCTTCCGTAAAGTGCGGACTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAATGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGc                                                                                                                                                                                                                                                    <  2:259808/300‑1 (MQ=255)                                        tCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTACACGTCGTAACGc                                                                                                                                                                                                                                                    <  1:475798/300‑1 (MQ=255)                                         cTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGca                                                                                                                                                                                                                                                   >  1:609696/1‑300 (MQ=255)                                          tttGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTNCCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcac                                                                                                                                                                                                                                                  <  1:495123/300‑1 (MQ=255)                                              aCGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccaca                                                                                                                                                                                                                                              >  2:322433/1‑300 (MQ=255)                                                 tggcgcTGCCACTCCACCTGGTACAGTAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCTGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccacacc                                                                                                                                                                                                                                            <  2:329433/297‑1 (MQ=255)                                                     gcTGCCACTCCACCTGGTACAGGAATTCTTCCTTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccacacc                                                                                                                                                                                                                                            <  2:90585/295‑1 (MQ=255)                                                                aCCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACt                                                                                                                                                                                                                            >  1:636194/1‑300 (MQ=255)                                                                aCCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACt                                                                                                                                                                                                                            >  2:402276/1‑300 (MQ=255)                                                                                        aaGTGTGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGTGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCAACGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGCCTTCCGc                                                                                                                                                                                                           >  2:19640/1‑293 (MQ=255)                                                                                               gATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCAACGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGa                                                                                                                                                                                             >  2:631147/1‑300 (MQ=255)                                                                                                       aaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGTTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATAc                                                                                                                                                                                     <  2:400277/300‑1 (MQ=255)                                                                                                          aaccgCCAGTTTTTACCTTGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGTGCCAATCATAATCACCTGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAACTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGTTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGa                                                                                                                                                                                  <  2:129101/296‑1 (MQ=255)                                                                                                          aaCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGTAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGa                                                                                                                                                                                  <  1:264902/300‑1 (MQ=255)                                                                                                          aaCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGa                                                                                                                                                                                  <  1:265175/300‑1 (MQ=255)                                                                                                          aaCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGACTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACTTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACGTCCGCCTGCGAGGAGGCGATGGAATACAGa                                                                                                                                                                                  <  1:659007/300‑1 (MQ=255)                                                                                                          aaCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGACTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGa                                                                                                                                                                                  <  1:659000/300‑1 (MQ=255)                                                                                                             aGCCAGTTTTTACCTGGCGCTTCGTCGGCTGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCTGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAAGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcg                                                                                                                                                                                <  1:348918/299‑1 (MQ=255)                                                                                                             aGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGTCGGAAGTTATCGTTATGTTCGATAAATAAGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcg                                                                                                                                                                                <  1:348510/299‑1 (MQ=255)                                                                                                             aGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgc                                                                                                                                                                               >  2:685400/1‑300 (MQ=255)                                                                                                                       tACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCg                                                                                                                                                                     <  1:98802/300‑1 (MQ=255)                                                                                                                                                ttGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTc                                                                                                                                            >  1:118748/1‑300 (MQ=255)                                                                                                                                                     gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCAt                                                                                                                                       >  1:637023/1‑300 (MQ=255)                                                                                                                                                                       gcgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAAGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGcc                                                                                                                            >  2:607518/1‑293 (MQ=255)                                                                                                                                                                       gcgcAATACCGGTGCCTGGGCCAATCATAATCAACGGGGTTTATGGATTGGCTGGCAGGGGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTACTTAACGCTTCGGCATCAAGCTGTGc                                                                                                                             >  2:193675/1‑292 (MQ=255)                                                                                                                                                                         gcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGTAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGaaa                                                                                                                   <  1:607518/300‑1 (MQ=255)                                                                                                                                                                                  gTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGccacat                                                                                                           >  2:87892/1‑294 (MQ=255)                                                                                                                                                                                             actcctAATCCCCGGGGTTTCTGGATTGGCTGGCCGACGGAAGTTATCGTTATGTTCGATAAATACGCGGGCTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTAGTAACGCACCACACCAACGATGACGTGTACTTCGTTCTCGAATTCCACCTGCGAGGAGGCGATGGAATACAGACGCGGCGTAAGCGGGCGCAGCAGATTAATTAGCGATTCGGCATCAAGCTGTGCAGGGGAGAAACGCAAAATGTCAACAATc                                                                                                 >  2:277737/6‑298 (MQ=255)                                                                                                                                                                                             actCATAATCAACGGGATTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAACTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTATACTTCGTTCTCGCCTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGACGCAGCAGcttaatt                                                                                                                                                    >  2:277683/3‑247 (MQ=255)                                                                                                                                                                                              aTCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCAGGGGAGAAACGCACCATGTCAACAATCGGc                                                                                              >  1:353274/1‑300 (MQ=255)                                                                                                                                                                                              aTCATAATCAACGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGAGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCACGTCAACAATCGGc                                                                                              >  2:143204/1‑300 (MQ=255)                                                                                                                                                                                                              tttCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAAt                                                                              >  1:561317/1‑300 (MQ=255)                                                                                                                                                                                                                tCTGGATTGGCTGGCAGGCGGAAGTGATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGTAGGCACCACCGGCACGGGCGCGGCCTTCCACGGCGTAACGCACCACACCAACGTTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGc                                                                            <  1:470578/300‑1 (MQ=255)                                                                                                                                                                                                                        ggctggcAGGCGGAAATTATCGTTATGTTCGATAAATACGTGGACTTCGCCCTCTTCTTCCCCGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGAGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCAGGGGAGAAGCGCACCATGTCAACAATCGGCGTCGCCGAGGTGTAATGCTGTAACtt                                                                    >  2:128943/1‑300 (MQ=255)                                                                                                                                                                                                                                     aaGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCAACGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGccc                                                       >  1:279293/1‑300 (MQ=255)                                                                                                                                                                                                                                             gTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTAGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAAAAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCccacca                                                   >  2:308710/1‑296 (MQ=255)                                                                                                                                                                                                                                                     cGATAAATCCGCGGACTTCGCCCTCTTCTCACACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCAACGTCGTACCGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGAGATGGAATACAGACGCGGCGTCAGCGGGTGCAGCAGATTAATTAGCGCTTAGGCATCAAGTTGTGGCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGTTGTAACttagatttctcgaacacaagtggcggccgt                                        >  2:277815/1‑271 (MQ=255)                                                                                                                                                                                                                                                       aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAgtgt                                      <  1:68842/299‑1 (MQ=255)                                                                                                                                                                                                                                                        tAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCTAAGAAGCTGGATGCACCACCGGCGCGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAACGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGttt                                    <  2:486437/300‑1 (MQ=255)                                                                                                                                                                                                                                                           aTACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCAc                                 <  1:530891/300‑1 (MQ=255)                                                                                                                                                                                                                                                            tACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGANGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTAGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACt                                >  2:513493/1‑300 (MQ=255)                                                                                                                                                                                                                                                                        gCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCg                     <  1:394090/299‑1 (MQ=255)                                                                                                                                                                                                                                                                        gCCCTCGTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCTGCACGGGAGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATNGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCg                     <  2:228405/299‑1 (MQ=255)                                                                                                                                                                                                                                                                           ctctTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGc                 >  1:126610/1‑300 (MQ=255)                                                                                                                                                                                                                                                                           ctctTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGc                 >  2:491925/1‑300 (MQ=255)                                                                                                                                                                                                                                                                           ctctTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTAGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGt                    >  2:445973/1‑297 (MQ=255)                                                                                                                                                                                                                                                                             cttcttCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCAGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTAGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGATTTATCGCCCACCAACGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGt               >  2:58886/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                       gcgGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTcaac     <  1:469980/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                          gTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTNATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAAt  >  1:504772/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                          gTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATAGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAAt  >  2:159965/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                    TCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAAT  >  NC_000913/2890337‑2890917

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