Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,871,153	T→G	N51T (AAC→ACC)	truD ←	tRNA pseudouridine(13) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,871,153	0	T	G	100.0%	149.8 / NA	42	N51T (AAC→ACC)	truD	tRNA pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (21/21);  total (21/21)

AGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTT  >  NC_000913/2870907‑2871393                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 ggggTAAAGGCGTTACCTGTCAGCGCGCCAAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCCCCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTGCGCCAGTGAATCCGCCACAAAACGGGTATTGCAGCCGgtttt                                                                                                                                                                                                                                              <  1:88798/250‑1 (MQ=255)   ggTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATAGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTg                                                                                                                                                                                                                                             >  1:155356/1‑250 (MQ=255)     tAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGAGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGCCCCGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTTCATCCGCCACAAAACGGGTATTGAGGCCGGTTTTGAgg                                                                                                                                                                                                                                          >  2:221990/1‑251 (MQ=255)             ttACCTTTCAGCGCGCCTAAACGCATCTTGCGCTTGTGCCTCGCATCCCCCAGCCCCTGGCATCCTTCCATTTGAAAGGCGCTCAGAGCGGGCATTTCCTTGCCCGGCACGCGAGCGCATACCCACTTTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCTGGAATTT‑TCAGGAATTTCGCCGGGGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTaa                                                                                                                                                                                                                                   <  1:372571/250‑1 (MQ=255)                        cgcgcCAAAACGCAGTTTGCGCTTGTGCCGCGCCTACTCCAGCCCCGGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTGGCCCGGCACGCGAGCGCATACCCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCAGGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCt                                                                                                                                                                                                                       <  1:384581/251‑1 (MQ=255)                            ccTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCAcc                                                                                                                                                                                                                   >  2:262797/1‑251 (MQ=255)                                         tGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCAAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGc                                                                                                                                                                                                      >  2:206507/1‑251 (MQ=255)                                         tGCGCTTGGGCCGCCCATACTCCACCCCCTGGCATCCTTCCAGTTGAAAGGCGCTCAGATCGGGCTTTTCCTTGCCCGGCACGCGAGCGAATACCCACTGTTCCGAAACAGCATTTTTGTCTTTTTCCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTGCGCAAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCa                                                                                                                                                                                                         <  2:330330/248‑1 (MQ=255)                                               ggtgCCGCGCATACCCCAGCCCCTGGCAGCCTCCCAGTTGAAAGGCGCTAAGATCGGGCATTTCCTTGCCGGGACCGCGAGCGAATAACCACTGTTCCGAAACAGCATTTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TAAGGAATTTCGCCAGGGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAg                                                                                                                                                                                                <  1:19399/250‑1 (MQ=255)                                                    cgcgCATACTCCACCACCTGGCAGCCTCCCAGTTGAAAGGCCCTCAGACCGGGATTTTCCTTGCCCGGCACGCGAGCGAATAACCATTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCaa                                                                                                                                                                                           <  1:81629/251‑1 (MQ=255)                                                     ccgcATACTCCAGCCCCTGGCAGCCTCCAAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCTCATAACCCTTGTTCCGTAACAGCATTTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCaaa                                                                                                                                                                                          <  1:261487/250‑1 (MQ=255)                                                                aGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTccaccac                                                                                                                                                                               >  1:80997/1‑251 (MQ=255)                                                                    ccTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGACTTT‑TCAGGAATTTCGCCAGTGCCTCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCa                                                                                                                                                                                                   >  1:373109/1‑227 (MQ=255)                                                                        gCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGAATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGt                                                                                                                                                                       >  1:404883/1‑251 (MQ=255)                                                                                ccAGTTGAAGGGCGCTCAGAGCGGTCATTTCCTTGCCCGGCACGCGAGCGAATAACCACTGTTCCGTAACGGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTTCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCTCCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCgg                                                                                                                                                                <  1:29942/250‑1 (MQ=255)                                                                                  attgtaaaGCCGCTCAGATCGGTCATTTCCTGGCCCGTCACGCGAGCTCATACCCACTGTCCCTTAACAGCATGATTGTCTTTTTCCCCAGCGAATCTGACTTCACGGGCA‑GAATTTGTCGGGAATTTCGCTAGGGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGTATTCTAACCAGAATATGCTCACCTTCTCCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCGTCCGGAtt                                                                                                                                                             <  1:388208/246‑1 (MQ=255)                                                                                   gTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAAACACTGTTCCGTAACAGACTGCTTGTCTTCTTGACCCGCGAAGCTGACCTCACGGGAATGAATTT‑ACAGGAATTCCGCCACAGCATCCGACACAAAACGGGTCTTGAAGACGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACAACCACAAAGTCTTCCGGATTg                                                                                                                                                            >  1:193326/1‑251 (MQ=255)                                                                                            cgcTCCGATCGGGCCTTCCCTTGCCGGGCACGCGCGCGAATAACCTCTTTTCCGACACCGCATGTTGTTCTTTTTCCCCAGCAAGGCGGACTACACGGGCATTAATTT‑CCAGGAATTTCGCCAGTCCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCAGCAGGCTCAAAGCCCAAATCTGCCACCACCACAAAGTCTTCCGGATTGGCTTTcagc                                                                                                                                                   <  2:46236/251‑1 (MQ=255)                                                                                                  gATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCg                                                                                                                                             >  1:194042/1‑251 (MQ=255)                                                                                                         cATTTCCTTGCCCGGCACGCGAGAGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGCATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACACAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCtt                                                                                                                                      >  2:68424/1‑251 (MQ=255)                                                                                                            ttCCTTGCCCGGCCCGGGAGCGCAAAACCCCTGTCCCGTAACAGCATGTTTGTTTTTTTGCCCAGCGAAGCTGACTCCACGGGAATGAATTT‑TCAGAAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGc                                                                                                                                    <  1:252028/250‑1 (MQ=255)                                                                                                             tCCTTGCCCGGCCCGCGAGCGCATAACCCCTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCg                                                                                                                                   <  1:256481/250‑1 (MQ=255)                                                                                                               cTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGtt                                                                                                                                >  1:135653/1‑251 (MQ=255)                                                                                                                 tGCCGGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCTTGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCCCGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGttt                                                                                                                               <  1:371091/250‑1 (MQ=255)                                                                                                                 ggccgtgCACGCGAGCCCATACCCTCTTTTCAGTAATATCAGTTTTGTTTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATGT‑ACAGGAATTGCGCCAGGGTATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCCCCTCCACAAAGTCTTCCGGATTGGCTTTCAGCAGCTCGGTGCCTTGCGGttt                                                                                                                               <  2:28225/244‑1 (MQ=255)                                                                                                                           gcgACCGCATACCCCCTGTTCCGTAACAGCATGTTTTTCTTTTGCCCCAGCGAAGCTGACTTCACGGGCATAAATTT‑TCAGGAATTTCGCCTGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCGTGCGGTTTACCGTGGAGGt                                                                                                                    <  1:364387/251‑1 (MQ=255)                                                                                                                                cccataccctcTGTTCCGTAACAGCATGTTTGTCTTTTGGCCCAGCAAAGCTGACTTCCCGGGCATGAATTT‑TCAGGAATTTCGCCAGGGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTg                                                                                                               <  1:127351/241‑1 (MQ=255)                                                                                                                                                   aaCAGCATGTTTGTCTTTTTCCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCa                                                                                            <  1:102139/251‑1 (MQ=255)                                                                                                                                                           gTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTttgttg                                                                                    >  1:98838/1‑251 (MQ=255)                                                                                                                                                               gTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCttt                                                                                                                                                       >  2:155787/1‑180 (MQ=255)                                                                                                                                                               gTCTTTTTGCCCAGCGAAGCTGACTTCACGGGAATGAATTT‑TCAGGAATTTCGCCAGGGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCttt                                                                                                                                                       <  1:155787/180‑1 (MQ=255)                                                                                                                                                                     ttGCCCCGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAACCGGGAATTGCAGCCGGTTTTGAGGCTTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCCAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAAGTCAATCATTTTGTTGCCTTAATg                                                                            >  2:55518/1‑249 (MQ=255)                                                                                                                                                                                 gCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGACACCGc                                                              >  2:308060/1‑251 (MQ=255)                                                                                                                                                                                         cACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAgg                                                      >  2:126622/1‑251 (MQ=255)                                                                                                                                                                                                       t‑tCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATAATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATACCTTccccc                                        >  1:388290/1‑250 (MQ=255)                                                                                                                                                                                                          cAGGAATTTCGCCATTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCAc                                       <  1:339863/250‑1 (MQ=255)                                                                                                                                                                                                                   cGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCa                                                                  >  1:319558/1‑214 (MQ=255)                                                                                                                                                                                                                      cAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGCGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATcc                          >  1:297861/1‑251 (MQ=255)                                                                                                                                                                                                                           cATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGtt                     >  2:170134/1‑251 (MQ=255)                                                                                                                                                                                                                            cgccGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGt                      <  2:289000/247‑1 (MQ=255)                                                                                                                                                                                                                               cGCCACAAAACGGGTATTGCTGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTc                 <  2:236617/251‑1 (MQ=255)                                                                                                                                                                                                                                    caAAACGGGTATTGCAGCCGGTGTTGAGGATTCTAACCAGAATATGCTCACCGTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAACTCCCAGTTTCTCCgtag            <  2:89574/251‑1 (MQ=255)                                                                                                                                                                                                                                              aTTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGCGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGTGTAGCGCGACCGCTTCACAGGCACTCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTCGTAGTGGCttt  >  2:98710/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 AGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTT  >  NC_000913/2870907‑2871393

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