Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	64,883	G→A	L300L (CTA→TTA)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	64,883	0	G	A	92.5%	110.2 / ‑0.7	40	L300L (CTA→TTA)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (21/16);  total (23/17)
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.84e-01

GCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCA  >  NC_000913/64643‑65125                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   gaggAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCg                                                                                                                                                                                                                                          <  2:276758/249‑1 (MQ=255)                                     cgccgTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTAGCGTAAAAATGGCATGATTTCAGTTTTGTGGAGGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGGCCATTCGGTCCCACGGGGTATCGATTGATTTTCCTTCGCCTAATAGCTCCTGCGGGACAGGTTCCAGCGAGAaagaaggataactccagaa                                                                                                                                                                                                       >  2:43893/1‑237 (MQ=255)                                                     cGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGCCAATTTCGTCCCTTCGATCCCCCGGGTGTTCGATAAATTTTCCTTTGCCTAATAACTCCTGGGCGACGGGTTCCCGAGAGAAAGAAGAGAAAATCCAAAAcgagggttttagacccc                                                                                                                                                                                      >  2:246427/1‑236 (MQ=255)                                                                    cgtcgctcGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATaa                                                                                                                                                                       <  1:41186/244‑1 (MQ=255)                                                                               cGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCGCGCAATCTTTCAGGGTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCt                                                                                                                                                            <  2:68618/251‑1 (MQ=255)                                                                                 aGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTTAGGTTATAAGGCGCCAGCGCATGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCttt                                                                                                                                                          <  1:268797/251‑1 (MQ=255)                                                                                     aaaaTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAAATCCTGAGCGACGGTTTCCAGCGCGAATGAAGAGAAATTCCAGAACGCGGATTTAAGCGCCTCGCAACCGTCACTAATTAGCCGCCCTTTaacc                                                                                                                                                      >  1:330391/1‑248 (MQ=255)                                                                                      aaaTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTc                                                                                                                                                                                             <  1:187497/211‑1 (MQ=255)                                                                                        aTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTgg                                                                                                                                                   <  1:341097/251‑1 (MQ=255)                                                                                         tGGCATGGTTTCAGATTTATGGCAAATCTGCGTCACCAGCTCGCAATATTTCAGTTTATAAGTCGCCAGCGAAGGTTTATCTTCGGCGAAACGGCGGTCAATTTAGTCCCTTCGATCCCACGGGTTCTCGATCGATTTTCCTTTGCCTAATTACACCTGTGCGCCGGTTTCCAGGGAGACTGTAGATAAAATCCAGGAGGCGGATTTCAGCGCCACGATACCGTTCATGATTACACGACCTTTATCCtagg                                                                                                                                                  >  2:373119/1‑248 (MQ=255)                                                                                          ggCATGATTTAAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAACCGGCGGTCACTTTCGTCCATTCGCTCCCACGGGGTATCGGAAGATTTTCCTTCGCCTAATAACTCCTGCGCGACGGTTTCCCGCGAGAATGAAGAGAAATTCCTGAACGCGGAGTTCATCGCATCGAGACCGTCGATAATTAACCGTCCtttcgcatgcgc                                                                                                                                                 >  2:81881/1‑242 (MQ=255)                                                                                               gATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTAGATAATTAGCCGACCTTTAGCCTGGGCaaaa                                                                                                                                             <  1:45512/250‑1 (MQ=255)                                                                                                       ttttATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGATAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGTGCAAAAAAGACGcc                                                                                                                                     <  1:137962/250‑1 (MQ=255)                                                                                                        tttatgGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGCCGGTTTCCAGCGCGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAACAAGACGCCGt                                                                                                                                   >  1:124193/5‑251 (MQ=255)                                                                                                                       cGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGACACGGCGGGCAATTTACTCCACTCGATCCCACGGGTTATCGTAAGATTTTCCTTCGCCTCATAACTCCTGAGCGCCGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTTAGCGCCCCGTTACCGTCGATAATTCTCCGCCCTTTATCCTGGGCCGACAAGACGCCGTTTTTAAAGCCGTGcc                                                                                                                    >  2:315894/1‑250 (MQ=255)                                                                                                                                 tCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCCTTCGATCCCACGGGGTATCGTTAGATTTTCCTTCGCCTAATAACTCCTGAGCGAAGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGCTTTCAGAGCCCCGATACCGTCGATAATTAACCGACCCTTAGCCTGGGCAAAAAAGACTCCGTTTTTAAAACCCTGCTCGCGCCAccc                                                                                                          >  2:32868/1‑249 (MQ=255)                                                                                                                                                              cgcAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAAAGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGccccagccgc                                                                              >  1:386715/1‑243 (MQ=255)                                                                                                                                                               gcAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAg                                                                            <  1:43893/251‑1 (MQ=255)                                                                                                                                                                     ttATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCCCACTGTTATCGCGCCCCAGCCGCAGCGGaat                                                                      >  1:64654/1‑250 (MQ=255)                                                                                                                                                                           tCGGCGAAACGGCGGTCAATTTCGTCCATTCGAACCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTGTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGt                                                                 <  1:315894/250‑1 (MQ=255)                                                                                                                                                                            cGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCCCCCCCGCCGCAGCGGAATACGGTa                                                                >  1:229418/1‑250 (MQ=255)                                                                                                                                                                               cGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTc                                                                                                                                                                            <  1:364612/139‑1 (MQ=255)                                                                                                                                                                               cGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGGTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTc                                                                                                                                                                            >  2:364612/1‑139 (MQ=255)                                                                                                                                                                                gAAACGGCGGTCAATTTCCTCCCTTCGTTCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATTGAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAACCGACCTTTAGCCTGGGCAAAAAAAACGCCGTTTTTAAAGCCGTGCTCGCGCCCCTCCCGCTCACTGTTATCGCGCCCCCGCCGCCGCGGAATACGGTAACGcc                                                           >  2:58536/1‑250 (MQ=255)                                                                                                                                                                                 aaaCGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTc                                                          <  1:248861/251‑1 (MQ=255)                                                                                                                                                                                     ggcggTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGc                                                       >  1:384918/1‑250 (MQ=255)                                                                                                                                                                                            aaTTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGttttt                                               >  2:375269/1‑251 (MQ=255)                                                                                                                                                                                              tttCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATCCGGTAACGCTCGGCATGTTTTTgc                                             >  1:167010/1‑251 (MQ=255)                                                                                                                                                                                               ttCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTgc                                             <  1:159008/250‑1 (MQ=255)                                                                                                                                                                                                         cGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGCCCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGc                                    >  1:68758/1‑249 (MQ=255)                                                                                                                                                                                                                aCGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCCGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAAc                           >  1:224245/1‑251 (MQ=255)                                                                                                                                                                                                                 cGGGTTATCGCTAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGCAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTACCGCTCGGCCTGTTTTTgaagaat                                        >  2:436815/1‑231 (MQ=255)                                                                                                                                                                                                                 cGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTCCCAGCGAGAATGAAGAGAAATTCCAGAACGCTGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCAt                                        <  1:436815/237‑1 (MQ=255)                                                                                                                                                                                                                  agattATAGATAGATTTTCCTTCGCATAATAACTCCGGGGCGACGGTTTACAGAGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGGGCCTCGAGACCGTCGATAATTAGACGACCTTGAGACTGGGCAACAAAGACGGCGTTTTTAAAGCCGTGCTCGCGCTACTCCAGCACACTGATATCGCGCACCAGCCGCAGCGGAATCCGGTAACGCTAGGCATGTTATTGCAGACTTCGCAGCTAGAACtc                         >  2:276648/4‑250 (MQ=255)                                                                                                                                                                                                                     ttATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGcac                      >  2:86426/1‑251 (MQ=255)                                                                                                                                                                                                                       aTCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTGGCAGCATGCGCAGATCGAACGGcacc                     <  1:129183/250‑1 (MQ=255)                                                                                                                                                                                                                                      ttCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCAAACTGCACCACGTTCCACCCgatg      >  2:31334/1‑250 (MQ=255)                                                                                                                                                                                                                                       tCGCCTAATAACTCCTGAGCTACTGTTGCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATCCCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACTCCGTTTTTAAAGCCGTGCTCGCGCCACTCTAGCTCACTTTTAGCGCGCCCCAGCCGCAGCGGAATACGGTAATGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCgatga     <  1:373119/250‑1 (MQ=255)                                                                                                                                                                                                                                         gCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATcc  >  1:27931/1‑250 (MQ=255)                                                                                                                                                                                                                                         gCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATca  <  1:246427/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   GCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCA  >  NC_000913/64643‑65125

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