Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,761,230	G→A	S7F (TCC→TTC)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,761,230	0	G	A	84.9%	123.9 / 11.3	54	S7F (TCC→TTC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base A (25/20);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

CGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGC  >  NC_000913/3760989‑3761475                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       cGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGttt                                                                                                                                                                                                                                               <  1:406974/251‑1 (MQ=255)        ttGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGt                                                                                                                                                                                                                                        >  1:382317/1‑251 (MQ=255)                     gCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCAACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTg                                                                                                                                                                                                                           >  2:343510/1‑251 (MQ=255)                      cAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTgg                                                                                                                                                                                                                          >  1:422549/1‑251 (MQ=255)                                      tcttTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGTT‑CAAGGCCctct                                                                                                                                                                                                         >  2:126865/1‑251 (MQ=255)                                               aaCCGGGCATCGACTTCTTGCGCCCAGTTTTCCCACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGgag                                                                                                                                                                                                <  2:29182/251‑1 (MQ=255)                                               aaCCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGgag                                                                                                                                                                                                <  1:219009/251‑1 (MQ=255)                                                 ccGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGaa                                                                                                                                                                                              <  1:383246/251‑1 (MQ=255)                                                  cGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGCCGCACCAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAACGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATCGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGgaga                                                                                                                                                                                               >  2:104083/1‑249 (MQ=255)                                                     gCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGc                                                                                                                                                                                          >  1:184829/1‑251 (MQ=255)                                                        tCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGtt                                                                                                                                                                                       >  1:73204/1‑251 (MQ=255)                                                                   cgcCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCa                                                                                                                                                                            <  1:108336/251‑1 (MQ=255)                                                                    gcCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGACAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCAc                                                                                                                                                                           >  1:245425/1‑251 (MQ=255)                                                                       cAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGATAAGGGCGTTTAACATAACCACgg                                                                                                                                                                         >  1:1739/1‑250 (MQ=255)                                                                              cacaCCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGAGCTCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAAc                                                                                                                                                                 >  2:400044/1‑251 (MQ=255)                                                                              cacaCCACGCAGGCAGCGTCTGGCTGCCACGAATCACCTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAAc                                                                                                                                                                 <  1:374395/251‑1 (MQ=255)                                                                              cacaCCACGCAGGCAGCGTCGGGCTGCCACGAATCCCTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGGTGACTATAGAGG‑AACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCATGGATTGTAACg                                                                                                                                                                <  2:284648/251‑1 (MQ=255)                                                                                                                    tctcGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGcc                                                                                                                            <  2:427651/250‑1 (MQ=255)                                                                                                                       cgcgCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCttt                                                                                                                         <  1:280618/250‑1 (MQ=255)                                                                                                                             tcgtcgAGCATCTGACGCAACAAATCCCCCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCGATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGt                                                                                                                   <  2:73204/250‑1 (MQ=255)                                                                                                                              cgtcgAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTc                                                                                                                 >  1:310821/1‑251 (MQ=255)                                                                                                                               gtcgAGCATCTGACGCAACAATTCCCCCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCg                                                                                                                <  2:451655/251‑1 (MQ=255)                                                                                                                               gtcgAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATCGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGAGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCg                                                                                                                >  2:90110/1‑251 (MQ=255)                                                                                                                                  gAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTg                                                                                                             >  1:55368/1‑251 (MQ=255)                                                                                                                                          gACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAg                                                                                                         <  2:379640/247‑1 (MQ=255)                                                                                                                                            cgctacaattccccccccggggggtgACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGGTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGTCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGTACAGATCGCGCATCAAGCCTTTGGCGGTTCGGTTCGAGCAAAAc                                                                                                   <  2:442974/228‑1 (MQ=255)                                                                                                                                              caacaaTTCCACCACGCGGGTGTGACCATAATCATCACGCAAAGAAAGGAAGGAGAGATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACTCGTTTCGGTTTTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACTGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGAGGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCATAACtt                                                                                                 <  2:324013/251‑1 (MQ=255)                                                                                                                                               aacaaTTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAAACACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTAGGGGCGGGCAAAACttt                                                                                                >  2:80169/1‑251 (MQ=255)                                                                                                                                                  aaTTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACAt                                                                                                                                                                                 <  2:237449/167‑1 (MQ=255)                                                                                                                                                  aaTTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACAt                                                                                                                                                                                 >  1:237449/1‑167 (MQ=255)                                                                                                                                                                    tgACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTGACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCgg                                                                           >  1:170468/1‑250 (MQ=255)                                                                                                                                                                     gACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGATAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGa                                                                           <  2:1739/250‑1 (MQ=255)                                                                                                                                                                       ccATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGATCGGTGCGGGCAAAACTTTCGCGGCTAAACCGATTTTTGAcc                                                                         >  2:7735/1‑250 (MQ=255)                                                                                                                                                                            aGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGCTTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAAATTTCGCGGCTAAACAGGTTTTCGACCAGTTTg                                                                   >  1:287081/1‑251 (MQ=255)                                                                                                                                                                                    gCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAAAAGGTTTTCGACCAGTTTGACTAAgtg                                                           >  1:275435/1‑251 (MQ=255)                                                                                                                                                                                          aaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCCGTTTGACTAAGTGTGGGCg                                                     >  1:22893/1‑251 (MQ=255)                                                                                                                                                                                            aggaaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGt                                                   >  1:340587/1‑251 (MQ=255)                                                                                                                                                                                             ggaaggaGCAATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGTAATCGCGTTTCGGGTGCCATAGGACAGCTGGCGGATACAAT‑CCCTCTCACACGGAGAAGGGCATTTAACATAACCACGACTTGTAACATGAGCTAGGTTCGGAGGACATCTCCCGCATCACACCTTTGGCGTGTCGGTGACAGCAAGACTTTCACGGCGCAACAGGATTTCGAACAGTTTGACTatgagggggggcac                                                  >  2:254860/1‑238 (MQ=255)                                                                                                                                                                                                 ggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATCGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGAGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGCCATATCGCGCATCAAGCCTTTGGCGGTTCGGGGCGAGATAAGCTTTCGCGGATAACCAGGTTTTAGGACAGTTTGGCTAAGAGTGGGCgggctgc                                               >  2:260125/1‑245 (MQ=255)                                                                                                                                                                                                    gCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGGCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGAGTGGGCGGTCGaagcac                                           >  1:186176/1‑246 (MQ=255)                                                                                                                                                                                                      tATCGCGCAATAAGAGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGCGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGAAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGAATAAGGTTGGGCGGTCGACACAcc                                          >  2:261098/1‑250 (MQ=255)                                                                                                                                                                                                                   gCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAAcg                            <  1:120678/251‑1 (MQ=255)                                                                                                                                                                                                                      aTCAATAGTCGGAAGTTGACTAAAGAGGAAACGCGTTACGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGGTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACgcgg                         <  2:186176/251‑1 (MQ=255)                                                                                                                                                                                                                      aTCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTg                                                                                                             <  1:23921/167‑1 (MQ=255)                                                                                                                                                                                                                            aGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGaa                    <  1:110166/250‑1 (MQ=255)                                                                                                                                                                                                                            aGCCGGAAGTTGACTATAGAGGAAACGCGTTTAGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTAACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGAGAGCCAAACTTTCGCGGCTAAACAGGTTTTCGACCCGTTTGACTAAGTGTGAGCGGTCGACACCCCCACCTGGCGCACCGCGGCGGaa                    >  2:153128/1‑250 (MQ=255)                                                                                                                                                                                                                            aGCCGGAAGTTGACTATAGAAGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGTAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGaaa                   <  2:417268/251‑1 (MQ=255)                                                                                                                                                                                                                             gCCGGACGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCg                                                                                                                >  2:283109/1‑157 (MQ=255)                                                                                                                                                                                                                             gCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATGGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCg                                                                                                                <  1:283109/157‑1 (MQ=255)                                                                                                                                                                                                                              ccGGAAGTTGCCTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAAAAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGCCACACCAACCTGGCGCAACGCCGCGGGAAtt                 >  2:221103/1‑251 (MQ=255)                                                                                                                                                                                                                              ccGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGAAACACCAACCTGGCGCAACGCGGCGGAAAtt                 >  2:210349/1‑251 (MQ=255)                                                                                                                                                                                                                                 gAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGACGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGATCGGTGCGGGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGCAATTg                >  1:261890/1‑249 (MQ=255)                                                                                                                                                                                                                                             agagGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACgg   <  1:343510/250‑1 (MQ=255)                                                                                                                                                                                                                                             agagGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGc  <  1:427265/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       CGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAT‑CAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGC  >  NC_000913/3760989‑3761475

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