Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	780,624	C→T	noncoding (71/76 nt)	lysT →	tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,624	0	C	T	96.5%	197.7 / ‑6.8	59	noncoding (71/76 nt)	lysT	tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (29/27);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AACGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGC  >  NC_000913/780378‑780867                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      aaCGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACtcacc                                                                                                                                                                                                                                                   >  1:266182/1‑251 (MQ=255)     cgATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATc                                                                                                                                                                                                                                               <  1:334037/251‑1 (MQ=255)          tAATGCATAACACCCGACTCGACGTCGCATTATTTCTGGTCGTGTCGCGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCACACAAACCGAAAGCAACGAAACAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACGCACCAATCGCTa                                                                                                                                                                                                                                           >  2:73321/1‑250 (MQ=38)             tGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGt                                                                                                                                                                                                                                     >  1:206396/1‑251 (MQ=255)                   acacacGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGc                                                                                                                                                                                                                               <  1:20346/251‑1 (MQ=255)                     acacGACCAGAAGTCGAATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCAA‑A‑GGTGGAATCgg                                                                                                                                                                                                                             >  2:55893/1‑251 (MQ=255)                      cacGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccaccagtctcaacggtgggagaggg                                                                                                                                                                                                                              >  2:448857/1‑231 (MQ=255)                           ccAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGgtagta                                                                                                                                                                                                                        >  1:585688/1‑250 (MQ=255)                              gAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAAc                                                                                                                                                                                                                    <  2:479239/251‑1 (MQ=255)                                 gTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTg                                                                                                                                                                                                                 >  1:27309/1‑251 (MQ=255)                                       ttattTCTGGTCGTGTCGCGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGAt                                                                                                                                                                                                           >  1:461001/1‑251 (MQ=255)                                           ttCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGGCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAAc                                                                                                                                                                                                        <  2:450863/250‑1 (MQ=255)                                            tCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACg                                                                                                                                                                                                       <  1:25646/250‑1 (MQ=255)                                               ggTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCTACGGGATATCAAACAAACCGAAAGCAACGAAAAATTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGTCTTGTAATCAATTGTTCGTAGGTTCGAATCCGGCACGATTCACCAATCGCTG‑A‑GGGGGAAGCGGTAGTAAAACGTGAAGGATAACGtt                                                                                                                                                                                                     <  1:29245/249‑1 (MQ=255)                                               ggTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTg                                                                                                                                                                                                    <  1:416617/250‑1 (MQ=255)                                                gtcgtgtcgtGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGCGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGc                                                                                                                                                                                                   >  1:239525/1‑250 (MQ=255)                                                  cgtgtagtgagcctCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTAGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCTATA‑A‑GGTGGAGGCGGTAGTAAAACGTGAAGGATAACTTTGCAtt                                                                                                                                                                                                >  2:488467/14‑250 (MQ=38)                                                   gtgtcgtgCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATg                                                                                                                                                                                                >  1:375209/1‑250 (MQ=255)                                                        gtgCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAAc                                                                                                                                                                                          >  1:580041/1‑251 (MQ=255)                                                                                      aTCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCAt                                                                                                                                                            >  1:288175/1‑251 (MQ=255)                                                                                      aTCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGAAGTAAAACGGGAAGGATAACGGTGCATGAGCAACGGCACGAAGGGTGAAACGAAGTCGAGTCa                                                                                                                                                             >  2:585955/1‑250 (MQ=38)                                                                                        cGAAATTTTTGTTGCGCTCAAGGCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATcc                                                                                                                                                          <  1:364954/251‑1 (MQ=255)                                                                                           aaTTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGTATGAGCAACGGTCCGAAGGGCGAGACGAAGTCGAGTCATCCttc                                                                                                                                                       <  1:210755/251‑3 (MQ=255)                                                                                              ttttGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCAc                                                                                                                                                     <  1:424153/250‑1 (MQ=255)                                                                                                 tgttgCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGAc                                                                                                                                                  <  2:27309/250‑1 (MQ=255)                                                                                                        cTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAACCAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGACGTCGAGTCATCCTGCACGACCCACCACt                                                                                                                                          >  2:525990/1‑251 (MQ=255)                                                                                                              tCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACAt                                                                                                                                     >  1:449915/1‑250 (MQ=255)                                                                                                               cTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGAGGAAGTCGAGTCATCCTGCACGACACACCACTAACAt                                                                                                                                     >  2:67132/1‑249 (MQ=255)                                                                                                                      cAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATagttagtt                                                                                                                             <  1:198952/250‑1 (MQ=255)                                                                                                                       aGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTgt                                                                                                                           <  1:7725/251‑1 (MQ=255)                                                                                                                                 gccgccGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGc                                                                                                                 >  1:384098/1‑251 (MQ=255)                                                                                                                                       gTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCgtagta                                                                                                           <  1:105660/251‑1 (MQ=255)                                                                                                                                           ccACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCgg                                                                                                       >  2:473054/1‑251 (MQ=255)                                                                                                                                            cACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACGACGTG‑TAGA‑GGAGGAGGCGGGAGTTAGGCGTGCAGGGTAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTGAGTTGTAGTATCCAGCGTAGTAGCggg                                                                                                      >  2:226088/1‑251 (MQ=21)                                                                                                                                               gggATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGa                                                                                                   <  1:537472/251‑1 (MQ=255)                                                                                                                                                     tCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGGAGTATCCAGCGTAGTATCGGGTGATTAGcc                                                                                             >  2:579555/1‑250 (MQ=255)                                                                                                                                                      caaacaaaCCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCTATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCt                                                                                             <  1:317223/250‑1 (MQ=255)                                                                                                                                                               cGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCAGCAGTCGATA‑A‑GTTGGAAGCGGGAGTGAAACGGGATGAATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACACACCACTAACATAGTTAGTTGTAGTATCCAGCGTATTATAGGGTGATTAGCTCCGCTggg                                                                                     >  2:414656/1‑249 (MQ=25)                                                                                                                                                               cGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGGGAAGGCTAACGTTGCATGAGCAACGGCCCGAAGGGCGAGAAGAAGACGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGAGTAGTATAGGGTGATTAGCTCAGCTggg                                                                                     >  1:287308/1‑249 (MQ=255)                                                                                                                                                                gAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGgaga                                                                                  <  1:421817/251‑1 (MQ=255)                                                                                                                                                                gAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCAAGAACGGCTA‑A‑GGAGGAAGCGGAAGTAAAACGAGAAGGATAACGTTGCATGAGCAACGGGCCGAAGGGAGAGACGAAGTCGAGACATCCTGCACGACACACCAAAAACATAGTTAGTAGGAGTATCCAGAGTAGTAACGGGAGAGTAGCTCAGCTggg                                                                                     >  1:37466/1‑248 (MQ=21)                                                                                                                                                                    gCAACGAAAAAGGGGGTCGGTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAACCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACAGAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCAc                                                                              <  2:35583/251‑1 (MQ=255)                                                                                                                                                                     cAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGTGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGTTTGATTAGCGCAGCTGGGAGAGCAcc                                                                             >  1:10759/1‑251 (MQ=255)                                                                                                                                                                      aaCGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCt                                                                            <  2:21031/251‑1 (MQ=255)                                                                                                                                                                               gTGGGTCGTTAGCTCAGTTGGTAGATCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACaa                                                                   <  2:267628/251‑1 (MQ=255)                                                                                                                                                                               gTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACaa                                                                   <  2:288175/251‑1 (MQ=255)                                                                                                                                                                               gTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAAT‑TTTA‑ACGGTGGAGACGGTTGTAAAGCGTAAAGCGTAACGGTGCCGGATCGGCAGCCAGAGGAGCGCTAGGAAGTCGCGTCATCATGAACGACACACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACa                                                                    >  2:223145/1‑250 (MQ=12)                                                                                                                                                                                 gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGAGTAGCTCAGCTGGGAGAGCACCTCCCTTACAagg                                                                 <  1:525990/251‑1 (MQ=255)                                                                                                                                                                                      gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAgggg                                                            <  1:289113/251‑1 (MQ=255)                                                                                                                                                                                             cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAAGGCGAGAAGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTcggcg                                                     <  2:405698/251‑1 (MQ=255)                                                                                                                                                                                                   gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGa                                               <  2:188166/251‑1 (MQ=255)                                                                                                                                                                                                               aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTTCAAGGAGGGGGTCGGCGGTTCGATCCCGTCATcac                                   >  1:69513/1‑251 (MQ=255)                                                                                                                                                                                                               aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGCGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGGAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATcac                                   >  1:237491/1‑251 (MQ=255)                                                                                                                                                                                                                cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGGTCGATCCCGTCATcacc                                  >  1:243789/1‑251 (MQ=255)                                                                                                                                                                                                                         aaTTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTAACCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGGTCCCGTCATCACCCACCACCgg                         >  1:465117/1‑251 (MQ=255)                                                                                                                                                                                                                             ggTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGAGTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTc                      <  2:585688/250‑1 (MQ=255)                                                                                                                                                                                                                                  cAGGTTCGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGTGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTccgtagc                >  1:509483/1‑245 (MQ=255)                                                                                                                                                                                                                                        cGAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGtt          <  1:473054/251‑1 (MQ=255)                                                                                                                                                                                                                                         gAATCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTg         >  1:426261/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTCCTGCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCCCCACCGGGTCGTTAGCTCAGTTGGt       >  1:129789/1‑251 (MQ=255)                                                                                                                                                                                                                                                gCACGACTCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCCCCACCGGGTCGTTAGCTCAGTTGGTAGAGc  >  1:142148/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      AACGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAAGTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTA‑A‑GGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGC  >  NC_000913/780378‑780867

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