Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	696,470	C→T	noncoding (35/75 nt)	glnX ←	tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	696,470	0	C	T	97.9%	141.1 / NA	47	noncoding (35/75 nt)	glnX	tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/21);  minor base G (0/1);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGA  >  NC_000913/696232‑696711                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  ggctggtTGAAACGTCTGTTGCTGAAATATTCGTGATGGTTGTAGCCTTACCACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCACATCCTACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAATCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑ACCCTCGGAATGCCGGAATTAGAATCCGGTGc                                                                                                                                                                                                                                       <  1:278637/246‑1 (MQ=255)       ttGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCCACCGTAGGCCGGATACGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGACTTAGAATCCGGTGCCTTAcc                                                                                                                                                                                                                                 >  1:294660/1‑251 (MQ=255)              gTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAA‑CCCGCCCGAGCGAGCTTTTTTAATTGGGCTGGGGTAAGAGGGTTCGGAACC‑CGGAATGCCGGGATTAGaaatcggggcctttaccccttgcg                                                                                                                                                                                                                         >  1:536346/1‑229 (MQ=255)                         aaaTATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAAc                                                                                                                                                                                                               >  1:8910/1‑251 (MQ=255)                          aaTATGCGTTAGGGTTGTAGCCTTATAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGGCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAATCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGTTTCGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACt                                                                                                                                                                                                              <  1:623198/251‑1 (MQ=255)                                cGTTAGGGTTGTAGCCTTACAACATTGCGGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACAGTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCTAAGCGAGCGTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGc                                                                                                                                                                                                        <  1:265218/251‑1 (MQ=255)                                cGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGGCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGc                                                                                                                                                                                                        <  1:114625/251‑1 (MQ=255)                                  ttAGGGTTGTAGCATTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGTGTGccg                                                                                                                                                                                                      >  2:110757/1‑249 (MQ=255)                                          gTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGGCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTGGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAgg                                                                                                                                                                                              <  1:608928/251‑1 (MQ=255)                                           tAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCCCTTACAAGGt                                                                                                                                                                                             >  2:540526/1‑251 (MQ=255)                                              cctgacaacaTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTGTTTAATGTGGCTGGGGTACGATGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAg                                                                                                                                                                                          <  1:560915/247‑1 (MQ=255)                                                       aTTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCATAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGaa                                                                                                                                                                                 <  1:383875/251‑1 (MQ=255)                                                         tGCCGGATGCGTGCCACCCGTAGGTCGGCTAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGCATTAGAATCCGGTGCCTTACCGCTTGGCGAt                                                                                                                                                                                                                       >  1:201836/1‑211 (MQ=255)                                                         tGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGAt                                                                                                                                                                                                                       <  2:201836/211‑1 (MQ=255)                                                             ggATGCGGGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCTAAATTGAAGGATAAGACGTGTCATCATCGCATTTGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCAGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTTGCGATACCCCAACTGGGTGCCCTTGCAAGGGAAGCGTCTTGAATAAAtt                                                                                                                                                                           <  2:108595/251‑1 (MQ=38)                                                                     gCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTggg                                                                                                                                                                    <  1:481160/250‑1 (MQ=255)                                                                              aGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGCCGGATAAGACGTGCCAACATCGCATTCGACATTGAATGAAAGCAGAAAAGCAAAAAGCTCGTCGAAGCGAGCATTTTTAATGTGGCTGGGGAACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGCTACCCTAACTGGGGGCAGTTACATGGTGAGCGCCATGACTAAATAGGGGGGGGGACGAGGAt                                                                                                                                                          >  2:232641/1‑251 (MQ=25)                                                                                        aGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGAACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCg                                                                                                                                                >  1:504093/1‑251 (MQ=255)                                                                                                cAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGGGCCCTTTCAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTAGGAATGCCg                                                                                                                                        >  1:465406/1‑251 (MQ=255)                                                                                                 aGCGGCGCATCCGACATTGAAGGATAAGACGTGTCAACATAGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTGTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAAGGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCgg                                                                                                                                       <  1:379996/251‑1 (MQ=255)                                                                                                   cGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGCATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCAGGGGTACGAGGATTCG‑AACATCGGAATGCCGGAATTAGAATCCAGTACCTTACCGCTTGGCGATACCCCAAATGAGTGACCTTAAAAGGTAAGCGTATTGAATAAATTGGATAGGTTGCGAGGATTCGAACATCGTAAAGCCGGa                                                                                                                                      >  1:621352/1‑250 (MQ=255)                                                                                                    gTCGCCTCCGACATTGAAGGATACGACTTGGCAACAGCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGGTCGGCGAAGTGAGCTTGGTTAATGTGGCTGCGGTACGAGGAGCCG‑AACATCGGAATGCCGGAAGTAGAATCCGGAGCCTTAACCCTTGGCGATGCCCCAACTGCGTGCACTTACACGGTAAGCGTCTTGAATAAAGTGGCTGGTGTACGAGGATTCGAAccccgtaccgctggaat                                                                                                                                    <  2:277590/251‑16 (MQ=255)                                                                                                                               gtgtCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCg                                                                                                         <  1:483560/251‑1 (MQ=255)                                                                                                                                    aacgtCGCATTCGGCATTGAAGTAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGTAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGACTATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCGGGGGTACGAGGATTCGAACCTCGTAATGCCGGAATCAGAATCCGGTGCCTGACCGCTTGGCGATAcc                                                                                                    <  1:232641/247‑1 (MQ=255)                                                                                                                                        tCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGACTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGACTAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCca                                                                                                 >  1:371322/1‑250 (MQ=255)                                                                                                                                                 aCATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGGACTTACACGGTAAGCGGCTTGAATAAATAGGCTGTGGTCCGAGGAATAGAACCTAGGAATGCCGGTATCAACATACGGTGCATTACCGCATGGCGATACCCCACCAAAtc                                                                                          >  2:478911/1‑247 (MQ=25)                                                                                                                                                  cagtgaatgaaCGCAGAAAAGCAAAAAGCAGGCCGAAGCGAGCTTGTTTAATGTGGCGGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATGAGAATCCGGGGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGtt                                                                                      <  1:310671/248‑1 (MQ=38)                                                                                                                                                     tgaatgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGtttt                                                                                    <  1:220808/250‑1 (MQ=255)                                                                                                                                                         tgaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGCATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACACATTGGTTTTGAAtt                                                                               >  1:238775/1‑251 (MQ=255)                                                                                                                                                          gaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACATGGTGAGCGTCTTGAATAAATTGGCTGGGGCTCGAGGATTCGAACGTCGGAAAGCCGGAATCAGGATCCGTCGCCTTACCGCTTGTCGATACTCCAACAAATGGGTTTTgtcttt                                                                              >  1:7282/1‑246 (MQ=37)                                                                                                                                                          gaaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAGATTGGTGTTGAAttt                                                                              >  1:177832/1‑251 (MQ=255)                                                                                                                                                           aaCGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGGAAGCGCCTTGAATAAATTGGCTGGAGTTCGCGGATTCGAACCTCGGAATGCCGGCATCAGACTCCGGTTCCTTACCGCTTCGCGATACCACAAAAACTTGGTTTTGAATTTg                                                                             >  2:325578/1‑251 (MQ=37)                                                                                                                                                               cAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGACGGAGTTAGAATCAGGTGCCTTACCGCTTGGCGATACCCCAACTGGCCGCCCTTCCAAGGCAACCGCATTGAATAAATTGGGTGGAGTACGAGGATACGAACCTCGGACTGCCGGAATTAGAAACCGGTGCCTTACCGCTTGGCGCCACCCCAACAAACTGGCTTTGGATCTGCCg                                                                          >  1:282080/1‑250 (MQ=21)                                                                                                                                                                aGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGa                                                                         >  1:523017/1‑250 (MQ=255)                                                                                                                                                                 gAAAAGCAAAAAGCTCGACGAAGCGAGCTATTTTAATGTGGCTGGGGTACGAGGATTCG‑AGCCTCGGAATGCAGGAATTAGAATCCGGTGCCTTACCGCATGGCGATACCCCAGCTGGGTGCACTTACAAGGTAAGCGGCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCCGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAAc                                                                       <  1:46438/251‑1 (MQ=38)                                                                                                                                                                   aaaGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAAGTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACat                                                                     <  1:241301/251‑1 (MQ=255)                                                                                                                                                                                    gAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGCTTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGCACATATTCGATCCATTCAGa                                                     >  1:306342/1‑250 (MQ=255)                                                                                                                                                                                        cGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGGAAGCGTCTTGTATAAATTGTCTGGGGTACGAGGCTTCGAACCTCGGAATGCCGGAATCAGAATCCGGGGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGATTTTGCCGCACCTATTCGATACATTCAAAATTtg                                                >  2:471729/1‑251 (MQ=38)                                                                                                                                                                                                         gCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCTACTGGGTGCCCTTGCAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGTGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTTCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGAt                                <  1:237925/250‑1 (MQ=255)                                                                                                                                                                                                          cTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGAACTTACAAGGTAAGCGTCTTGACTAAATTGTCTGGGGTACGAGGATTCGAACCTCGGAGTACCGGAATCAGAATCCGGTGCCTTTCCGCTTGGCGATACACCAAAAAATTGGGGCTGATTTTGGCGTACATCTTCGAGACATTCAGATTTTGGTGGTTACGTCGGGAtc                               >  1:359698/1‑249 (MQ=25)                                                                                                                                                                                                           tGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGCATTTGGTGGCTACGACGGGATTCg                             >  1:536655/1‑251 (MQ=255)                                                                                                                                                                                                                         tCG‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGAcccc                <  1:471729/250‑1 (MQ=255)                                                                                                                                                                                                                           g‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTCAGCGTCTTGAATAAATTGGCTGGGTTACGCGGATTCGACCCTCGGAATGCCGGAATCAGAATCCGGTGCCTTCCCTCTTGGCGATCCCCCAACAACTTGGTTTTGAATTTTCCGCACATCTTCGCTACATTCAGACTTTGGTTGCTACGTCTGGATTCGCACCTGTGACCCcatc             >  1:444647/1‑251 (MQ=255)                                                                                                                                                                                                                           g‑AACCTCGGAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGAAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGCTTCGAACCTCGGAATGCCGGAATCAGAATCCGGGGCATTACCGCTTGGCGATACACCAACAAATAGGATTAGAATTTGTCGAACAGAGTCGATATCTTCAGAATTTGGTGGTTACTACGGGTGTTGAACTTGTGACCCcatt             >  1:534838/1‑250 (MQ=25)                                                                                                                                                                                                                                   ggAATGCCGGAATTAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATg       >  1:350602/1‑250 (MQ=255)                                                                                                                                                                                                                                      aTGCCGGAATTAGAAGCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTg   <  1:512836/251‑1 (MQ=255)                                                                                                                                                                                                                                        gCCGGAATTAGAATCCGGTGCCTTAACGCTTGGCGATTCCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGa  <  2:536655/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  GGCGTGTTGAAACGTCAGGCGCTGAAATATGCGTTAGGGTTGTAGCCTTACAACATTGCCGGATGCGTGCCAACCGTAGGTCGGATAAGACGCGGCAGCGTCGCATCCGACATTGAAGGATAAGACGTGTCAACATCGCATTCGACATTGAATGAACGCAGAAAAGCAAAAAGCTCGCCGAAGCGAGCTTTTTTAATGTGGCTGGGGTACGAGGATTCG‑AACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACTGGGTGCACTTACAAGGTAAGCGTCTTGAATAAATTGGCTGGGGTACGAGGATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAACAAATTGGTTTTGAATTTGCCGAACATATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGA  >  NC_000913/696232‑696711

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