Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,515,348	A→C	L444V (TTG→GTG)	fecA ←	ferric citrate outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,515,348	0	A	C	100.0%	165.9 / NA	49	L444V (TTG→GTG)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/23);  total (26/23)

GGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGG  >  NC_000913/4515101‑4515586                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               ggAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGGTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTGGATTCATTCACATa                                                                                                                                                                                                                                             <  1:54423/251‑1 (MQ=255)    aGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCg                                                                                                                                                                                                                                          <  1:267581/251‑1 (MQ=255)     gctgTGCGTTATAGCTCTCTTCTTCGTGCGTGCCTGGGATGGCGTTGTTCTTGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCTATGTCGATTTTGTCATCCAGATACCAGTCGTGCGCCTCGGTTCCGGATCTCGTATCGCGGTCGTAATGGCTTTTGCCGGCCGGCAACTGCTCGCTTCTGGTGTCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGa                                                                                                                                                                                                                                         <  1:207911/248‑1 (MQ=255)       ggTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATa                                                                                                                                                                                                                                       <  1:239084/251‑1 (MQ=255)          gCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGGAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTTGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAgc                                                                                                                                                                                                                                     <  1:140264/250‑1 (MQ=255)              tataGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAgcccacg                                                                                                                                                                                                                                >  1:96186/1‑251 (MQ=255)               ataGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAgcccacgc                                                                                                                                                                                                                               >  1:241656/1‑251 (MQ=255)                 aGCTCACTTCTGCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGGCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAgcccacgccc                                                                                                                                                                                                                             <  1:103276/251‑1 (MQ=255)                     cACTTCTTCGTGCGTGCCTGTGATGGCGTTGGTCTGGTATGACTCGATATGGTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAgcccacgcccac                                                                                                                                                                                                                           <  1:118528/249‑1 (MQ=255)                            tcGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCACGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTTAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTg                                                                                                                                                                                                                     <  1:234920/248‑1 (MQ=255)                             cgtgcgtgCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGcgcg                                                                                                                                                                                                                 >  1:92547/1‑251 (MQ=255)                              gtgcgtgcCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCgg                                                                                                                                                                                                                <  1:166762/251‑1 (MQ=255)                                     ccTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGCAGGGcc                                                                                                                                                                                                         >  1:122882/1‑251 (MQ=255)                                                gttgttCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAgc                                                                                                                                                                                                                                     >  1:56799/1‑212 (MQ=255)                                                gttgttCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAgc                                                                                                                                                                                                                                     <  2:56799/212‑1 (MQ=255)                                                        ggTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGGGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTg                                                                                                                                                                                      >  2:249511/1‑251 (MQ=255)                                                                              aCGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACcgagcc                                                                                                                                                                >  1:24459/1‑247 (MQ=255)                                                                                  aTACCCGGTGTGATGTTCCAGTTGCCGATGTCGATTTTGTCATCCCGATACCAGGCGTGCGCCTCGGTGCCGGCACGCGTATCGCGGTCGTAAGGGCTTGCGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGCTTCTTTCACCTAGCGATAGCCCACGCCCACTTCGTGCGCTGAAGGGCCGATCATAAAGATCTGTCTGTAGCCTGGCTCAATACCGCGCAcccc                                                                                                                                                            >  2:223479/1‑250 (MQ=255)                                                                                        ggCGTGATGGTCCAGTTGCCGATGTCGATTTTGGCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGtt                                                                                                                                                      <  2:122882/251‑1 (MQ=255)                                                                                              aTGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGCGCGCGGAAGGGCCGATCATAAAGCTCTGGCTGTAGCCTGGCCCAATACCGCGCACCCATCAGTTACgcgg                                                                                                                                                >  2:2651/1‑251 (MQ=255)                                                                                                   ccAGGTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGGGGGCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCAGTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGAGCATAAAGATCTGGCTGTAGCGTGACTCAATACCGCGCTCCCAGTAGTTACGAGGCgaga                                                                                                                                           <  2:44796/251‑1 (MQ=255)                                                                                                     aGTTGACGATGTCGATTTTGTAATCCAGATACCAGGCGTGCGCATCGGGGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGGGCCGGACGGCAACTGCCAGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTAACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGCTCTGGCTGAAGCGTGGCACAATACCGCGAACCCAGTAGTTACTCGGCGAGAgg                                                                                                                                         >  1:62742/1‑251 (MQ=255)                                                                                                               gTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGtt                                                                                                                               >  1:169271/1‑251 (MQ=255)                                                                                                               gTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGGTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCCATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTTATGAGtt                                                                                                                               >  2:193543/1‑251 (MQ=255)                                                                                                                   aTTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCCTTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGTTGCGTTTGcc                                                                                                                           >  2:112648/1‑251 (MQ=255)                                                                                                                   aTTTTGTAATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGACGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCCCGCCCACTTCGTGCGCGGAAGGGACGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGcc                                                                                                                           >  2:240705/1‑251 (MQ=255)                                                                                                                    ttttGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGGGAGGGTGATGCGTTTGCCt                                                                                                                          >  2:23026/1‑251 (MQ=255)                                                                                                                      ttGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCCCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCtt                                                                                                                         <  1:88715/250‑1 (MQ=255)                                                                                                                       tGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTg                                                                                                                        <  2:96186/250‑1 (MQ=255)                                                                                                                                     cAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCGTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGTGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCt                                                                                                         <  2:170535/251‑1 (MQ=255)                                                                                                                                      aGGCGTGCGCCTCGGTGCCGGAAAGCGTCTCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACCTAGCGATAGCCCCCGCCCACTTCGTGCGCGGACGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATAACGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGAGTTTGCCTTTCTCCAGGTAGCCGCTg                                                                                                        >  2:22589/1‑251 (MQ=255)                                                                                                                                            gcgcCTCGGGGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTCACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAgg                                                                                                  <  2:86347/251‑1 (MQ=255)                                                                                                                                                cTCGGTGCCGGAACGCGTATCGCGGTCGTAAGTGCTAGAGCCGGACGGCAACTGCCCGCTGCTGGAGGCGGAGTAGTAACGCATTTCACGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGTCGATCCTAAAGTTCTGGCTGTAGCGTGGCACAATACCGCGCACCACGTCGGTACGCGGCGCGGGGGTGATGCGGTTGCCTTGCTCCAGGAATCCGCTgcgacgggctt                                                                                              >  2:145103/1‑243 (MQ=255)                                                                                                                                                  cGGTGCCGGAACGCGTATCGCGGGCGTAAGGGCTTGATCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGAGAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCCCCCAGTAGTTACGCTGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGttt                                                                                              <  1:174661/249‑1 (MQ=255)                                                                                                                                                   ggTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTggg                                                                                           >  2:59441/1‑251 (MQ=255)                                                                                                                                                        cGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCCGGTAGCCGCTg                                                                                                        >  2:186534/1‑233 (MQ=255)                                                                                                                                                        cGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTg                                                                                                        <  1:186534/233‑1 (MQ=255)                                                                                                                                                                gTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGGAGCCGCTGCGCAGGGGTTGGGTGTAGCACCCCTg                                                                              >  2:157570/1‑251 (MQ=255)                                                                                                                                                                                    tGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATgctggc                                                          >  1:214764/1‑251 (MQ=255)                                                                                                                                                                                     gAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATgctggct                                                         >  1:203515/1‑251 (MQ=255)                                                                                                                                                                                     gAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATgctggc                                                          <  1:249511/250‑1 (MQ=255)                                                                                                                                                                                          ggACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTctg                                                    >  1:216267/1‑251 (MQ=255)                                                                                                                                                                                                            gctgGTGGCGGTGTAGTAACGAATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGTCTCAATACCGCGCAACCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTGATGCTGGCTGTCTGGCTGGAACTGGTAGCCCa                                  <  1:223479/251‑1 (MQ=255)                                                                                                                                                                                                               ggtggCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAgg                                <  2:169271/250‑1 (MQ=255)                                                                                                                                                                                                                      gtgtAGTAACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCAATCATAAAGATCTGGCTGTAGCGAGGCTCAATACCGCGCACCCAGTAGTTACGAGGCGAGAGGGCGATGCGTTTGCCTTGATCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGcc                          >  1:34596/1‑249 (MQ=255)                                                                                                                                                                                                                          agtaACGCATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCtt                     <  1:47231/250‑1 (MQ=255)                                                                                                                                                                                                                                 cATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCgacgac             >  1:44426/1‑251 (MQ=255)                                                                                                                                                                                                                                          cGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGc    <  1:60447/251‑1 (MQ=255)                                                                                                                                                                                                                                            ttGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGGAGCgg  >  1:7002/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               GGAAGCGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATGACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGG  >  NC_000913/4515101‑4515586

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