Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,760,996	T→C	N85S (AAC→AGC)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,996	0	T	C	98.4%	218.7 / NA	64	N85S (AAC→AGC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (40/23);  minor base A (0/1);  total (40/24)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGA  >  NC_000913/3760754‑3761231                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            tGGCAGTGGCCGCCAACATCCATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGCTGGTATGCAGCCCGGTTCCCGTCAGGCTGATCAcc                                                                                                                                                                                                                                     >  2:361793/1‑251 (MQ=255)  ggCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCg                                                                                                                                                                                                                                    >  1:32823/1‑251 (MQ=255)     aGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGAc                                                                                                                                                                                                                                 >  1:328811/1‑251 (MQ=255)      gTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCTTCTTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACg                                                                                                                                                                                                                                <  1:17208/251‑1 (MQ=255)        ggaagaaaagaacAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGCTCACCGGACGc                                                                                                                                                                                                                               >  1:350956/13‑250 (MQ=255)          ccgccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAgc                                                                                                                                                                                                                            <  2:126948/251‑1 (MQ=255)            gccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAgcgc                                                                                                                                                                                                                          >  1:56026/1‑251 (MQ=255)                 aaccaatacccccgccgccgCATTGTTATTGCCGATCCAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGAGGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGctct                                                                                                                                                                                                                     <  2:296062/240‑1 (MQ=255)                     aataacaaacccgccgcATTGTTATTGACGATACACGCATCTTCCGACCCCGTACTACGGCACAGCAGCTGCGCCAACGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGCCGAACGCATACCCTGCGCAAAGGCTTCCACCGCGGCTTCCGCCTGTAACGATCGCCCAAGGTTGGAATGCAGCTCGATTCCCGCCAGGCTAGTCACCATACGCACCGCGCTCTccgc                                                                                                                                                                                                                 >  2:297206/11‑251 (MQ=255)                        aaCACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGcttc                                                                                                                                                                                                              >  2:173768/1‑251 (MQ=255)                          cacCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGGATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGcttctt                                                                                                                                                                                                            >  2:44340/1‑251 (MQ=255)                            ccgccgccgCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTc                                                                                                                                                                                                          >  1:271797/1‑251 (MQ=255)                                            aTTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCCAAGCCGCTGCGCCAGCGCCCGATCGCGATGTCCGCGCCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTACAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCgc                                                                                                                                                                                          >  1:227136/1‑250 (MQ=255)                                            aTTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCg                                                                                                                                                                                           >  2:283416/1‑250 (MQ=255)                                                 cGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCAGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGActtct                                                                                                                                                                                     >  1:36333/1‑251 (MQ=255)                                                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTg                                                                                                                                                                                   <  2:28170/251‑1 (MQ=255)                                                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTg                                                                                                                                                                                   >  2:63998/1‑251 (MQ=255)                                                     aCAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTgc                                                                                                                                                                                  <  1:123971/250‑1 (MQ=255)                                                         gCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTACAGCTCGCCCAAGGTTGGTATGCAGCACGGGTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGccc                                                                                                                                                                              >  2:314795/1‑250 (MQ=255)                                                         gCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCCAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGCCTTCTTGCGccc                                                                                                                                                                              >  2:195005/1‑250 (MQ=255)                                                         gCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGcccc                                                                                                                                                                             >  2:207215/1‑250 (MQ=255)                                                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGCGGTCCGCGTCCGGCGTCGTCCACATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGAGGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCAcacc                                                                                                                                                                  <  2:251898/251‑1 (MQ=255)                                                                      ccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCAcaccac                                                                                                                                                                >  1:188722/1‑251 (MQ=255)                                                                         taataaGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATAATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACgc                                                                                                                                                              >  2:131915/1‑250 (MQ=255)                                                                          aaTACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACgc                                                                                                                                                              >  2:351675/1‑249 (MQ=255)                                                                          aaTACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACgc                                                                                                                                                              >  1:218025/1‑249 (MQ=255)                                                                             aCGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACgcaggc                                                                                                                                                          >  1:70229/1‑250 (MQ=255)                                                                                     gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGTCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAACGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCCGCGTCTgg                                                                                                                                                 >  2:150339/1‑251 (MQ=255)                                                                                      cagcTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGc                                                                                                                                                <  2:369979/251‑1 (MQ=255)                                                                                         cTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGc                                                                                                                                             >  1:53451/1‑251 (MQ=255)                                                                                           gcgcCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGccc                                                                                                                                           >  2:272725/1‑250 (MQ=255)                                                                                            cgcCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCAc                                                                                                                                          >  1:130636/1‑251 (MQ=255)                                                                                               cAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGa                                                                                                                                        >  2:284101/1‑250 (MQ=255)                                                                                                      cgagcgcgATGTCCTCGTCCGGCGTCGTCCAGATCATACTCGAGGTTCCCTGGCGAAGGTATAGCCTGCGCACCGGCTTCCCCGGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGTGCATCGACTTCTTGCGCCCAGTTTTCAGACCGCGCTGTCATCGGCTGGCTGCCGCGAATTACTtc                                                                                                                                <  2:27744/247‑1 (MQ=255)                                                                                                        aTCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGACTCACTtctc                                                                                                                              >  2:222685/1‑251 (MQ=255)                                                                                                           gcgATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTcgcg                                                                                                                           >  1:294577/1‑251 (MQ=255)                                                                                                             gATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCt                                                                                                                         >  1:80976/1‑251 (MQ=255)                                                                                                                   cgcgTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTtcgtcg                                                                                                                   >  1:335047/1‑251 (MQ=255)                                                                                                                     cgTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTtcgtcg                                                                                                                   >  1:269199/1‑249 (MQ=255)                                                                                                                                      aTCATACTCGAGGGTCACTGGCGACCGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCGCGCCCAAGGTGGGTATGCAGCACGGTTCCCTTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCGTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTTGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAAtt                                                                                                <  1:199986/251‑1 (MQ=255)                                                                                                                                       tCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTc                                                                                               >  1:234751/1‑251 (MQ=255)                                                                                                                                        cATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCGAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTcc                                                                                              >  2:189317/1‑251 (MQ=255)                                                                                                                                           aCTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTccacc                                                                                           <  1:64124/251‑1 (MQ=255)                                                                                                                                              cgtgggTCTCTGGCGAACGCAGAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCGGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCGTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCTTCGAGCATCTGACTCAACAATTccccccc                                                                                         <  2:210595/247‑6 (MQ=255)                                                                                                                                                  gggcACTGGCGAATGAATTGCCTGCTCACCGGCGTCCCCCGGGGTTTCCGCCTGGAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCGCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACGCCACGCAGGCAGCGTCTGGCTGCCACGAATCGCTTCTCGCGCGTCGTCGAGCATCTGACGCAACAATTCCACCACGCggg                                                                                    <  2:372991/248‑1 (MQ=255)                                                                                                                                                         ggCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGAcc                                                                             <  1:258963/251‑1 (MQ=255)                                                                                                                                                                    aGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACGCCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCa                                                                   <  1:357875/250‑1 (MQ=255)                                                                                                                                                                               cGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAaagg                                                       >  1:86936/1‑251 (MQ=255)                                                                                                                                                                                    tCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAaaggaagg                                                   <  1:100175/250‑1 (MQ=255)                                                                                                                                                                                     cccccGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAATAAAGGAAGGAg                                                 <  2:334905/248‑1 (MQ=255)                                                                                                                                                                                      cACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGc                                                >  1:108060/1‑251 (MQ=255)                                                                                                                                                                                          gcgGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATc                                            <  1:284101/251‑1 (MQ=255)                                                                                                                                                                                           cgGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCGGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTAGCACGCAAAGAAAGGAAGGAGCTATcg                                           <  1:288782/251‑1 (MQ=255)                                                                                                                                                                                                   gCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATaa                                   <  1:222685/251‑1 (MQ=255)                                                                                                                                                                                                     cTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGc                                 >  1:335083/1‑251 (MQ=255)                                                                                                                                                                                                               cGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATa                         <  1:97097/249‑1 (MQ=255)                                                                                                                                                                                                                gCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGATGATCACCGGACGCATCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCATCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGCAAGGAAGGAGCGATCGCGCAATAAGCGATCAATAGcc                      <  2:161196/251‑1 (MQ=255)                                                                                                                                                                                                                 cccAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCtt                                                                                                                        <  1:50988/152‑1 (MQ=255)                                                                                                                                                                                                                 cccAAGGTTGGTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCtt                                                                                                                        >  2:50988/1‑152 (MQ=255)                                                                                                                                                                                                                          ggTATGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGAc            >  1:344056/1‑251 (MQ=255)                                                                                                                                                                                                                             aTGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTg                                                                                                           <  2:351248/153‑1 (MQ=255)                                                                                                                                                                                                                             aTGCAGCACGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTg                                                                                                           >  1:351248/1‑153 (MQ=255)                                                                                                                                                                                                                                    aCGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAggaa  <  1:207215/251‑3 (MQ=255)                                                                                                                                                                                                                                    aCGGTTCCCGTCAGGCTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAgga   >  2:287650/1‑249 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            TGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAAGCTCGCCCAAGGTTGGTATGCAGCACGGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGA  >  NC_000913/3760754‑3761231

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