Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,517,355	C→T	Q121Q (CAG→CAA)	fecR ←	ferric citrate regulator FecR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,517,355	0	C	T	100.0%	145.0 / NA	43	Q121Q (CAG→CAA)	fecR	ferric citrate regulator FecR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/22);  total (21/22)

CACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCA  >  NC_000913/4517112‑4517598                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    cacaGCGTGCTGCTGCAGGTCAAGCTGCGAGAAATTATCCTGCCGGCGGAAGGTAAATTCTGTCCCTAAAGCGGGGAGCTGTACATGACGGGTCAAGACACGGAAGGGGCGTTGCAGGGCATCTTTCGTGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACCGTGCGCTGCTGCGCATCAAAACGCAAACCCGCCGCGCTTTGGGTATTCAGCGCGAGCAGGGAGCAATCTTCCAGTTTTTGACGGc                                                                                                                                                                                                                                              >  2:203666/1‑251 (MQ=255)         gctgctgcACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACtgtg                                                                                                                                                                                                                                      >  1:121549/1‑251 (MQ=255)            gctgcCCGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCt                                                                                                                                                                                                                                   >  2:197489/1‑251 (MQ=255)                    tCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggt                                                                                                                                                                                                                           >  1:274046/1‑251 (MQ=255)                     cAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACACACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtg                                                                                                                                                                                                                          >  1:172650/1‑251 (MQ=255)                       aGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgcg                                                                                                                                                                                                                        <  1:173681/251‑1 (MQ=255)                       aGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtg                                                                                                                                                                                                                          >  2:297508/1‑249 (MQ=255)                         cTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGgcggtgcggg                                                                                                                                                                                                                      >  1:257212/1‑250 (MQ=255)                             gTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGGAATc                                                                                                                                                                                                                  >  2:275724/1‑251 (MQ=255)                               gAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTGGGCGGTGCGGTAATCt                                                                                                                                                                                                                 <  1:256550/250‑1 (MQ=255)                                       ccTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAg                                                                                                                                                                                                         <  1:224076/250‑1 (MQ=255)                                               cGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACACACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGcc                                                                                                                                                                                                <  2:172650/251‑1 (MQ=255)                                                   cGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGACCGCAGACCTTCGCCGGtt                                                                                                                                                                                            >  1:204720/1‑251 (MQ=255)                                                      tAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTcc                                                                                                                                                                                         >  2:275394/1‑251 (MQ=255)                                                                  cTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGc                                                                                                                                                                             >  2:108670/1‑251 (MQ=255)                                                                                tGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGcc                                                                                                                                                               >  1:190637/1‑251 (MQ=255)                                                                                 ggCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGcc                                                                                                                                                               <  1:123588/250‑1 (MQ=255)                                                                                     cTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGAAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAATCGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAgcgc                                                                                                                                                          >  1:141606/1‑251 (MQ=255)                                                                                         cGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAg                                                                                                                                                      <  1:64741/251‑1 (MQ=255)                                                                                           ggTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCCCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCa                                                                                                                                                    <  1:82319/251‑1 (MQ=255)                                                                                               aGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCACTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGCCTGCCAGAGCTGCCAGCCTCCGCCCGCGCCGAGCAAc                                                                                                                                                  >  2:147553/1‑249 (MQ=255)                                                                                                      ggAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGGTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCCGTaa                                                                                                                                          >  2:232097/1‑250 (MQ=255)                                                                                                         aGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCCCCGCCAGCGCCTAGCAACAGCAGTAATCCt                                                                                                                                      >  2:221996/1‑251 (MQ=255)                                                                                                           gggggtggCAGGGCATCTTTCGCGGTAGTAATGGCGATTTACCCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCttt                                                                                                                                    <  2:245235/244‑1 (MQ=255)                                                                                                                  gCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCAc                                                                                                                              <  1:139710/250‑1 (MQ=255)                                                                                                                  gCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTGCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCAc                                                                                                                              <  1:139703/250‑1 (MQ=255)                                                                                                                     gggCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGa                                                                                                                          <  1:42619/251‑1 (MQ=255)                                                                                                                       gCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGAcg                                                                                                                        >  1:99432/1‑251 (MQ=255)                                                                                                                            tttCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCggg                                                                                                                   >  1:112987/1‑251 (MQ=255)                                                                                                                             ttCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAACCGCCCATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCCTCTTCCAGTTGTTGACGGCTGCCTGTGCCTTTGGCGGAGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCCTCACGTCCCGGCgggg                                                                                                                  >  2:41955/1‑250 (MQ=255)                                                                                                                                 cgGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTCATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCggggggcg                                                                                                              >  2:24553/1‑246 (MQ=255)                                                                                                                                         aTGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGGGGGGCggggatg                                                                                                      >  2:70929/1‑247 (MQ=255)                                                                                                                                                             gagCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCaa                                                                                  <  1:127430/251‑1 (MQ=255)                                                                                                                                                              ggCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTGGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAAc                                                                                 <  2:274046/250‑1 (MQ=255)                                                                                                                                                                  ggACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCg                                                                             <  1:147553/251‑1 (MQ=255)                                                                                                                                                                     ggttGCGCTGATTCTCATCAAACCCGACATCCGCCGCGCTTTGGGTATTCTGCGTGGGTAGGGGGCCATCTTCCAGTTGTTGACGGCTGACGGTGCCGTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACAGCCATCGCTGCCAGTCTCCGCCGGCGCCGAGCAACAGCAGTAACCCTTTCATGACGTGATGGCGGGTGAGTCGGGTATCGTGCCACGCCCGGTTCGCAACGTCTCCa                                                                          <  2:263009/248‑1 (MQ=255)                                                                                                                                                                       gTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAgg                                                                        <  1:211130/251‑1 (MQ=255)                                                                                                                                                                                        aaaaCGCACATCCGCCGCGCTTTGGGTATTCAGCGGGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTTACTGTGCCTTTGGCGGTGCGGTAAGCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCGGCCAGCCTCCGCCAGCTCCGAGCAACAGCAGTAATCCTTTCATCGCGGGTCGGCGGGTGAGGCGGGTAGCGTGCAGCGCCCGGTGGGCAACGTCGCCAGGCACACCGCCAAGCTGGt                                                       <  1:178928/251‑1 (MQ=255)                                                                                                                                                                                                  tCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGGGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGTACTCCGCCAAGCTGGTTGTGCAgg                                               <  2:294599/249‑1 (MQ=255)                                                                                                                                                                                                    cgccgcGCTTGGGGTATTCAGCGTGAGCAGGGAGCCCTCTTCCAGTTGTTGACGGCGGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTGCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGtttt                                           <  1:181147/251‑1 (MQ=255)                                                                                                                                                                                                       cgcgCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCaa                                        <  1:76861/251‑1 (MQ=255)                                                                                                                                                                                                                       aGCGTGAGCAGGGAGGCATCTTCCAGTTGTTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGCCCTTCGCCGGTTTCCGACTGCCAGGGCTGCCAGCCTCCGCCAGCGCCGCGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCa                        <  2:80981/251‑1 (MQ=255)                                                                                                                                                                                                                                             ccAGTTGTTGACGGTGGACGGTGCCTTTGGCGGTTGGGTAATCTGCCCGCAGACCTTCGGCGGGTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCAGCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCGGCTGCCAGGCCCACTGGTTATACTGATCCTGTTCa  <  2:184259/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    CACAGCGTGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGCCGGTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCA  >  NC_000913/4517112‑4517598

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