Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,512,293	T→C	H41R (CAC→CGC)	fecD ←	ferric citrate ABC transporter membrane subunit FecD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,512,293	0	T	C	100.0%	200.0 / NA	60	H41R (CAC→CGC)	fecD	ferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/33);  total (27/33)

GCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCC  >  NC_000913/4511997‑4512584                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                     gcgcCAGCTTCATCGGCTGGTGGGGCTTTGCCAGCATCTTCAGCAATATCGACCCCGCCATGCCGCCCGCAAATGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATAGCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTc                                                                                                                                                                                                                                                                                                   <  2:488572/300‑1 (MQ=255)      aGCTGCATCGGCTGGTGGGTCTTTGCCAGCATCTTCATCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACATAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGcc                                                                                                                                                                                                                                                                                              <  2:139669/300‑1 (MQ=255)                        gTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCagcagc                                                                                                                                                                                                                                                                            >  1:872147/1‑300 (MQ=255)                            ttGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTTCCTCAATCCATAACAATAATAGCATACGGCATGCCAGTCGGTCAGCAgcgcgc                                                                                                                                                                                                                                                                        >  2:372823/1‑300 (MQ=255)                                           gcaacaccAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGATGCCATGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGa                                                                                                                                                                                                                                                         <  2:346538/293‑1 (MQ=255)                                           gCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGACAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACATAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGa                                                                                                                                                                                                                                                         <  2:42971/300‑1 (MQ=255)                                                 tCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCAGAAGAAGTAGAGCCCCCACAGAGTCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTc                                                                                                                                                                                                                                                   <  2:349851/300‑1 (MQ=255)                                                 tCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCAGAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGTTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTc                                                                                                                                                                                                                                                   <  2:350116/300‑1 (MQ=255)                                                 tCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCACTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTc                                                                                                                                                                                                                                                   <  2:820959/300‑1 (MQ=255)                                                       ccGCCATGCCGTCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGATGCCATGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCatat                                                                                                                                                                                                                                             <  2:151387/300‑1 (MQ=255)                                                           cATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCACCACGTGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAg                                                                                                                                                                                                                                         <  1:195863/300‑1 (MQ=255)                                                                  cccgcAAATGCCAGCAGCGGCAGCACCATCACGTGCGGTGACGGCATAAGAAGTAGAGCCCCCACAGATGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAAc                                                                                                                                                                                                                                  <  2:288705/300‑1 (MQ=255)                                                                  cacgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCNCGCCACGGCGAGGGCTGCACCGACAANCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAAc                                                                                                                                                                                                                                  <  2:580494/298‑1 (MQ=255)                                                                       aaaGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGTTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAgcgc                                                                                                                                                                                                                             <  1:5535/300‑1 (MQ=255)                                                                            ccAGCAGCGGCAGCACCATCACGGGCAGTGACGGCAGAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGc                                                                                                                                                                                                                        <  1:240374/300‑1 (MQ=255)                                                                            ccAGCAACTGCAGCACCATCACGTGCATTTACTGCATAAGAAGTAGGCCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCTGGTTATTCCATAGGGTTGCGCACAATCCCCTGTATCAACACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTACCATCAGCAAGCGCGGCATTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGACAGTCGTTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGcacacc                                                                                                                                                                                                                        <  2:510949/300‑3 (MQ=255)                                                                                       agcaCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGTCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCa                                                                                                                                                                                                             <  2:663166/300‑1 (MQ=255)                                                                                           ccATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAggg                                                                                                                                                                                                         >  1:337319/1‑300 (MQ=255)                                                                                             atcccGGGCAATGACGGCATAAGAAACAGAGCCCCAACCGCAACCCGGCTGACGGCATAGCTAACACCGAGAATATCCGGTGACGCAAGAAGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCCCGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCGGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATCTAGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtg                                                                                                                                                                                                       >  2:775562/5‑300 (MQ=255)                                                                                             aTCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtg                                                                                                                                                                                                       >  1:161839/1‑300 (MQ=255)                                                                                                  gggCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaa                                                                                                                                                                                                  >  1:860054/1‑300 (MQ=255)                                                                                                                gcgagatgtgtataagagacagggCCAGGCTGGCGGCATGGTTCACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCAAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGTGCGCGCCACGGCACGGGGATCACTCACATATGGAGTGATAACAGCGCACAGCATGCCAGGGCAAGGGTGATGAAAATAGCCAGCGCag                                                                                                                                                                                    >  1:263791/23‑299 (MQ=255)                                                                                                                gcgagatgtgtataagagacagggCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCaa                                                                                                                                                                                    >  1:263225/23‑300 (MQ=255)                                                                                                                gcgagatgtgtataagagacagggCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAAGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCaa                                                                                                                                                                                    >  1:263293/23‑300 (MQ=255)                                                                                                                   gcagTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCAATGCACCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACCGTGCCCGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCCCTCCCATATGGAGTGATAACCGCGCACAGCCCGCCAGGGCCAGGGTGATGAAGATCACCAGCGCactgt                                                                                                                                                                                 >  2:18138/3‑296 (MQ=255)                                                                                                                           gCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑Tc                                                                                                                                                                        >  1:190015/1‑300 (MQ=255)                                                                                                                                     agGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTcac                                                                                                                                                              <  1:456380/300‑1 (MQ=255)                                                                                                                                      ggCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTcaca                                                                                                                                                             >  1:808494/1‑300 (MQ=255)                                                                                                                                      ggCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAAGCCACTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACTGACAAACAGTGCCAGCAGCACGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTAAGCAGCGCGCGCCACGGCACGAGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCCGGGCACGGGTGATGACACTAACCAGAGCAATTTTCATCC‑TCGcctactaaca                                                                                                                                                             >  2:465067/1‑293 (MQ=255)                                                                                                                                      ggCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGATGATGCCAGAGGGTTGCGCATAATCCACTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCAACAAACAGTGCCAGCAGCAAACGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGCCCGGCCTGCAAATCGGCCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGGTAACAGCGCATAGCCCGCCCGGGCAAGGTTGATGAAAATAACCAGCGCAATTTTCAtactagcccactaaca                                                                                                                                                              >  2:491691/1‑285 (MQ=255)                                                                                                                                      cccccggCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAAGGTTGAGCACAATCCACTGCATCAGCACGCCCGCCACGGCGAGGGCTGCACCAAAAAACAGCGCAAGCAGCAAGGGCGGCAGTCGGTACTCCATCAATACATAATCATGCTCGCGTACGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGCTCACTCACATATGGAGTGATAACAGCGCACAGCCCGACAGGGCCAGGGTGATGAAAATAACCAGCAGAATTTT‑AGCCATCGCCTCCTcaaa                                                                                                                                                             >  2:176820/6‑298 (MQ=255)                                                                                                                                                       ggTTAACACCGAGAATGTCCGGTGATGCCAGAGGGTTGAGCACAATCCACTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATTAATACATAATAATGTTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCAATATGGAGTGGTAAAAGTGCATGGCCCGCCAGGGCAAGGATGATGAAAATACACAGTGCAATCTTCATTC‑TCGCATCATCACAAGCCAGACAAAGAAAgg                                                                                                                                            >  2:260070/1‑300 (MQ=255)                                                                                                                                                       ggTTAACACCGAGAATATCCGGTGATGCCCGAGGGTTGCGCACAATCCCCTGTATCAGCCCGCCCGCCCCGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCACGCACGGCAGTCGGTACTCCCTCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATAGAGTGATAACAGCGCAAAGACCGCCAGGGCAGGGGTGATGNAAATAACCAGCGCAGTTTGCATCC‑TCGACTACTCAAAAGCCCGACAANGAAAgg                                                                                                                                            >  2:709345/1‑300 (MQ=255)                                                                                                                                                       ggTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAgg                                                                                                                                            >  1:390184/1‑300 (MQ=255)                                                                                                                                                       ggTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCACTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACAGGAAAACAGTGCAAGCAGCAAGCGCGACAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATTGAGTACTAACAGCGCACAGCCCGCCAGGGCAAGGGTGATAACAATACCCAACGCAATTTCCATCC‑TCGCATAATAACAAGCCAGACCAAGCAAgg                                                                                                                                            >  2:147808/1‑300 (MQ=255)                                                                                                                                                       ggTTAACACCGAGAATATCCGGTAATGCCAGAGGGTTGCGCACAATCCACTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCATGCCAGTCGGTCAGTAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGTACAGCCCGCCAGGGCAAGGGTGATGNAAATAACCAGCGCAATTTTCATAC‑TCGCATCCTCACAAGACAGACAANGCAAgg                                                                                                                                            >  2:712243/1‑300 (MQ=255)                                                                                                                                                        gTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAggg                                                                                                                                           <  1:536481/300‑1 (MQ=255)                                                                                                                                                                            gTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCAc                                                                                                                       <  1:217858/300‑1 (MQ=255)                                                                                                                                                                            gTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCAAGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCAc                                                                                                                       <  2:408941/300‑1 (MQ=255)                                                                                                                                                                                                     ccccTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGTGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCgggg                                                                                              <  1:258771/300‑1 (MQ=255)                                                                                                                                                                                                        cTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAg                                                                                           >  1:390381/1‑300 (MQ=255)                                                                                                                                                                                                                gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAgtgc                                                                                   <  2:234501/300‑4 (MQ=255)                                                                                                                                                                                                                gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAgtgc                                                                                   <  1:379996/300‑4 (MQ=255)                                                                                                                                                                                                                            cGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTAc                                                                       <  1:807690/300‑1 (MQ=255)                                                                                                                                                                                                                             ggCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAACGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAACGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTACGGGGAAGGCCAGTGCGCGTGCGAGTACa                                                                      >  1:218497/1‑300 (MQ=255)                                                                                                                                                                                                                                   ggCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCGCGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGcc                                                                <  1:146271/300‑1 (MQ=255)                                                                                                                                                                                                                                                 aCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGt                                                  <  1:457920/300‑1 (MQ=255)                                                                                                                                                                                                                                                      gCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGTTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGccc                                             <  2:826649/300‑1 (MQ=255)                                                                                                                                                                                                                                                           cagcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCcagc                                        <  1:122141/300‑1 (MQ=255)                                                                                                                                                                                                                                                             gcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGANTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCcagcag                                      <  2:768941/300‑1 (MQ=255)                                                                                                                                                                                                                                                             gcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCcagcag                                      <  2:664171/300‑1 (MQ=255)                                                                                                                                                                                                                                                                    cgGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTc                               <  1:768911/300‑1 (MQ=255)                                                                                                                                                                                                                                                                        aGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTg                           >  2:499228/1‑300 (MQ=255)                                                                                                                                                                                                                                                                        aGTCGGTACTCCATCAATACATAATAATGCTCGCGGCCGGGCTGCCAGTCGGGCAGCAGCGCGCGCCACGGCACGGGGATCAGTCCCATATTGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGTGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGTGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTg                           <  2:746607/300‑1 (MQ=255)                                                                                                                                                                                                                                                                             gTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTACGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGt                      >  2:679836/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                        ggagataataATGCTCGCGTCCGGCCTGCCAGTCGGTCGGCAGCGCGCGCCAAGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCATAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAACTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGAAGGGGCTGCCAATCAGCGCCAGCACTACGCCTGTGGGCAGATCTAAGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCAAGCAGCACCAGCGTGGCAACCAGCAGCATGCTGACTGGCAGTACGTTGCGCta           >  2:145918/5‑299 (MQ=255)                                                                                                                                                                                                                                                                                        aTACATAATAATGCTCGCGTCCGGCCTGCCANTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACNAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCACAGCAGCATGCTCACTGGCAGTACGTTGCGCTg           >  2:617383/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                          aCATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGAt         >  1:302701/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                           cATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATc        <  1:348572/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                 aatGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATTGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGcc  <  1:233268/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                     GCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCC‑TCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCC  >  NC_000913/4511997‑4512584

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: ‑