Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,512,501	C→T	A303A (GCG→GCA)	fecC ←	ferric citrate ABC transporter membrane subunit FecC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,512,501	0	C	T	100.0%	197.9 / NA	65	A303A (GCG→GCA)	fecC	ferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/37);  total (28/37)

GGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGT  >  NC_000913/4512257‑4512745                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      ggCAGTCGGTACTCCATCAATCCATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCACGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATACCCAGCGCAATTTTCATCCTCGCCTCCTCACACGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGt                                                                                                                                                                                                                                                 >  2:125286/1‑251 (MQ=255)              ccATCAATACATAATAATGCTCGCGCCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAATCAAGGGCTGCCAATCAGCGCCAGCACTTCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCt                                                                                                                                                                                                                                    <  1:179774/250‑1 (MQ=255)              ccATCAACACATAATAACGCTCGCGTCCGGCCTGCCAGTCGGTCACCAGCGCGCGCCACGGCACGGGGATCACCCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGTGGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCTCAAGCCAGACAAAGCAAGGGCTGACCATCAGCGCCAGCACTGCGCCTGCGTGCGGATCTCCGGGGAAGGCCAGTGCGCGTGTG‑AGTACATCt                                                                                                                                                                                                                                    <  2:192365/250‑1 (MQ=255)                  cAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCca                                                                                                                                                                                                                                <  1:101267/250‑1 (MQ=255)                    aTACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCcagca                                                                                                                                                                                                                             <  1:239488/251‑1 (MQ=255)                      aCATAATAATGCTCGCGTCCGGCCTGCCCGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCagcagc                                                                                                                                                                                                                           <  1:197489/251‑1 (MQ=255)                           ataatGCTCGCGTCCGGCCGGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCTCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCcatcagcatcag                                                                                                                                                                                                                      <  2:251109/251‑1 (MQ=255)                                 cTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGc                                                                                                                                                                                                                >  1:210442/1‑251 (MQ=255)                                        ccGGCCTGCCAGGGGGTCAGCAGCGCGCGCCACGGCAGGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCGTCGCCTCCTCACAAGCCAGGCAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGTCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCcag                                                                                                                                                                                                          <  1:159626/250‑1 (MQ=255)                                                ccAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATGACGCCCATATGGAGTGATAACAGCGCATAGCCCGCCAGTGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCGCCGGGGAAGGCCATTGCTCGTGCG‑AGTACATCTGCCAGCAGCATGAGCATGGCCCCCAGGGGCATGCt                                                                                                                                                                                                 <  2:32564/251‑1 (MQ=255)                                                   gTCGGACAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCAc                                                                                                                                                                                              >  2:75104/1‑251 (MQ=255)                                                      ggTCAGCAGCGCGCGCCACGGCACGGGGATCCCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCATCAGCATCAGCGTGGCCCACAGCAGAATGCTCAATgg                                                                                                                                                                                           <  2:61838/251‑1 (MQ=255)                                                       gTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGc                                                                                                                                                                                          <  1:27622/251‑1 (MQ=255)                                                                   cgcCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTgcgc                                                                                                                                                                              >  2:133205/1‑251 (MQ=255)                                                                       ccggcacggGGATCACTCCCATAGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCGCACTGGCAGTACGTTGCGCTGAt                                                                                                                                                                          <  2:58109/250‑1 (MQ=255)                                                                        cggcacggGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATc                                                                                                                                                                         <  2:203689/251‑1 (MQ=255)                                                                              ggggATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGCTGAAAATAACCAGCGCACTTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTCCATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGcc                                                                                                                                                                   >  1:240478/1‑251 (MQ=255)                                                                                  accACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGTCAATCAGCGCCAGCCCTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTCCGTTGAGCTGATCGAAGCCTGcc                                                                                                                                                               <  1:56111/249‑1 (MQ=255)                                                                                    cACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGccc                                                                                                                                                              >  2:214170/1‑250 (MQ=255)                                                                                     aCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATACTAGCCTCCTCACAAGCCAGACACAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGccc                                                                                                                                                              >  1:39955/1‑249 (MQ=255)                                                                                        cccATAGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGaa                                                                                                                                                          <  1:54069/250‑1 (MQ=255)                                                                                             atGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAgcgcg                                                                                                                                                     <  1:63956/250‑1 (MQ=255)                                                                                                 ggTGATACCAGCGCCAAGCCCGCCAGGGCAAGGGGGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGa                                                                                                                                                <  2:149937/250‑1 (MQ=255)                                                                                                    gATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGAtgt                                                                                                                                             <  2:257250/251‑1 (MQ=255)                                                                                                      tAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTgg                                                                                                                                           >  1:98851/1‑251 (MQ=255)                                                                                                        aCAGCGCCCAGCCCGCCAGGGCAAGCGTGATGAAAATAACAAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCACGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGc                                                                                                                                          >  2:265133/1‑250 (MQ=255)                                                                                                          aGCGCACAGCACGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGACAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCAATGCGACTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAACGCGCCAGATGTGGCAcc                                                                                                                                       >  2:125528/1‑251 (MQ=255)                                                                                                             gcACAGCCCGCGAGGGCAAGGGTGATGAACATAACCAGCGCAATTTTCACCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGGCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCTCTGGCAGTACGTGGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACcagc                                                                                                                                    <  1:36978/251‑1 (MQ=255)                                                                                                                 aGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGAc                                                                                                                                <  1:230692/251‑1 (MQ=255)                                                                                                                      gCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGAt                                                                                                                            >  1:210118/1‑250 (MQ=255)                                                                                                                           ggCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGc                                                                                                                      <  2:210442/251‑1 (MQ=255)                                                                                                                                    gatgaAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACGGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGAcc                                                                                                             <  2:15786/251‑1 (MQ=255)                                                                                                                                     aggaAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCATCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCTGACCt                                                                                                            <  2:183096/249‑1 (MQ=255)                                                                                                                                          aaTAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCGCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCg                                                                                                                                                                        <  1:182802/186‑1 (MQ=255)                                                                                                                                          aaTAAACAGCGCACTTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCg                                                                                                                                                                        >  2:182802/1‑186 (MQ=255)                                                                                                                                                 aGCGCCATTTTCCTCCTCGCCTCCTCTCAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCGCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCATCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCCCCAGCAGACCGATAACCGCTTCCGGACCTGCCTCGCTGacg                                                                                                <  2:214652/251‑1 (MQ=255)                                                                                                                                                               ccTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCaccaacca                                                                                   <  1:87561/250‑1 (MQ=255)                                                                                                                                                                  cGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGa                                                                                 >  1:224572/1‑249 (MQ=255)                                                                                                                                                                    cctcctCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAagc                                                                              >  1:94214/1‑250 (MQ=255)                                                                                                                                                                       cctcACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAagcagcc                                                                          >  1:101190/1‑250 (MQ=255)                                                                                                                                                                         ccacaAGCCAGAAAAATCAAGGGCCTCCAAGCAGCGCCAGCACGGCGCCTGCGGGCAGATCCCCGGGGAAGGCCAGTGCGCG‑GCGGAGTCCATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCGCCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACt                                                                        <  2:118847/250‑1 (MQ=255)                                                                                                                                                                          cccaAGCCAGACAAAGCAAGGGCTGCCAACCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACGGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCAc                                                                         <  2:180205/247‑1 (MQ=255)                                                                                                                                                                                   gACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCGCTAACAGATtga                                                              <  1:226145/251‑1 (MQ=255)                                                                                                                                                                                    aCAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgat                                                             <  1:133205/251‑1 (MQ=255)                                                                                                                                                                                     aaaaGCAAGGGCTGCCCATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTCACAAATtgatg                                                            >  2:809/2‑251 (MQ=255)                                                                                                                                                                                      aaaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgatgc                                                           >  1:139916/1‑250 (MQ=255)                                                                                                                                                                                      aaaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgatg                                                            >  1:83656/1‑250 (MQ=255)                                                                                                                                                                                      aaaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACAAGAAGCAGCACTAACATATtgatg                                                            >  2:90797/1‑250 (MQ=255)                                                                                                                                                                                      aaaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGAGAAACGCCACCGGTCCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgatga                                                           <  1:265093/251‑1 (MQ=255)                                                                                                                                                                                           aaGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACCTCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCCCTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCCGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATCTTGATGACCaa                                                       >  1:20579/1‑250 (MQ=255)                                                                                                                                                                                               gCTGCCAAGCAGCGCCAGCCCTGGGCCGGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCAGCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTa                                                  <  2:79024/251‑1 (MQ=255)                                                                                                                                                                                                      aTCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACCTATTGATGACCAAACGTAGCCt                                              >  2:169856/1‑248 (MQ=255)                                                                                                                                                                                                      aTCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCt                                              <  1:169856/248‑1 (MQ=255)                                                                                                                                                                                                        cAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTc                                          >  1:63846/1‑250 (MQ=255)                                                                                                                                                                                                        cAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTc                                          >  1:14165/1‑250 (MQ=255)                                                                                                                                                                                                             gcAGCACTGCGCCGGAGGGCAGATCTCCGGGGAAGGCCCGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCCAAATGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGAGGCACCAGCAGACCGATAAACGCCACCTGACCTGCAAAGCTGACGCACGCACAAACCAGAAGCGGCACTAACATATTGATGACCAAACGAAGCCACGTCCGGttt                                    >  2:164866/2‑250 (MQ=255)                                                                                                                                                                                                                    tGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTccc                              >  1:28557/1‑250 (MQ=255)                                                                                                                                                                                                                    ggcgcCTGCGGGCAGATCTCCGGGGAAGGCCGGTGCGCGGGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACGGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCa                             <  2:14165/250‑1 (MQ=255)                                                                                                                                                                                                                    ggcgcCAGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCa                             <  2:66426/250‑1 (MQ=255)                                                                                                                                                                                                                       gcCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCg                          <  2:210118/251‑1 (MQ=255)                                                                                                                                                                                                                         cTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGt                         >  1:104308/1‑250 (MQ=255)                                                                                                                                                                                                                                           ggAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTg      >  1:192851/1‑251 (MQ=255)                                                                                                                                                                                                                                            gAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACTCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGa     <  1:182766/251‑1 (MQ=255)                                                                                                                                                                                                                                             aaGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGAGCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCTCTGAg    <  1:214170/251‑1 (MQ=255)                                                                                                                                                                                                                                               ggCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGt  >  1:246877/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      GGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGT  >  NC_000913/4512257‑4512745

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