Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,514,053	T→G	T86P (ACG→CCG)	fecB ←	ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,514,053	0	T	G	100.0%	213.2 / NA	59	T86P (ACG→CCG)	fecB	ferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (34/25);  total (34/25)

TCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGC  >  NC_000913/4513824‑4514273                                                                                                                                                                                                                                        |                                                                                                                                                                                                                             tctTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACgg                                                                                                                                                                                                                      >  1:747041/1‑238 (MQ=255)                 gACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACCGCCTGCCACGGTTTCAgg                                                                                                                                                                                                     >  1:260624/1‑239 (MQ=255)                      gCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGACCAGGTCGGGTTTCAGAGCGGCAATGGCCTCCAGGCCCGGCTGCGCGCGCGGTCCGACGGCCTGCCACGGTTTCCGgcgcgc                                                                                                                                                                                                >  2:704165/1‑239 (MQ=255)                         tGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGcgcacg                                                                                                                                                                                             <  1:366783/239‑1 (MQ=255)                           cATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCACGCGCGGTCCGAAGGACTGCCACGGTTTCAGGTGcgcacgc                                                                                                                                                                                            >  1:1165757/1‑238 (MQ=255)                                      ttttACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAgg                                                                                                                                                                                <  2:747041/239‑1 (MQ=255)                                             caccaTTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGtt                                                                                                                                                                         >  2:1104619/1‑239 (MQ=255)                                               ccaTTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGtttt                                                                                                                                                                       >  2:227984/1‑239 (MQ=255)                                                caTTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTg                                                                                                                                                                      <  2:456889/239‑1 (MQ=255)                                                          atgatAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGcgcacg                                                                                                                                                                                             <  1:359679/206‑1 (MQ=255)                                                          atgatAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGcgcacg                                                                                                                                                                                             >  2:359679/1‑206 (MQ=255)                                                          atgatAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtc                                                                                                                                                            >  1:678218/1‑239 (MQ=255)                                                           tgatAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtc                                                                                                                                                            <  2:678218/238‑1 (MQ=255)                                                             atAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGCTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtcgcc                                                                                                                                                         >  1:1180297/1‑237 (MQ=255)                                                                                   cAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGcgcacg                                                                                                                                                                                             <  2:607729/181‑1 (MQ=255)                                                                                   cAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGcgcacg                                                                                                                                                                                             >  1:607729/1‑181 (MQ=255)                                                                                   cAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATAACGATCGGGATGACg                                                                                                                                   >  2:835371/1‑239 (MQ=255)                                                                                        tAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCAc                                                                                                                              >  1:696336/1‑239 (MQ=255)                                                                                                                  cagTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTg                                                                                                                                                                      >  1:539773/1‑173 (MQ=255)                                                                                                                  cagTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTg                                                                                                                                                                      <  2:539773/173‑1 (MQ=255)                                                                                                                   agTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGcgcacg                                                                                                                                                                                             <  1:762374/149‑1 (MQ=255)                                                                                                                   agTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGcgcacg                                                                                                                                                                                             >  2:762374/1‑149 (MQ=255)                                                                                                                    gTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGt                                                                                                  >  1:67262/1‑239 (MQ=255)                                                                                                                      aCCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGt                                                                                                  <  2:905942/239‑1 (MQ=255)                                                                                                                       ccGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCAt                                                                                                                                                                   >  2:972264/1‑171 (MQ=255)                                                                                                                       ccGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTcc                                                                                               >  1:985234/1‑239 (MQ=255)                                                                                                                       ccGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACGGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCAt                                                                                                                                                                   <  1:972264/171‑1 (MQ=255)                                                                                                                          gcgcGATTTGCTGCAAGGGGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGcc                                                                                                                                                                                                                <  1:1011596/123‑1 (MQ=255)                                                                                                                          gcgcGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGcc                                                                                                                                                                                                                >  2:1011596/1‑123 (MQ=255)                                                                                                                          gcgcGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTc                                                                                                <  2:206660/235‑1 (MQ=255)                                                                                                                          gcgcGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTc                                                                                                >  1:206660/1‑235 (MQ=255)                                                                                                                                        aaGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtc                                                                                                                                                            >  1:396900/1‑161 (MQ=255)                                                                                                                                        aaGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTAtc                                                                                                                                                            <  2:396900/161‑1 (MQ=255)                                                                                                                                               tGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGCTGTGGCGtt                                                                       >  1:971263/1‑239 (MQ=255)                                                                                                                                                           cGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCa                                                                                                                                                                    >  1:1120235/1‑136 (MQ=255)                                                                                                                                                           cGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCa                                                                                                                                                                    <  2:1120235/136‑1 (MQ=255)                                                                                                                                                           cGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGGGGCGTTTTTTCg                                                                 >  1:72764/1‑235 (MQ=255)                                                                                                                                                                        ctgTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCgtg                                                >  1:993807/1‑239 (MQ=255)                                                                                                                                                                         tgTCGGCAATGGTCAGGTCGGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCgtgt                                               <  2:1159673/239‑1 (MQ=255)                                                                                                                                                                         tgTCGGCAATGATCAGGTCTGGTTTCAGGGCGGCTATGGCTGCGAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGGCACGGTTTCAGGTGCGCACGCACTTCGGGGCGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCTAGCACCACAATTCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCTgtgt                                               <  2:1180297/239‑1 (MQ=255)                                                                                                                                                                            cGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCg                                            >  2:219484/1‑239 (MQ=255)                                                                                                                                                                            cGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGGCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCCCCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTc                                             <  2:67262/238‑1 (MQ=255)                                                                                                                                                                                  tGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAgg                                                                                                                                                                                                     <  1:1034368/80‑1 (MQ=255)                                                                                                                                                                                  tGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAgg                                                                                                                                                                                                     >  2:1034368/1‑80 (MQ=255)                                                                                                                                                                                   gATCAGGGCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGaa                                     <  1:835371/239‑1 (MQ=255)                                                                                                                                                                                            tGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTggcg                             >  1:1127447/1‑238 (MQ=255)                                                                                                                                                                                              gTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATAc                                                                                                                                                 <  2:301612/120‑1 (MQ=255)                                                                                                                                                                                              gTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATAc                                                                                                                                                 >  1:301612/1‑120 (MQ=255)                                                                                                                                                                                              gTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCgtgt                                               <  1:783854/218‑1 (MQ=255)                                                                                                                                                                                              gTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCgtgt                                               >  2:783854/1‑218 (MQ=255)                                                                                                                                                                                                   agagCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGc                     >  2:645606/1‑239 (MQ=255)                                                                                                                                                                                                   agagCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGGCCTGAACCGGGGCGGCCAAGGc                     >  2:85821/1‑239 (MQ=255)                                                                                                                                                                                                      gCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGt                   <  2:993807/238‑1 (MQ=255)                                                                                                                                                                                                          cAATGGCTTCCAGGCTCGGCTGCGCGCGCGGTCCGACg                                                                                                                                                                                                                       >  1:121303/1‑38 (MQ=255)                                                                                                                                                                                                                cTTCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCAc         <  2:696336/238‑1 (MQ=255)                                                                                                                                                                                                                 ttCCAGGCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCAcc        <  1:645606/238‑1 (MQ=255)                                                                                                                                                                                                                      ggCTCGGCTGCGCGCGCGGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACcagcagc  >  2:442733/1‑239 (MQ=255)                                                                                                                                                                                                                                 gcgcgcgGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGAt                                                                                                                                       <  2:351545/95‑1 (MQ=255)                                                                                                                                                                                                                                 gcgcgcgGTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGAt                                                                                                                                       >  1:351545/1‑95 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                             TCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAA‑‑CCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGC  >  NC_000913/4513824‑4514273

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