Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,174,041	A→T	E320D (GAA→GAT)	birA →	DNA‑binding transcriptional repressor/biotin‑[acetyl‑CoA‑carboxylase] ligase BirA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,174,041	0	A	T	100.0%	287.2 / NA	73	E320D (GAA→GAT)	birA	DNA‑binding transcriptional repressor/biotin‑[acetyl‑CoA‑carboxylase] ligase BirA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (34/39);  total (34/39)

CGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCT  >  NC_000913/4173797‑4174280                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                cggcCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCGTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGTGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGGATTGCAGataaataa                                                                                                                                                                                                                                           <  2:395875/251‑1 (MQ=255)     cATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGaaa                                                                                                                                                                                                                                       <  1:468108/251‑1 (MQ=255)           aaTACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAgggg                                                                                                                                                                                                                                   >  2:674043/1‑249 (MQ=255)            aTACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGCTTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGt                                                                                                                                                                                                                                >  2:56845/1‑251 (MQ=255)               cGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCAGTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTTAGCAGGGTGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTAtt                                                                                                                                                                                                                             <  2:548824/251‑1 (MQ=255)                  gAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGc                                                                                                                                                                                                                          <  1:534085/251‑1 (MQ=255)                    aTTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTc                                                                                                                                                                                                                        <  1:219713/251‑1 (MQ=255)                         gTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAA‑CAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGc                                                                                                                                                                                                                  >  1:20903/1‑251 (MQ=255)                          tgctgcGTTGGAACTCTTCGAACAAGAAGGATTGGCCCCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGCTAAAGAAATCTTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGc                                                                                                                                                                                                                  >  2:818596/1‑251 (MQ=255)                                  tGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAattat                                                                                                                                                                                                           <  1:674043/250‑1 (MQ=255)                                    gAACTCTTCGAACAAGAAGGATTTGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCAGTTCACGCGGAATAGACAAACAGGGGGCTTTATGACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATAGCCCTGCGGAGGGCAGATAAATAAGAAAGGGGAGTAGTCGCTCCCCTGCAAATTATTTg                                                                                                                                                                                                        <  2:853367/251‑1 (MQ=255)                                            cGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATGAGTTGCGTAGTCt                                                                                                                                                                                                <  1:314796/251‑1 (MQ=255)                                             gAACAAGAAGGCTTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATACTTTTATTAATCGCCCACTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGAGGTGAAATATCCCTGCGTAGAGCAGATAAATAAGAAAGGGCAGAATTAGCTCCCCTGCAAATTATTTGCGTAGTCTg                                                                                                                                                                                               >  1:164267/1‑251 (MQ=255)                                                     agGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCtcttc                                                                                                                                                                                       >  1:236153/1‑251 (MQ=255)                                                     agGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCtcttc                                                                                                                                                                                       <  2:477486/251‑1 (MQ=255)                                                     agGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCGGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCtcttc                                                                                                                                                                                       <  1:920246/251‑1 (MQ=255)                                                           ggCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGc                                                                                                                                                                                 <  1:332964/251‑1 (MQ=255)                                                           ggCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGAGAAACTTCTAATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATA‑GAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGcc                                                                                                                                                                                >  1:854188/1‑250 (MQ=255)                                                             cACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCAt                                                                                                                                                                               >  2:842900/1‑251 (MQ=255)                                                              aCCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATg                                                                                                                                                                              <  2:769772/251‑1 (MQ=255)                                                                              tGGGAAAAGCTGGATCATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTc                                                                                                                                                              >  2:94508/1‑251 (MQ=255)                                                                              tGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAACCTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTc                                                                                                                                                                 >  2:97590/1‑248 (MQ=255)                                                                              tGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTc                                                                                                                                                                 <  1:97590/248‑1 (MQ=255)                                                                                  aaaaGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGCGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAgg                                                                                                                                                           >  2:589659/1‑250 (MQ=255)                                                                                   aaaGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGAt                                                                                                                                                         <  1:810914/251‑1 (MQ=255)                                                                                             aagtttATTAAGCGCCCAGTGAAACTTATCAGTGGTGATAAAGAAAGATTTGGCATTTAACGCGGAATAGACAAACAGGGGGCTTTATTTCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCATGATTAAACTGgc                                                                                                                                                <  2:542159/247‑1 (MQ=255)                                                                                             aaTTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAA‑GAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGGCCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGgcgc                                                                                                                                              <  1:95132/251‑1 (MQ=255)                                                                                                   attaATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAA‑CAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCa                                                                                                                                         <  2:20903/250‑1 (MQ=255)                                                                                                       aTCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACg                                                                                                                                       >  1:826274/1‑249 (MQ=255)                                                                                                         cGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGGCAGTATTAAACTGGCGCGCTCACGAgt                                                                                                                                    <  1:714258/250‑1 (MQ=255)                                                                                                           cccAGTGAAAATTATCATCGGTGATAATGATATATTTGGCATTTCAAGCGGTATAGCCTAACAGGGGGCTTTATTAATTGAGCAGGATGGAATAATAAAAACCTGGCTGGGTGGTGAAATATCGCTGCGTCGTGCAGATAAATACGACAGGGGAGTATTCGCTACCCTGCCAATTATTTGCGTAGTCTGACCTCTTCTACAGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCACACGAgtagg                                                                                                                                 >  2:589335/1‑251 (MQ=255)                                                                                                                  aaaCTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATg                                                                                                                                                                              <  1:718116/199‑1 (MQ=255)                                                                                                                  aaaCTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATg                                                                                                                                                                              >  2:718116/1‑199 (MQ=255)                                                                                                                       tATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATAtttt                                                                                                                      <  1:22629/250‑1 (MQ=255)                                                                                                                        aTCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACACACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGc                                                                                                                    >  2:722422/1‑251 (MQ=255)                                                                                                                        aTCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGATTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCGCACGAGTAGGTAGAATATTTTg                                                                                                                     <  2:826274/250‑1 (MQ=255)                                                                                                                                      aGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGcc                                                                                                       >  1:546824/1‑250 (MQ=255)                                                                                                                                        aaaTATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTTACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGAGTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGt                                                                                                    <  2:546824/251‑1 (MQ=255)                                                                                                                                                ggCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGAtctc                                                                                             <  2:724297/250‑1 (MQ=255)                                                                                                                                                 gCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCtt                                                                                           >  2:318009/1‑251 (MQ=255)                                                                                                                                                       cACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCcac                                                                                      >  1:534400/1‑250 (MQ=255)                                                                                                                                                          gcgGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGCATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATg                                                                                  >  1:474947/1‑251 (MQ=255)                                                                                                                                                                tAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATgg                                                                            >  1:276147/1‑251 (MQ=255)                                                                                                                                                                    cAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGt                                                                        <  1:94508/251‑1 (MQ=255)                                                                                                                                                                    cAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGt                                                                        <  1:509589/251‑1 (MQ=255)                                                                                                                                                                          gggggCGTTATTTCTTGATCAGGATGGAATAAGAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTTCAAATTAGTTGCGTAGTCTGTCCTCTTCGACCGCATGGTGAGCGCTTTTCGGCAGGATTAAACCGACGCGCTCACGAGTAGGTAGAATATTTTGCTGTAAGGTCTGCGAGTGGATCCCTTTCCAGAAGGTCATGGCATTCTTAAt                                                                  <  2:929719/251‑1 (MQ=255)                                                                                                                                                                              gCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGct                                                              >  1:119175/1‑251 (MQ=255)                                                                                                                                                                                tttATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGcttc                                                            >  2:225533/1‑251 (MQ=255)                                                                                                                                                                                         tGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAACTATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGt                                                   >  1:901035/1‑251 (MQ=255)                                                                                                                                                                                         tGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTa                                                     <  1:301675/249‑1 (MQ=255)                                                                                                                                                                                                gATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAg                                                                                                              <  1:248424/185‑1 (MQ=255)                                                                                                                                                                                                gATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAg                                                                                                              >  2:248424/1‑185 (MQ=255)                                                                                                                                                                                                gATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTc                                            <  2:276147/251‑1 (MQ=255)                                                                                                                                                                                                 aTGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGCTTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTcg                                           >  2:82461/1‑251 (MQ=255)                                                                                                                                                                                                      ataataAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAg                                      <  2:119175/251‑1 (MQ=255)                                                                                                                                                                                                      ataataAAACCCTGGATGGGCGGAGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGt                                        >  2:102809/1‑249 (MQ=255)                                                                                                                                                                                                           aaaaCCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATg                                 >  2:563268/1‑251 (MQ=255)                                                                                                                                                                                                              aCCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATCAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGaa                               >  1:886319/1‑250 (MQ=255)                                                                                                                                                                                                               cccTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCAt                                                                              >  2:47141/1‑202 (MQ=255)                                                                                                                                                                                                               cccTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCAt                                                                              <  1:47141/202‑1 (MQ=255)                                                                                                                                                                                                                   ggATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTGAATTTCGCGTAGTTATGAAAATa                           <  1:818596/249‑1 (MQ=255)                                                                                                                                                                                                                       gggCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTGCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATc                     <  1:469078/251‑1 (MQ=255)                                                                                                                                                                                                                        ggCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATc                     <  2:46202/250‑1 (MQ=255)                                                                                                                                                                                                                        ggCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATc                     >  1:46202/1‑250 (MQ=255)                                                                                                                                                                                                                         gCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCg                   <  1:56845/251‑1 (MQ=255)                                                                                                                                                                                                                            gTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGt                <  1:842900/251‑1 (MQ=255)                                                                                                                                                                                                                                aaTATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTa                                   <  1:157311/228‑1 (MQ=255)                                                                                                                                                                                                                                aaTATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTa                                   >  2:157311/1‑228 (MQ=255)                                                                                                                                                                                                                                      ccTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATc                                                                 >  1:216334/1‑192 (MQ=255)                                                                                                                                                                                                                                      ccTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATc                                                                 <  2:216334/192‑1 (MQ=255)                                                                                                                                                                                                                                       cTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCcaca                                                                                     <  1:130215/171‑1 (MQ=255)                                                                                                                                                                                                                                       cTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCcaca                                                                                     >  2:130215/1‑171 (MQ=255)                                                                                                                                                                                                                                         gCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGcc     <  1:82461/249‑1 (MQ=255)                                                                                                                                                                                                                                          cGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCt  >  1:611871/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                CGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCT  >  NC_000913/4173797‑4174280

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