Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,183,258	G→A	E672K (GAA→AAA)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,183,258	0	G	A	100.0%	217.3 / NA	62	E672K (GAA→AAA)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/31);  total (31/31)

TCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCA  >  NC_000913/4182970‑4183544                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                cccGTTTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACg                                                                                                                                                                                                                                                                                     <  1:412427/300‑1 (MQ=255)  ccGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGc                                                                                                                                                                                                                                                                                    <  2:218340/301‑1 (MQ=255)           cacaGACTAACGAATACGGCGTCCTTGAGACTCCGTATCGGAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGc                                                                                                                                                                                                                                                                           <  1:543488/301‑1 (MQ=255)                  aaaCGACAACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATg                                                                                                                                                                                                                                                                    <  1:576491/300‑1 (MQ=255)                             gctgcatggagaCTCCGTATCTTAAAGGGAGCGACGGGGGTAGAACGGACGAAAGTCAAGACCGGGCTGCAATCGCAGAAGGAAACTACGTTATCGCCAAGGGGAAATCCAACGTGTGTGAAGAAGGCCACTGCGTGGTAGACCGGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGGTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAAc                                                                                                                                                                                                                                                          <  1:217470/291‑1 (MQ=255)                                                 cGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCCGGCCGTTCCGa                                                                                                                                                                                                                                    >  1:319964/1‑301 (MQ=255)                                                        tGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTgc                                                                                                                                                                                                                              <  2:175530/300‑1 (MQ=255)                                                               cGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCGGCGCGCTGAt                                                                                                                                                                                                                      >  1:462247/1‑301 (MQ=255)                                                                     tgtAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCg                                                                                                                                                                                                                >  2:74557/1‑301 (MQ=255)                                                                          cTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCAGTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCgggg                                                                                                                                                                                                           >  2:449605/1‑300 (MQ=255)                                                                           tGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt                                                                                                                                                                                                          >  2:46693/1‑301 (MQ=255)                                                                             aCGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttgg                                                                                                                                                                                                        >  1:66041/1‑301 (MQ=255)                                                                              cGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTGCGTAGAAGACCTGGTAACTTGCCGTAGCACAGGCGAATCCAGCATGTTCAGCCGCGACAAGCTTGACTACTTGGACGTATCCACCCAGCAGGTTGTAACCGTCGGTGCGTCCGTGACCCCGTTCCTGAAACACGATGACCCCAACCGCGCATTG‑TGTGTCTCGAACATGCAACGTCACACAGTTCCGACACTgcccgcattcaacactccggtggtt                                                                                                                                                                                                       >  1:14464/1‑282 (MQ=255)                                                                               gAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa                                                                                                                                                                                                      >  2:127472/1‑301 (MQ=255)                                                                                 caTTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGATGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCTAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCATCAGGTGGTCTCCGTCGATGCGTCCCTGATCCCGTTACGGAAACACGAAGACCACAATCGTGCATGGATGGGT‑GTGACCATGCCAAGTCAGGACGTTCCGACGCGAGGCGATGTTAAGCTGCTGGTTTGTaaa                                                                                                                                                                                                    >  1:569462/2‑299 (MQ=255)                                                                                 aaTTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATCGACGTATCCACCCAGCAGGTCGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACCCTGCGCGCTGCTAAGCCGCTGGTTGGTa                                                                                                                                                                                                      >  1:256290/1‑299 (MQ=255)                                                                                     cACTAACTGTCTGCAATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACATGGATGAAGAAGGACACTTCGTAGAAGACCTGGTAACTGGACGTAGCAAAGGCGAATCCTGCTTGTACAGCCGCGACCAGGTTGACTACATGGACCTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGG‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTa                                                                                                                                                                                                <  1:46650/301‑1 (MQ=255)                                                                                            tGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACg                                                                                                                                                                                         <  2:66127/301‑1 (MQ=255)                                                                                                           aagGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTg                                                                                                                                                                           <  2:462818/300‑1 (MQ=255)                                                                                                                  tctcGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTAATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGgt                                                                                                                                                                   <  2:159747/298‑1 (MQ=255)                                                                                                                     cGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTAATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGgtaga                                                                                                                                                                >  1:159550/1‑298 (MQ=255)                                                                                                                       ttatcactccgtcaaaCTTCGACTGGGATGAAGACGTCCCCTTCGTTGAAGCCTTGTTTACTTGCCTGAACGAAGGCGATTCCAGCTTGTTCTCTCTCGCCCCACTTGACTCCATGGATGTCTGCACCCATCATGTGGTATTCGGTGATGGGTCCATGATCCCTCTCCTGAATCGCGCTGACGTCAAACGGGCGTTGAGGGGT‑GCGAACATGCAACGTCAGGCCGTTCGGACCCTGCTCGGTGATGAGGCGCTGGTTGGCACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAAc                                                                                                                                                              <  1:155002/287‑1 (MQ=255)                                                                                                                           cGCCCCGGCGAAATCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACTTATCCCCCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCg                                                                                                                                                          <  1:42378/301‑1 (MQ=255)                                                                                                                               cAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGGTGCACCCGGCGGAATTGCGGTAg                                                                                                                                                      >  1:421556/1‑301 (MQ=255)                                                                                                                                 ggCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAg                                                                                                                                                      >  2:545557/1‑299 (MQ=255)                                                                                                                                            aCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGTCGAATACAGCTTGTTCAGCCACGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTACCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTTCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACgtggtgg                                                                                                                                         <  2:402747/301‑1 (MQ=255)                                                                                                                                                 gATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCg                                                                                                                                    >  2:350471/1‑301 (MQ=255)                                                                                                                                                 ccgaggctcttccgatcggTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCg                                                                                                                                    <  1:527594/283‑1 (MQ=255)                                                                                                                                                        aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAg                                                                                                                               <  2:325181/299‑1 (MQ=255)                                                                                                                                                          ggCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACg                                                                                                                           >  2:84502/1‑301 (MQ=255)                                                                                                                                                              aCTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTggg                                                                                                                       >  2:169225/1‑300 (MQ=255)                                                                                                                                                                  cGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATg                                                                                                                     >  2:344229/1‑299 (MQ=255)                                                                                                                                                                   gTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCgagatcggaagagcacacg                                                                                                                  >  2:528259/1‑283 (MQ=255)                                                                                                                                                                   gTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCtt                                                                                                                  >  2:146993/1‑301 (MQ=255)                                                                                                                                                                            cTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCg                                                                                                         >  2:363125/1‑301 (MQ=255)                                                                                                                                                                                tAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTa                                                                                                      <  2:36366/300‑1 (MQ=255)                                                                                                                                                                                   cTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCaa                                                                                                  >  2:22813/1‑301 (MQ=255)                                                                                                                                                                                      gCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGt                                                                                               >  2:539922/1‑301 (MQ=255)                                                                                                                                                                                              aaaGGCAAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAg                                                                                       <  1:127309/301‑1 (MQ=255)                                                                                                                                                                                                    gAATCCAGCTGGGTTAAGCGGGACCATGTTGACTACATGGACGTATCCACCGAGCAGGTGGTACCCGTCGGTGCGTCCCTGGTCCCGTTCCTGAAACACGATGACGCCAACAGGTCATTGATGGGT‑GCGAACATGCAACGTGAGGCCGTTCCGACTCTGCGCGCTGAAAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGTCGTTGACTCCGGTGTAACTGCGGTAGCTACACGTGGTGGTGTCGTTCAGTACGTGGATGCTGCCCGTATCGTTATCAAAGTTAACGAAGACgag                                                                                  <  2:570173/300‑1 (MQ=255)                                                                                                                                                                                                     aaTCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAATTGCGGTAGACAAACGTGGTGGTTCCGTTCAGTACGTGGGTGCTTCCCGTTTCGGTAtatacttcaatgaagactgaga                                                                                >  1:484112/1‑280 (MQ=255)                                                                                                                                                                                                     aaTCCAGCTTCTTCAGCCGCGACCAGGTTGACTACATGGACGTAGCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGAt                                                                                <  1:22790/301‑1 (MQ=255)                                                                                                                                                                                                          aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTAt                                                                            >  2:323748/1‑300 (MQ=255)                                                                                                                                                                                                                ccccgacacgacgctcttccgatctttgACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGCTACCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTg                                                                     <  1:125650/274‑1 (MQ=255)                                                                                                                                                                                                                      cgcgACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAg                                                               <  1:146808/301‑1 (MQ=255)                                                                                                                                                                                                                      acgtgtgctcttccgatcttGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCa                                                                <  2:65207/281‑1 (MQ=255)                                                                                                                                                                                                                        cgACCAGGTTGACTAAATGGACGCATCCACCCAGCAGGTGGAACCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGAGGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGAGTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAgg                                                              <  1:350071/300‑1 (MQ=255)                                                                                                                                                                                                                             aGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGa                                                        >  2:274477/1‑301 (MQ=255)                                                                                                                                                                                                                                  gACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGTAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACAt                                                     >  2:289468/1‑299 (MQ=255)                                                                                                                                                                                                                                         tGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCACAGTTAACGAAGACGAGATGTATCCGGGTGACGCaagatcggaagggcggcgtgg                                            >  1:65121/1‑281 (MQ=255)                                                                                                                                                                                                                                         tGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGagatcggaagagcacacgtctgaaa                                            >  2:125813/1‑277 (MQ=255)                                                                                                                                                                                                                                          ccgacgctctgccgatctggtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTg                                           <  1:118901/283‑1 (MQ=255)                                                                                                                                                                                                                                           taCGTAGCCACGCAGCTGGGGGTATCCGTCGGTGCGTCCCTGATCCCGTGCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTg                                           <  1:74471/299‑1 (MQ=255)                                                                                                                                                                                                                                           gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGa                                          >  2:97228/1‑301 (MQ=255)                                                                                                                                                                                                                                                    ccccAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATaca                                 <  1:449058/300‑1 (MQ=255)                                                                                                                                                                                                                                                            ggtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGagatcggaagagcacacg                         >  2:119041/1‑284 (MQ=255)                                                                                                                                                                                                                                                             gtggtATCCGGCGGTGCGTCCCTGATCCCGGTCCGGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGtt                           <  1:169028/298‑1 (MQ=255)                                                                                                                                                                                                                                                               gggatacGTCGGGGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTaa                       <  1:84395/294‑1 (MQ=255)                                                                                                                                                                                                                                                                gtATCCGTCGTTGCGTCCCTGAGCCCGTTCCTGAAACACGATAACTCAAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGTCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAAc                      <  1:289137/300‑1 (MQ=255)                                                                                                                                                                                                                                                                  aTCCGTCGGTGCGTCCCTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAg                   <  2:422051/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                cccTGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTCTCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACATGACCAAATACACCCGTTCTAACCAGAACACCTGTATCaa     >  2:443065/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                   tGATCCCGTTCCTGAAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCa  <  2:320321/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                TCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAAGGCAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGT‑GCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAGCTAAACGTGGTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCA  >  NC_000913/4182970‑4183544

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