Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	448,621	G→A	H14Y (CAC→TAC)	cyoC ←	cytochrome bo3 subunit 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	448,621	0	G	A	100.0%	154.4 / NA	45	H14Y (CAC→TAC)	cyoC	cytochrome bo3 subunit 3
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/31);  total (14/31)

GTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACG  >  NC_000913/448388‑448863                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   gTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTcgcg                                                                                                                                                                                                                                    <  1:39800/250‑1 (MQ=255)      gggCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTc                                                                                                                                                                                                                              <  1:186082/250‑1 (MQ=255)       ggaagtagcgcccATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTGTCAACCAGAACTAACGGCAGTTCGAAAAGGTCCTTACCTGGCGGGCCGCCGGCGGTGCCGTTCGCCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAAAGCAGTCGCTCATCAGGTAGATCCAAAAACCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGTGTGGGCAGTCGCGTGCGTCa                                                                                                                                                                                                                             <  2:74486/240‑1 (MQ=255)           aTAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGt                                                                                                                                                                                                                         <  1:367257/251‑1 (MQ=255)                  ccATGCCGTAGGTGATGGAGCCGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGtt                                                                                                                                                                                                                  <  1:433367/251‑1 (MQ=255)                                    aGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGCAAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGCCAAAGTATCAGTTGCCATTTTTCAGCCCTGc                                                                                                                                                                                                >  1:156365/1‑251 (MQ=255)                                      cTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGcc                                                                                                                                                                                               >  1:389664/1‑250 (MQ=255)                                        gAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGcc                                                                                                                                                                                               <  2:389664/248‑1 (MQ=255)                                        gAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCtt                                                                                                                                                                                             <  1:155953/250‑1 (MQ=255)                                        gAACAAAAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCATTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGTTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCtt                                                                                                                                                                                             <  2:394817/250‑1 (MQ=255)                                         aacaacAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCCTTTTTCAGCCCTGCCtt                                                                                                                                                                                             >  2:48886/1‑249 (MQ=255)                                            aacaGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTaa                                                                                                                                                                                        <  2:393112/251‑1 (MQ=255)                                                gCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAAtctc                                                                                                                                                                                    <  1:166766/251‑1 (MQ=255)                                                     aaaGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTGACCTGTCGTGCCGCCGGCGGTGCCGGTCACCAGAACGGCATAGTTAGCAAACAAGATAGATAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCTTACGCGTGGGCAGTCGCGTTCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGa                                                                                                                                                                               <  2:347302/251‑1 (MQ=255)                                                          ttCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTGCCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATg                                                                                                                                                                           <  1:277543/250‑1 (MQ=255)                                                                    acgaacgGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCa                                                                                                                                                                <  1:152652/251‑1 (MQ=255)                                                                        acgGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGGATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAACTGCTGGTTTTCCAGttt                                                                                                                                                            >  2:326306/1‑251 (MQ=255)                                                                                tCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGAtt                                                                                                                                                     <  1:326567/250‑1 (MQ=255)                                                                                             tACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGc                                                                                                                                         <  1:193408/249‑1 (MQ=255)                                                                                                    cGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACgtagta                                                                                                                                 <  1:80624/250‑1 (MQ=255)                                                                                                     gggCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGGGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAGGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTATTAATCTCATCGAAATGCTGGGTTTTCAGTTTGTCGATTTCTGCCGCCGGCACGTAGTAAt                                                                                                                               <  2:40305/251‑1 (MQ=255)                                                                                                     gggCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAAt                                                                                                                               >  2:56659/1‑251 (MQ=255)                                                                                                       gccgccTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATcc                                                                                                                             <  1:56659/251‑1 (MQ=255)                                                                                                       gccgccTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATcc                                                                                                                             <  1:247961/251‑1 (MQ=255)                                                                                                        ccgccTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCa                                                                                                                            <  1:307410/251‑1 (MQ=255)                                                                                                         cgccTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAGCAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCAc                                                                                                                           <  2:166206/251‑1 (MQ=255)                                                                                                         cgccTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCAc                                                                                                                           >  1:62716/1‑251 (MQ=255)                                                                                                               cGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCtc                                                                                                                     <  1:52952/251‑1 (MQ=255)                                                                                                                         cACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCt                                                                                                            <  1:424537/250‑1 (MQ=255)                                                                                                                                  ggCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCATCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTcc                                                                                                                                                                 <  1:151908/188‑1 (MQ=255)                                                                                                                                  ggCATAGGTAGCAAACAAGATAGAGAACAAAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTcc                                                                                                                                                                 >  2:151908/1‑188 (MQ=255)                                                                                                                                         gTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGtga                                                                                           <  2:427937/251‑1 (MQ=255)                                                                                                                                          aaGCAAACAAGATAGAGAACAGAATGCAGCCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGtgat                                                                                          >  2:319791/2‑251 (MQ=255)                                                                                                                                              aaaCAAGTTAGATAACAGAATGCAGTCGCTCATTAGGTAGATCCAAAATCCGAAGAGTTTGGATCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCTTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAAGCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATc                                                                                      <  2:75089/251‑1 (MQ=255)                                                                                                                                                     aTAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTg                                                                               >  2:127915/1‑251 (MQ=255)                                                                                                                                                        gagaACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGa                                                                            <  2:259637/251‑1 (MQ=255)                                                                                                                                                                    aaGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTGTGTTTCCGGCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGACATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCg                                                                <  1:230525/250‑1 (MQ=255)                                                                                                                                                                     aGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCACATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCAc                                                                                                                           <  1:201161/191‑1 (MQ=255)                                                                                                                                                                     aGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGAGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCAc                                                                                                                           >  2:201161/1‑191 (MQ=255)                                                                                                                                                                            catcaGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGccac                                                        >  1:322665/1‑251 (MQ=255)                                                                                                                                                                                            aaTCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGACAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAAAAATCGCCAGCCACCAGATCTGCCAGATc                                         >  1:192127/1‑250 (MQ=255)                                                                                                                                                                                                              ccGCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGAtgg                      <  1:380259/251‑1 (MQ=255)                                                                                                                                                                                                                gCCTGCATCGTGGTGCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTGTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGAGCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGAtggtg                    <  2:129270/251‑1 (MQ=255)                                                                                                                                                                                                                            gtgCCCGTGTTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGc         >  1:91145/1‑250 (MQ=255)                                                                                                                                                                                                                                  gtgtTCGTACGCGTGGGCAGTCGCGTGCGTCAAAGTATCCGTTGCCATTTTTCAGACCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACg  >  2:275370/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   GTACATGGCGATAGCCGCCATGCCGTAGGTGATGGAGCTGAACAACAGCAAGAAAGTTTCAACCAGAACGAACGGCAGTTCGAAAATGTCCTTACCTGTCGGGCCGCCTGCGGTGCCGTTCACCAGAACGGCATAGGTAGCAAACAAGATAGAGAACAGAATGCAGTCGCTCATCAGGTAGATCCAAAATCCGAAGATTTTGGTTCCGCCTGCATCGTGGTGCCCGTGTTCGTGCGCGTGGGCAGTCGCGTGCGTCAAAGTATCAGTTGCCATTTTTCAGCCCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACG  >  NC_000913/448388‑448863

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