Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,377,545	C→T	V60I (GTC→ATC)	blc ←	outer membrane lipoprotein Blc

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,377,545	0	C	T	100.0%	143.2 / NA	42	V60I (GTC→ATC)	blc	outer membrane lipoprotein Blc
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/23);  total (19/23)

GGTAGTCGCGGTCCGGGCCGCAAACCAGCGCATGGCGGTATTCCCGATCGAGTGCAATAACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGACTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCC  >  NC_000913/4377311‑4377787                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   ggtggtcgcggcccggccgccacaCCAGCGCAGGGAGGTTTCCCCGATGGAGCGAAATAACTTTAAAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGt                                                                                                                                                                                                                                    <  2:240277/228‑1 (MQ=255)               gggCCGCAAACCAGCGCATGGCGGTATTCCCGATCGAGTGCAATAACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATc                                                                                                                                                                                                                      >  1:90348/1‑251 (MQ=255)               gggCCGCAAACCAGCGCATGGCGGTATTCCCGATCGAGAGCAATAACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGCGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGAt                                                                                                                                                                                                                       >  2:297280/1‑250 (MQ=255)                      aaaccggcgcggggCGTTATCCCCGATCGTGTACAATAAGTTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACgg                                                                                                                                                                                                               <  1:45131/239‑1 (MQ=255)                                                gggAGTGAAATAACGTTAAAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATa                                                                                                                                                                                     <  2:137568/249‑1 (MQ=255)                                                    gtacattatctttATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCATTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATa                                                                                                                                                                                     <  1:337557/236‑1 (MQ=255)                                                     gcaaaaactgtataatcctccATAGAAAGGACCAAAGAATGACTCTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAATTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGtt                                                                                                                                                                                <  1:116075/232‑1 (MQ=255)                                                        aatattagtaaacccGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTGCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCGTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTgg                                                                                                                                                                             <  2:36810/238‑1 (MQ=255)                                                            acgtaaaatccGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTc                                                                                                                                                                         <  1:276748/242‑1 (MQ=255)                                                                tataACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGaa                                                                                                                                                                      >  1:344060/1‑250 (MQ=255)                                                                   aaCCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAatta                                                                                                                                                                  <  1:166199/251‑1 (MQ=255)                                                                      cGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTAttt                                                                                                                                                               >  2:290631/1‑251 (MQ=255)                                                                               aaGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGGCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCAGACCAGGTACCGAGATAACGTTTGGCGTCGATATTATTTACTACGGTc                                                                                                                                                      <  1:150285/251‑1 (MQ=255)                                                                                aGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTca                                                                                                                                                     <  2:139849/251‑1 (MQ=255)                                                                                  gACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTcacg                                                                                                                                                   >  1:56567/1‑251 (MQ=255)                                                                                        aGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGGCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTcacgccacgc                                                                                                                                             >  2:394576/1‑251 (MQ=255)                                                                                          aaTGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACgcgg                                                                                                                                           >  1:13338/1‑251 (MQ=255)                                                                                                acTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCGCAGTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTTCCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGt                                                                                                                                     <  2:385020/251‑1 (MQ=255)                                                                                                     cagcgcagcgCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGgag                                                                                                                                <  2:179342/251‑1 (MQ=255)                                                                                                             cgcgAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAg                                                                                                                        <  2:397900/251‑1 (MQ=255)                                                                                                               cgAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGc                                                                                                                                                                            >  1:63758/1‑197 (MQ=255)                                                                                                               cgAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGc                                                                                                                                                                            <  2:63758/197‑1 (MQ=255)                                                                                                                                   gtgcGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTg                                                                                                  <  1:376031/248‑1 (MQ=255)                                                                                                                                           tCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCg                                                                                          >  1:368360/1‑251 (MQ=255)                                                                                                                                                 cACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTCATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCNACg                                                                                    >  2:408523/1‑251 (MQ=255)                                                                                                                                                           ccACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGc                                                                           >  1:321030/1‑250 (MQ=255)                                                                                                                                                             acaTTCCTCTGTCAGGGATATCGACATTATTAATGACATTCAGGCCGCCGTCATCACGCAGACTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAGGCGCTGATCAAAACGGGCAATCACATACCAGGTACGGAGATCACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCCCGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAACAgg                                                                        >  2:60761/1‑251 (MQ=255)                                                                                                                                                              caTTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGAAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGgc                                                                       >  2:192212/1‑251 (MQ=255)                                                                                                                                                                                       ttattaATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCACTCCTCGTTCAAAGCGGTGATCAAATCGGGCAATCTCATACCAGGTGCCGACATAACGTCTGGCGTCGAAATTATTTACTACGGTCACGCCACGCTGCGGCGTAGGAGAACTGCAGGCCACGACCAGAAATGCAGCTTTCGCTGCTGCAACGAGAGGGAGCAGGCGCATAAATGTTGCCTTACTgttt                                                <  2:154237/249‑4 (MQ=255)                                                                                                                                                                                           tactGACATTACGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGGACCGAGATAACGTTTGGCGTCGAAATTATTTACTACAGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAggg                                                                              >  2:64257/4‑215 (MQ=255)                                                                                                                                                                                           taATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAggg                                                                              <  1:64257/215‑1 (MQ=255)                                                                                                                                                                                           taATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGttttttt                                           <  2:90348/250‑1 (MQ=255)                                                                                                                                                                                                        gccgccGTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGa                             <  1:408523/251‑1 (MQ=255)                                                                                                                                                                                                           gccgTCATCACGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGCGTAGATGAc                            >  1:192813/1‑249 (MQ=255)                                                                                                                                                                                                              gTCATCACGGAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCATCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCaa                       <  1:79170/251‑1 (MQ=255)                                                                                                                                                                                                                  tcaCGCAGGATGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGa                    >  2:253721/1‑250 (MQ=255)                                                                                                                                                                                                                   caCGCAGGCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGaaa                  <  1:282723/251‑1 (MQ=255)                                                                                                                                                                                                                         ggCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCgg             >  2:147973/1‑250 (MQ=255)                                                                                                                                                                                                                          gCTGTATGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCgg             >  1:169255/1‑249 (MQ=255)                                                                                                                                                                                                                               aTGTTGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGGGTAGATGACAGCAAGGAAAAGCGGAGaaa        >  2:220037/1‑249 (MQ=255)                                                                                                                                                                                                                                   tGCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTc   <  2:169255/250‑1 (MQ=255)                                                                                                                                                                                                                                    gCGGTGATTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTcc  <  1:394576/250‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   GGTAGTCGCGGTCCGGGCCGCAAACCAGCGCATGGCGGTATTCCCGATCGAGTGCAATAACGTTATAACCGCCATAGAAAGGACCAAAGAATGACACTTTCAGCGCAGCGCGAGTTGGTGCGCCGGTAAAGTACGCTTTCCCTTCACTCTGCTGCCACATTCCTCTGTCAGGGTTATAGCCTTTATTAATGACATTCAGGCCGCCGTCATCACGCAGGCTGTATGTTGCGGTGACTTTTTCCAGTCCACGTTCAAAGCGGTGATCAAAACGGGCAATCTCATACCAGGTACCGAGATAACGTTTGGCGTCGAAATTATTTACTACGGTCACGCCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCC  >  NC_000913/4377311‑4377787

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: ‑