Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,760,132	(G)6→7	coding (1118/1392 nt)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,132	1	.	G	89.7%	115.3 / 7.8	39	coding (1118/1392 nt)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (16/19);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TTGCTCATCTTCAAGGCAGCGTAAATCCAGCCACAATCGTCCGTCATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑G‑CCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTA  >  NC_000913/3759885‑3760360                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                     ttGCTCATCTTCAAGGCAGCGTAAATCCAGCCACAATAGTCCGTCATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGAATCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCCGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAg‑g‑ggg‑ggcc                                                                                                                                                                                                                                    >  2:53807/1‑248 (MQ=255)  tGCTCATCTTCAAGGCAGCGTAAATCCAGCCACAATCGTCCGTCATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAg‑g‑ggg‑ggcc                                                                                                                                                                                                                                    <  2:346511/250‑4 (MQ=255)   gCTCATCTTCAAGGCAGCGTAAATCCAGCCACAATCGTCCGTCATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAg‑g‑ggg‑ggcctg                                                                                                                                                                                                                                  >  2:313043/1‑246 (MQ=255)   gCTCATCTTCAAGGCAGCGTAAATCCAGCCACAATCGTCCGTCATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAg‑g‑ggg‑ggcct                                                                                                                                                                                                                                   <  1:327914/250‑5 (MQ=255)                          tCCAGCCACAATCGTCCGTCATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGa                                                                                                                                                                                                            >  1:70475/1‑250 (MQ=255)                                 acaATCGTCCGTCATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑G‑CCTGTAAACGTTGTGCCTGGATTTGAATGACctct                                                                                                                                                                                                   >  1:320521/1‑251 (MQ=255)                                      cgtccgtcATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCCTGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCCGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAACGTGCGGCAAG‑G‑GGG‑GGCCTGTAACCGTTGTGCCTGGATTTGAATGACCTctgcgc                                                                                                                                                                                               >  2:69569/1‑251 (MQ=255)                                       gtccgtcATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcgct                                                                                                                                                                                              >  1:212679/1‑251 (MQ=255)                                       gtccgtcATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcgct                                                                                                                                                                                              <  2:351890/251‑1 (MQ=255)                                       gtccgtcATAAATACGACCAACCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGGGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑GTGGG‑G‑CCTGTAAACGTTGTGCCTGGATTTGAATGACCTctacgct                                                                                                                                                                                              <  2:263708/251‑1 (MQ=255)                                                    aCGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGCTTTGAATGACCTCTGCGCTGCGGGTAAGCAgg                                                                                                                                                                                 >  2:76921/1‑251 (MQ=255)                                                            tCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGt                                                                                                                                                                         >  1:88268/1‑251 (MQ=255)                                                                   cACTGGCACTTCCCGCCGGCGGTCGGCTAATGACTCAAGCTGGCTACCGCTTCAATCATGGGTTGTAATCGTGACTGCCGCGCTCTGCATGCGATCAACGGGCAGTGAACCACGGCCAATCTGCGAAAGACATGGAAAAACCTGTACCTAAAACTCCGCGCCGTAATGTGCGGCAAGTG‑GGG‑A‑CCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGTTAAGCAGGCGCGGCGTCGGTAAt                                                                                                                                                                  <  2:331702/251‑1 (MQ=255)                                                                    ccTGGCAAATCACGCCAGAGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCGGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAAGGTGCGGCAAG‑G‑GGGTG‑CCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGTGTAAGCAGGCGCAGGGTCGGTAAtt                                                                                                                                                                 <  2:116748/250‑1 (MQ=255)                                                                         cAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTc                                                                                                                                                             >  1:101163/1‑250 (MQ=255)                                                                               ccGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCaga                                                                                                                                                      <  1:338570/250‑1 (MQ=255)                                                                                                      cAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAc                                                                                                                                <  1:153408/250‑1 (MQ=255)                                                                                                            ggCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCGCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCGTCAGGGTGTAAATAAAGACGCAACGt                                                                                                                         <  2:228070/251‑1 (MQ=255)                                                                                                                     gTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTcc                                                                                                                  >  1:375354/1‑249 (MQ=255)                                                                                                                     gTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACAACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAACCTCCGCGCCGTAATGAGTGGCAAG‑G‑GGG‑GGCCTGTAAACGTGGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGGAGGGTCGGTAATTTTTc                                                                                                                                                             >  2:355154/1‑206 (MQ=255)                                                                                                                            atGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCcgcg                                                                                                         <  1:176428/251‑1 (MQ=255)                                                                                                                             tGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCcgcgg                                                                                                        >  2:126144/1‑250 (MQ=255)                                                                                                                                     aaaCGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCAt                                                                                                 <  1:276414/250‑1 (MQ=255)                                                                                                                                            aaTGCCGCGCTTGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCAGCGCCGCAATGTGCGGCAAG‑G‑GGG‑GGACTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGGAATTTTTCACTCAGAGCTTCTGGGTGTCAATACAGACGCAACGTGGATTCCAGCGCCGCGAGGGGCAGTTTATCcg                                                                                         >  2:193544/1‑251 (MQ=255)                                                                                                                                              tGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCAtt                                                                                                >  1:304819/1‑242 (MQ=255)                                                                                                                                              tGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCAtt                                                                                                <  2:304819/242‑1 (MQ=255)                                                                                                                                               gCCGCGCTCGGCAGGCTATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCcgcgc                                                                                      <  1:358233/251‑1 (MQ=255)                                                                                                                                               gCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCcgcgc                                                                                      >  1:330113/1‑251 (MQ=255)                                                                                                                                                    gcTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTaa                                                                                  <  2:326705/250‑1 (MQ=255)                                                                                                                                                    gcTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTaa                                                                                  >  1:326705/1‑250 (MQ=255)                                                                                                                                                       cGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCa                                                                              <  2:93059/251‑1 (MQ=255)                                                                                                                                                                                  aatgcaaaTCTGCGAAGAACATGGCAGAACCTGTCCCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTGCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTTGCTTTGCAggcgg                                                   <  2:126236/245‑1 (MQ=255)                                                                                                                                                                                     gCCAAGCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGTTCATTTTATCCGCGCGTAATGCACGCTTCATCGGGTGGCTTTGCAggcggg                                                  <  2:375354/249‑1 (MQ=255)                                                                                                                                                                                      ccAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGa                                               <  1:44004/251‑1 (MQ=255)                                                                                                                                                                                              cGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAggcgggcggacctctct                                       >  2:111461/1‑242 (MQ=255)                                                                                                                                                                                                    aCATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCttttttt                                 >  1:350494/1‑250 (MQ=255)                                                                                                                                                                                                        ggCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACca                              <  1:4062/250‑1 (MQ=255)                                                                                                                                                                                                         gCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGAGCATCTCTTTTTTACcaac                            <  1:368790/251‑1 (MQ=255)                                                                                                                                                                                                             aaCCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACaata                        <  1:313043/251‑1 (MQ=255)                                                                                                                                                                                                             aaCCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACaata                        <  2:101163/251‑1 (MQ=255)                                                                                                                                                                                                              aCCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACaataa                       <  2:116789/251‑1 (MQ=255)                                                                                                                                                                                                                               cgcgCCGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCGTTTTTACCAACAATAATTCCTGCCTGCGGCCCg      <  1:310022/251‑1 (MQ=255)                                                                                                                                                                                                                                    cGTAATGTGCGGCAAG‑G‑GGG‑GGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTa  >  1:190683/1‑250 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                     TTGCTCATCTTCAAGGCAGCGTAAATCCAGCCACAATCGTCCGTCATAAATACGACCAATCACCGGCACTGGCAATTCACGCCAGCGGGCGGCTAATGACTCAAGGTGGCTACCGCGTCCATCATGGGGTGTAAACGTTAATGCCGCGCTCGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAG‑G‑GGG‑G‑CCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTA  >  NC_000913/3759885‑3760360

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