Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,541,058	G→A	T511T (ACC→ACT)	narZ ←	nitrate reductase Z subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,541,058	0	G	A	100.0%	215.7 / NA	62	T511T (ACC→ACT)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/32);  total (30/32)

GCTGTTCATTTGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAG  >  NC_000913/1540770‑1541336                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                       gCTGTTCATTTGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACTAGCCGGTTTGTGGGCGCAGTTTTTACTGGCCGACATAGTGTGCCCAGCCGACACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCGATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCACCGAGCGCCCATGCGTTTCATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATTTACTGGTGAGGCACGCCAGTAGTTTGCTc                                                                                                                                                                                                                                                                             >  2:182565/1‑300 (MQ=255)           tGGCGCGGTGGTCGTTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCt                                                                                                                                                                                                                                                                   <  1:356401/300‑1 (MQ=255)           tGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCNGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCt                                                                                                                                                                                                                                                                   <  1:646667/300‑1 (MQ=255)           tGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCt                                                                                                                                                                                                                                                                   <  1:521264/300‑1 (MQ=255)               gcgGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTgggg                                                                                                                                                                                                                                                               <  1:276151/300‑1 (MQ=255)                  gtggtCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGgtgt                                                                                                                                                                                                                                                            >  1:778012/1‑300 (MQ=255)                    gggcgtttCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATTTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTAc                                                                                                                                                                                                                                                          <  2:567733/294‑1 (MQ=255)                      tCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACgg                                                                                                                                                                                                                                                        <  1:775691/300‑1 (MQ=255)                         gTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGttt                                                                                                                                                                                                                                                     <  2:401429/300‑1 (MQ=255)                         gTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGttt                                                                                                                                                                                                                                                     >  1:39298/1‑300 (MQ=255)                                        aaaGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATc                                                                                                                                                                                                                                      <  1:683997/300‑1 (MQ=255)                                           ggCCATTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCATGCTGGGGTGTACGGTTTGATTTCAGCGTAATCttt                                                                                                                                                                                                                                   <  1:246651/300‑1 (MQ=255)                                           ggCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGCTCATCCCGCGGTAGTTCATGTCCATGTGATAACAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCATTTTACGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGTGGCCCGCCAGTAATTTACTCACCCAAGGCTGGGGTGTACGGTTTGATTCCAGCGtcctcttt                                                                                                                                                                                                                                   >  2:849984/1‑293 (MQ=255)                                           ggCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGCTTTCAGCGTAATCttt                                                                                                                                                                                                                                   >  2:268841/1‑300 (MQ=255)                                           ggCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCttt                                                                                                                                                                                                                                   >  1:701851/1‑300 (MQ=255)                                                gTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCCGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCat                                                                                                                                                                                                                              >  2:427501/1‑299 (MQ=255)                                                          gCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATTTGATACCAGTGGTTAACACCTGCGCCGAGGNTAATCATCGAGCGCCCATGCGTTTTANGGGCGGTATCGGCAAATTCACGTGCGATGGTTTCGATGTACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATc                                                                                                                                                                                                                    <  2:613220/300‑1 (MQ=255)                                                                  gtgGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGcc                                                                                                                                                                                                            <  1:373056/300‑1 (MQ=255)                                                                                      ccGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTc                                                                                                                                                                                        >  1:858189/1‑300 (MQ=255)                                                                                        gACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGc                                                                                                                                                                                      >  1:430722/1‑300 (MQ=255)                                                                                        gACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGc                                                                                                                                                                                      >  2:686599/1‑300 (MQ=255)                                                                                        gACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGACCGCGATCGAGGCAGTAATTCGc                                                                                                                                                                                      >  2:68/1‑300 (MQ=255)                                                                                             aGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTa                                                                                                                                                                                 >  2:197348/1‑300 (MQ=255)                                                                                                gtgCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGTCGTAATTCGCCAATACCa                                                                                                                                                                              >  2:398018/1‑300 (MQ=255)                                                                                                   cccAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAAt                                                                                                                                                                           <  2:391994/300‑1 (MQ=255)                                                                                                                   ttGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCTGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGaccac                                                                                                                                                           <  1:275565/300‑1 (MQ=255)                                                                                                                                ccGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTg                                                                                                                                              >  1:421975/1‑300 (MQ=255)                                                                                                                                             gCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACg                                                                                                                                 >  1:461915/1‑300 (MQ=255)                                                                                                                                                                gTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCaa                                                                                                              >  2:870480/1‑300 (MQ=255)                                                                                                                                                                   gTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTg                                                                                                           >  1:798847/1‑300 (MQ=255)                                                                                                                                                                   gTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGACAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTCACGTCAGGTTTTTAACGTGCAATTg                                                                                                           >  2:391905/1‑300 (MQ=255)                                                                                                                                                                        tGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCACCGAGCGCCCANGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCNCCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCa                                                                                                      >  2:835678/1‑300 (MQ=255)                                                                                                                                                                                tgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCCATTGGCGCACCAGCAcc                                                                                              >  2:542226/1‑300 (MQ=255)                                                                                                                                                                                tgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGATTTCGATATACTGGCGCGGCACGCCAGTAGTTTGCTNACCCCAGGCTGGGGTGTACGGTTTGAGTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCNGCCCGCGATCGAGGCCGTAATTCGCAAGTACCAAATCATAAACGCTGAACACCGGACAGGTGTTCCCATCGACTCACATCAGGTTTTTAACGAGCAATTGGCGCACCAGCAcc                                                                                              >  2:722036/1‑300 (MQ=255)                                                                                                                                                                                    aCCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGAt                                                                                          <  1:304344/300‑1 (MQ=255)                                                                                                                                                                                      cAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGGGGGGCCAGAGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATtg                                                                                        <  1:155007/300‑1 (MQ=255)                                                                                                                                                                                          ggTTAACACCTGCGCCGAGGATAATCATCGANCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCNGGGGTGTACGGTTTGATTTCAGCGTACTCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGtt                                                                                    >  2:628109/1‑300 (MQ=255)                                                                                                                                                                                          ggTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGTGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACTGGCAATTGGCGCACCAGCACCGGATtgtgct                                                                                    >  2:312275/1‑298 (MQ=255)                                                                                                                                                                                           gTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCATCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGttt                                                                                   >  1:33149/1‑300 (MQ=255)                                                                                                                                                                                           gTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAATGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGttt                                                                                   >  2:232035/1‑300 (MQ=255)                                                                                                                                                                                              aaCACCTGCGCCTAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGTGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCATCGGATTGTGTTTTAc                                                                                <  1:290092/300‑1 (MQ=255)                                                                                                                                                                                                cacCTGNGCCGAGGATAATCATCGAGCGCCCATGCGTTGTATGGGCGGTATCGGNAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCATGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGc                                                                              <  2:740911/300‑1 (MQ=255)                                                                                                                                                                                                cacCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGc                                                                              <  1:110533/300‑1 (MQ=255)                                                                                                                                                                                                cacCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGc                                                                              >  2:873209/1‑300 (MQ=255)                                                                                                                                                                                                        ccGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTGTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGANGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCANCAGCACCGGATTGTGTTTTACGCTGCGaaaa                                                                      <  2:616349/300‑1 (MQ=255)                                                                                                                                                                                                                            cccATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATgccgcc                                                  >  2:550195/1‑300 (MQ=255)                                                                                                                                                                                                                                      ttATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAgggg                                        <  1:475170/300‑1 (MQ=255)                                                                                                                                                                                                                                           gtcGGTATCGGCAAGTTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGTCTGGGGTGTACGGTTTGATTTCAGCTTAATCTTTCGCACTGTTTTCATCTTCCAGCCCNCGATCGATGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGc                                   <  2:215815/298‑1 (MQ=255)                                                                                                                                                                                                                                               gTATCGNCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTANTTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACg                               <  2:727478/300‑1 (MQ=255)                                                                                                                                                                                                                                                tATCGGCAAATTCACGGTCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGc                              <  2:812931/300‑1 (MQ=255)                                                                                                                                                                                                                                                 aTCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGcc                             <  1:815603/300‑1 (MQ=255)                                                                                                                                                                                                                                                    ggCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGc                          >  1:717628/1‑300 (MQ=255)                                                                                                                                                                                                                                                       aaaTTCACGGTCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCNAATCATAAACGCTGACCACCGGACAGNTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAAc                       <  2:614272/300‑1 (MQ=255)                                                                                                                                                                                                                                                       aaaTTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAAc                       <  1:259441/300‑1 (MQ=255)                                                                                                                                                                                                                                                         gttCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCGCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAg                     <  2:295385/299‑1 (MQ=255)                                                                                                                                                                                                                                                                ggCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATg              >  1:765954/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 gCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTNGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGt             <  2:667763/300‑1 (MQ=255)                                                                                                                                                                                                                                                                   gATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTg           <  2:696277/300‑1 (MQ=255)                                                                                                                                                                                                                                                                      ggTTTCGATATACTGGCGCGGCACGCTAGTAATTTGCTCACCCCAGGCTGGGGGGTACGGTTTGATTTCAGCTTAATCTTTCACACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGAcc        <  2:306234/300‑1 (MQ=255)                                                                                                                                                                                                                                                                      ggTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGAcc        <  1:302664/300‑1 (MQ=255)                                                                                                                                                                                                                                                                      ggTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGACACGACTGCAACAGCGTCATGTTGAcc        >  2:533976/1‑300 (MQ=255)                                                                                                                                                                                                                                                                            gatatACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAg  <  1:333280/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                       GCTGTTCATTTGGCGCGGTGGTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAG  >  NC_000913/1540770‑1541336

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