Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	694,721	C→T	G79S (GGC→AGC)	miaB ←	isopentenyl‑adenosine A37 tRNA methylthiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	694,721	0	C	T	100.0%	205.6 / NA	58	G79S (GGC→AGC)	miaB	isopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/36);  total (22/36)

TATTTATTGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCT  >  NC_000913/694477‑694966                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      tatttattGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAAt                                                                                                                                                                                                                                                 >  2:105203/1‑251 (MQ=255)      attGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGt                                                                                                                                                                                                                                            >  1:94925/1‑251 (MQ=255)        tGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCGACAACCGGGCTGCGGTCGCCACGCACGGGGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATAGCTACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCTTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTc                                                                                                                                                                                                                                           <  1:268060/250‑1 (MQ=255)         gCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTg                                                                                                                                                                                                                                         >  2:330344/1‑251 (MQ=255)             ccTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGAtt                                                                                                                                                                                                                                     <  2:90194/251‑1 (MQ=255)             ccTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGGTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGAtt                                                                                                                                                                                                                                     <  2:329843/251‑1 (MQ=255)                       ggggtaCAAACACGATCGGCCAGTCGGCGCGCGGTTCAGGTAGACAGTCAAACTTCACGCGTTCCGGGAAGCTGATATCTACAACCGGGCTGCTGTCGCCACGCACTGAGTTGATCATCTCCGGCAGACGGTGCAGCGTGTGCGGCACAAAAATAATATCGCCATAGTGGGCGCGCTTGCGAATGTG‑CTCGCCTTCTTGCGTTGCCACGAAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCtttt                                                                                                                                                                                                                           <  1:246088/246‑1 (MQ=255)                         gagaCAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATACTCAGGTCTGGATTCTTCTCTTTTaa                                                                                                                                                                                                                         <  1:342460/251‑1 (MQ=255)                          agaCAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAAc                                                                                                                                                                                                                        <  2:38987/251‑1 (MQ=255)                            aCAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTGCTCTTTTAACAg                                                                                                                                                                                                                      <  2:133164/251‑1 (MQ=255)                             cAAACGCTGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGGTTGCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGTAGTTGATCATCTCCGGCAGACGGTGCAGCGTGTGCGGCCCAAAAATAATATCGACATAGTGGTCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGAGTCCTCGCAGCGACCGACGCTGATAATCAGGTTTGGATTCTTCTCTTTTAACAGt                                                                                                                                                                                                                     <  2:291280/251‑1 (MQ=255)                             cAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGt                                                                                                                                                                                                                     >  1:58169/1‑251 (MQ=255)                               aaCGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCCCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGttt                                                                                                                                                                                                                   <  1:10660/251‑1 (MQ=255)                                aCGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTc                                                                                                                                                                                                                  >  1:152292/1‑251 (MQ=255)                                       cGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCTAGCGAc                                                                                                                                                                                                           >  2:371895/1‑251 (MQ=255)                                          cccTTCGGCGCGCGGTTCCGGTCGCCGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGCGCAGCGTTTGCGGCCCACAAATAATATCGACATGGTGGGCGCGCTGGCGAATGTG‑CTGGCCTTCTTGCGATGCCACGCAGACACCGCCGCTGATAATCAGGTCTGGGTGCTTCTCTTTTAACCGTTTCACGCGAccc                                                                                                                                                                                                         >  2:203815/1‑250 (MQ=255)                                              tCGGCGCGCGGTTCCGTTAGACGGTCAAACTTCTCGATTTCCGGGAATCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACAGAGTGGGCGCGCTGTCGAATGTG‑CTCGCCTTCTTGCGATTCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACt                                                                                                                                                                                                     <  1:407039/250‑1 (MQ=255)                                                  cgcgcgGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATgg                                                                                                                                                                                                <  1:12863/251‑1 (MQ=255)                                                       ggTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTAAAACCGGGCTGCGGCCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCCGCGTTTGCGGCCCACAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGAc                                                                                                                                                                                           >  1:390721/1‑251 (MQ=255)                                                                 gACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATATCCAACCAAGCTCCGCTCGCCACGCACGGAGCTGATCATCTCCCGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGCGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATACTCAGGTCTGGATTCTTCTCTTTTAACAGCTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTg                                                                                                                                                                                  >  1:412002/1‑250 (MQ=255)                                                                                    ttCCGGGCAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGc                                                                                                                                                              >  1:118135/1‑251 (MQ=255)                                                                                            aGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTc                                                                                                                                                      <  1:379284/251‑1 (MQ=255)                                                                                                gagatCTACAACCGGGCTGCGGGCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTGTGCCGCCCAAAAATAATATCGACATAGTGGTCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCATGTCTGGATTCTTCTCTTTTAACAGGTTCCAGCTACCCAACTGATTGAAGACGTTTGCCTGAGCCCTCTCGCGTATTAAGCAGGTGTTcag                                                                                                                                                    <  2:390721/246‑1 (MQ=255)                                                                                                gATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcag                                                                                                                                                    <  1:216872/249‑1 (MQ=255)                                                                                                   atCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTcagcagc                                                                                                                                                <  1:326139/250‑1 (MQ=255)                                                                                                       aCAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCCAAACTAATATCGACATCGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGCTAATCAGGCCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCagccgcagccc                                                                                                                                           >  2:246960/1‑249 (MQ=255)                                                                                                       aCAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCagcagcagc                                                                                                                                             <  1:233727/249‑1 (MQ=255)                                                                                                               gCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGctt                                                                                                                                   <  1:313604/251‑1 (MQ=255)                                                                                                                  gCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGcttctt                                                                                                                                <  1:98570/251‑1 (MQ=255)                                                                                                                    ggTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGGGTTCAGCCGCAGCACATCCGCTTCTTcc                                                                                                                              >  2:380989/1‑251 (MQ=255)                                                                                                                      gcGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGc                                                                                                                            <  2:344715/250‑1 (MQ=255)                                                                                                                           acgcacgGCGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCACAAAATAATATCGACATAGTGGGCGCGCTGGCGACTGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGt                                                                                                                       >  2:158740/1‑251 (MQ=255)                                                                                                                               cctgAGTTGATCATCTCCGGCAGACGTTGCAGGTTTTGCGGCCCAAAAATAATATCGACATAGTGGGGGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCGGCAGCAGCACATCCGTTTCTGCCGCCACGTCGGt                                                                                                                   <  2:302509/248‑1 (MQ=255)                                                                                                                                                 ggCAGACGGTGCCGCGTTTGCGGCCCAAAAATAATATCGACAGAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtgggt                                                                                                 <  2:152292/251‑1 (MQ=255)                                                                                                                                                 ggCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGGCCCAGCTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtgggt                                                                                                 <  1:269689/251‑1 (MQ=255)                                                                                                                                                 ggCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtgggt                                                                                                 <  1:15819/251‑1 (MQ=255)                                                                                                                                                  gCAGACGGTGCAGCGTGTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTGTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtgggtg                                                                                                <  2:257076/251‑1 (MQ=255)                                                                                                                                                   cAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtgggtgg                                                                                               <  1:11776/251‑1 (MQ=255)                                                                                                                                                    aGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTGGATAGCCGTGGGTGGc                                                                                              <  1:371754/251‑1 (MQ=255)                                                                                                                                                    aGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCCCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGGCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGGCAGTTGATAGCCGTGGGTGGc                                                                                              >  1:1849/1‑251 (MQ=255)                                                                                                                                                        ggTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATcc                                                                                          <  1:246960/251‑1 (MQ=255)                                                                                                                                                           gCAGCGTGTGCGGCCCAAAAATAATATCGACATAGGGGGCGCGCGGGCGAATGTG‑CTCGCCTTCTGGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTGTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTGCCGCCACGTGGGTCAGTTGATAGCCGTGGGTGGCATCcagc                                                                                       <  2:174517/251‑1 (MQ=255)                                                                                                                                                             aGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCcagc                                                                                       >  2:61629/1‑249 (MQ=255)                                                                                                                                                                    gCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGcc                                                                              <  1:303900/251‑1 (MQ=255)                                                                                                                                                                         ccAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGGTTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTGTCCTGGGCCTTCTCGCGGATTGAGCAGGTGTTCTGCATCAGCAGATCCGCTTCTTCCGCCACGTCGGGCAGGGGATAGCCGTGGGTTGCATCCAGCAGATCGGCCATCtt                                                                         <  1:413505/251‑1 (MQ=255)                                                                                                                                                                         ccAAAAATAATATCGACATAGGGGGCGCGCTGGCTAAT‑TGTCTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCtt                                                                         <  2:244509/251‑1 (MQ=255)                                                                                                                                                                                  atatCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGCTGAATc                                                                >  1:220185/1‑251 (MQ=255)                                                                                                                                                                                       gACATAGAGGGCGCGCTGGAGAATGTG‑ATCGCCTTATTGCGATGCCAAGAAGCAAACGACGCTGATAATCAGGTCTGGAATCTTCTCTTTTAACAGTTTACAGCGCCCCCACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCAGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACt                                                           >  2:217011/1‑251 (MQ=255)                                                                                                                                                                                           tAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTCGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGAACTCGtt                                                       >  2:417897/1‑251 (MQ=255)                                                                                                                                                                                                     gcTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGAc                                               >  1:345220/1‑249 (MQ=255)                                                                                                                                                                                                         gCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGcccc                                         <  1:61629/251‑1 (MQ=255)                                                                                                                                                                                                          cGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCa                                        <  1:330344/251‑1 (MQ=255)                                                                                                                                                                                                               gtg‑CTCGCCTTCTTGCGATGCCAAGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGCTTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGCATCGTACTCGTTCATCTGACAGCCACAGGttt                                   >  2:257729/1‑251 (MQ=255)                                                                                                                                                                                                                gtgttggctttgttGCGATGCCCCGCAGCCCCCAACGGTGGTAATCAGGGCTGGATTCTTCTCTTTGAACAGTTTCCAGCGATCCAACGGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTTAGCAGGGGTTAAGCAGCAGCACCGCCGCTTCTTACGCCACGTCGGTCAGTTGATAGCCTTGGGTGGCATCCAGCAGATCGGCCATCTGCTATGAATCGTACTTGTTCATCTGACAGCCCCAGGttt                                   <  2:265270/239‑1 (MQ=255)                                                                                                                                                                                                                             tGCGAGGCCACGCAGCCACCGACGCTGGTAATCAGGGTTGGTTTCTTCTCTTTTAACAGTTTCCAGCGTCCCAACTGAGGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCTATGAATCGTTCTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAg                         <  2:374057/248‑1 (MQ=255)                                                                                                                                                                                                                                 aTGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTg                  <  2:118135/251‑1 (MQ=255)                                                                                                                                                                                                                                       cGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATc            >  2:110876/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGctct  >  2:334601/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      TATTTATTGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGTG‑CTCGCCTTCTTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCT  >  NC_000913/694477‑694966

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