Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,341,582	C→T	R16H (CGC→CAC)	melR ←	DNA‑binding transcriptional dual regulator MelR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,341,582	0	C	T	100.0%	132.9 / NA	42	R16H (CGC→CAC)	melR	DNA‑binding transcriptional dual regulator MelR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/25);  total (17/25)

AAATCGCCATGCTCTGACAGGTTCCGGTATCTGTTAGTTGGTGCGGTGTACAGGCCCAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGATCTGAGTTT  >  NC_000913/4341343‑4341825                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    aaaTCGCCATGCTCTGACAGGTTCCGGTATCTGTTAGTTGGTGCGGTGTACAGGCCCAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGtttt                                                                                                                                                                                                                                          <  1:124834/251‑1 (MQ=255)  aaTCGCCAAGCTCTGACAGGTGCCGGTATCTGTTAGTTGGTGCGGTGTACAGGCCCAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGttttt                                                                                                                                                                                                                                         <  2:318629/251‑1 (MQ=255)              ctGACAGGTTCCGGTATCTGTTAGTTGGTGCGGTGTACAGGCCCAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCgctgct                                                                                                                                                                                                                               <  1:233669/249‑1 (MQ=255)                      ttCCGGTCTCTGTTAGTTGGTGCGGTGTACAGGCCCAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCCTCCGCTGGTATTCCGAATACCGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATacac                                                                                                                                                                                                                     >  2:43134/1‑248 (MQ=255)                       gccGGGATCTGTTAGTTGGTGCGGTGTACAGGCCCAGAACAGCGGGATATGACCCTGATTGATAGTCGCTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGGGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACg                                                                                                                                                                                                                    <  2:135972/250‑1 (MQ=255)                        ccGGTATCTGTTAGTTGGTGCGGTGTACAGGCCCAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCAGCGAAAGGCACATTCACGTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGTTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACg                                                                                                                                                                                                                    <  1:270560/250‑1 (MQ=255)                         cGGTATCTGTTAGTTGGTGCGGTGTACAGGCCCAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGt                                                                                                                                                                                                                   >  1:324630/1‑250 (MQ=255)                             aTCTGTTAGTTGGTGCGGTGTACAGGCCCAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCt                                                                                                                                                                                                               <  2:314042/250‑1 (MQ=255)                                    aGTTGGTGCGGTGTACAGGCCCAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTGTCGTCGCTGCTGCACATAAACGTATCTGTATTCa                                                                                                                                                                                                        <  1:323443/250‑1 (MQ=255)                                         gttaggTGAACAGGCCCAGAACAGCGTGATATGACCCTGAATGATTTTCACATTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACCTTCAATTAGACCTGACCATGCCAGTGGCTGGTAGGCATGATATTCGGGGCGAGAAACTCACTATCCATCCTCTGGCATTACTAATAAAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATGAACGAATCTGTCTTTatggctt                                                                                                                                                                                                  >  2:17509/5‑248 (MQ=255)                                                gtACAGGCACAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGGATCTGTATTCATGGATGGctctct                                                                                                                                                                                           >  2:213542/1‑251 (MQ=255)                                                               aGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGGTTTTCGTCGCTGCTGCACAGAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTTGAATATCAg                                                                                                                                                                            <  2:195312/251‑1 (MQ=255)                                                                                                    tgttgATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTgg                                                                                                                                       >  1:217841/1‑251 (MQ=255)                                                                                                         aTCAGGTATTCCACATCGCAATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGGGAGGGAGGATGACTGCGAGTgggggc                                                                                                                                   >  2:29385/1‑247 (MQ=255)                                                                                                              gTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGtt                                                                                                                             <  1:16478/251‑1 (MQ=255)                                                                                                                aTTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGttt                                                                                                                            <  1:240711/250‑1 (MQ=255)                                                                                                                aTTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCGCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGtttt                                                                                                                           <  1:44294/251‑1 (MQ=255)                                                                                                                             aTCGAAAGGCACATTCACTTCGACCTGACAATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCCCGGTTTTCACCCTCTTcccc                                                                                                              >  2:278482/1‑250 (MQ=255)                                                                                                                                gAAAGGCACATTAACTTCGACCTGACCATGCCAGTGGCTGGTGGGCCTGATATGCGGTGCGCGAAACTCAATCTACATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGGGTGGGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCag                                                                                                             >  2:40714/1‑249 (MQ=255)                                                                                                                                  aaGGCAAATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAggg                                                                                                         >  1:309862/1‑251 (MQ=255)                                                                                                                                                 cGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTcc                                                                                          <  2:324630/251‑1 (MQ=255)                                                                                                                                                 cGACCTGACCATGCAAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGATGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTACCAGAGGGGCGAGGGGGCTCTcc                                                                                          >  2:398848/1‑251 (MQ=255)                                                                                                                                                  gACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCg                                                                                         <  1:223867/251‑1 (MQ=255)                                                                                                                                                                       ggtggGCATGATATGCGTGGAGCGAATCTCAATCCCCATCCGATGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAAct                                                                    <  1:328076/251‑1 (MQ=255)                                                                                                                                                                                 atatGCGGTGCGCGAAACGCAAGCTCCAGCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCa                                                          <  2:417759/251‑1 (MQ=255)                                                                                                                                                                                   atGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAg                                                         <  1:900/250‑1 (MQ=255)                                                                                                                                                                                     gCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTa                                                      <  1:29385/251‑1 (MQ=255)                                                                                                                                                                                          gcgcGAAACTCAATCTCCTTCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTACCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAATCTCTGCTTTTCTGGTAATTTa                                                 <  2:377841/251‑1 (MQ=255)                                                                                                                                                                                              gAAACTCAATCTCCATCCGATGGTATTCCGACTACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGCGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCCCCGAGTATCCTGAGGCCGAAAACTCTGCTTTTACGGTAATTTATTcc                                             >  2:53901/1‑251 (MQ=255)                                                                                                                                                                                                  cTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCCCATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCCCCCTCTTCCCAGAGGGGCGCGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTCTTCCCAt                                          >  1:257448/1‑250 (MQ=255)                                                                                                                                                                                                          ccATCCGCTGGTATGCCGAATACAGCGACAGCGGGTTGTGGGTCTGTTTTTCGTCTCTGCTGCTCATAAACGTATCTTTATTCATGTATGGCTCTCTTTCCTGGAATATCAGAATTAGGGCAGGAGTGAGGGATGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCTAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCa                                   <  1:54510/249‑1 (MQ=255)                                                                                                                                                                                                               cGCTGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGTCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTAc                            >  1:253274/1‑251 (MQ=255)                                                                                                                                                                                                                  tggaaTTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGGCTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCCGGTAATTTATTCCCATAAAc                                      >  2:37983/5‑238 (MQ=255)                                                                                                                                                                                                                  tGGTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAAc                                      <  1:37983/238‑1 (MQ=255)                                                                                                                                                                                                                    gTATTCCGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTActgctg                       <  1:403683/251‑1 (MQ=255)                                                                                                                                                                                                                          cGAATACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCAc                 <  1:43134/251‑1 (MQ=255)                                                                                                                                                                                                                             aTACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCa              <  1:260947/251‑1 (MQ=255)                                                                                                                                                                                                                             aTACAGCGACAGCGGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGGGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGc               >  1:134605/1‑250 (MQ=255)                                                                                                                                                                                                                                  gcgacagcgGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGATGCTTCACGCAGGATc         >  1:291749/1‑251 (MQ=255)                                                                                                                                                                                                                                   cgacagctGTCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGATCt        <  2:158442/251‑1 (MQ=255)                                                                                                                                                                                                                                    gacagcgGGCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGATCTg       <  1:249606/251‑1 (MQ=255)                                                                                                                                                                                                                                          gggCTGTGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGATCTGAGttt  >  1:293340/1‑250 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    AAATCGCCATGCTCTGACAGGTTCCGGTATCTGTTAGTTGGTGCGGTGTACAGGCCCAGAACAGCGTGATATGACCCTGATTGATATTCACTTTTTCATTGTTGATCAGGTATTCCACATCGCCATCGAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTCTGCTTTTCAGGTAATTTATTCCCATAAACTCAGATTTACTGCTGCTTCACGCAGGATCTGAGTTT  >  NC_000913/4341343‑4341825

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