Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	64,665	T→C	A372A (GCA→GCG)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	64,665	0	T	C	97.7%	134.5 / NA	43	A372A (GCA→GCG)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (14/28);  minor base A (1/0);  total (15/28)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAAT  >  NC_000913/64419‑64911                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       aaaaCCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACgccg                                                                                                                                                                                                                                                    <  1:124809/251‑1 (MQ=255)          ggTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGatc                                                                                                                                                                                                                                           <  1:16988/251‑3 (MQ=255)                gtgCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGCTGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAAAGCCGCCACCGATCCCCCgt                                                                                                                                                                                                                                     >  2:188806/1‑251 (MQ=255)                gtgCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGatcccccgg                                                                                                                                                                                                                                     >  2:324313/1‑243 (MQ=255)                gtgCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCgt                                                                                                                                                                                                                                     <  1:281431/251‑1 (MQ=255)                            aGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCg                                                                                                                                                                                                                         >  1:91273/1‑251 (MQ=255)                                       tGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTc                                                                                                                                                                                                              >  2:56025/1‑251 (MQ=255)                                        gCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTCCCGGCAGGCCGTTct                                                                                                                                                                                                             >  2:2661/1‑250 (MQ=255)                                          cTTCCACCAGACCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGAACCAGTACCGAATCATAAAGCACTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACAAAACGCCGCCACCGATCCCCCGTGGCGGTATACCAGCATGCCCTTCAcc                                                                                                                                                                                                           >  2:301592/1‑251 (MQ=255)                                             ccaccaGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGtc                                                                                                                                                                                                        <  2:121029/251‑2 (MQ=255)                                              caccaGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGCGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTAACGGCAGGCCGTTCACCGtcc                                                                                                                                                                                                       >  2:221298/1‑249 (MQ=255)                                                ccaGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGtcgct                                                                                                                                                                                                     <  1:150772/251‑5 (MQ=255)                                                     ccGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCTGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGTGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGTCGGTCTACCGGCATGCCGTTAACCGTCGCTCGTTc                                                                                                                                                                                                <  2:156553/251‑1 (MQ=255)                                                         cGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGt                                                                                                                                                                                            <  2:57103/251‑1 (MQ=255)                                                           ggATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCCCGTTCGAGTa                                                                                                                                                                                           >  1:312932/1‑250 (MQ=255)                                                                 atcaGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATgg                                                                                                                                                                                    <  2:67851/251‑1 (MQ=255)                                                                      aaaGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGa                                                                                                                                                                               <  1:38525/251‑1 (MQ=255)                                                                               gatgatCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGttt                                                                                                                                                                      <  1:241729/251‑1 (MQ=255)                                                                                        cGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGaaaa                                                                                                                                                             <  1:324313/251‑3 (MQ=255)                                                                                        cGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGaaa                                                                                                                                                              >  1:357006/1‑249 (MQ=255)                                                                                         gggTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAGTCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGaaaat                                                                                                                                                            >  1:217919/1‑248 (MQ=255)                                                                                                       ataGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCTTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGc                                                                                                                                               <  1:17064/250‑1 (MQ=255)                                                                                                         aGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTc                                                                                                                                             <  1:195345/250‑1 (MQ=255)                                                                                                               gcaccagtacCGAATCATAAATCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCt                                                                                                                                      <  1:186065/251‑1 (MQ=255)                                                                                                                                              ggAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAGCGCCGCCGCCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCTCTCGTTCGAGTAAAAATTGCATGATTGCAGTTTTATGGAAAATCTGCGTCACCAGCTCGTAGTCTTTCAGGTTAGAAGTTGCCAGCGCAGGTTTAt                                                                                                       <  2:164920/250‑1 (MQ=255)                                                                                                                                                          gTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGGTATAAGTCGACAGCGAAGGTTTATCTTCGGCGAAAc                                                                                           <  2:263240/251‑1 (MQ=255)                                                                                                                                                                     ggCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAAtt                                                                                <  2:222898/251‑1 (MQ=255)                                                                                                                                                                         ggCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGt                                                                            <  1:56025/251‑1 (MQ=255)                                                                                                                                                                                     cACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCCGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCCTTCGATccc                                                                >  2:15857/1‑251 (MQ=255)                                                                                                                                                                                                 aTTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGa                                                    <  2:91273/251‑1 (MQ=255)                                                                                                                                                                                                       cgcgACATAGCCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCATTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCTATCCCACGGGTTATCGATAGAtt                                              <  2:113841/251‑1 (MQ=255)                                                                                                                                                                                                                  cAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCTCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAATGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCt                                   <  1:207536/251‑1 (MQ=255)                                                                                                                                                                                                                       cggtGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATa                               >  1:135159/4‑250 (MQ=255)                                                                                                                                                                                                                         gtGCATACGAGGAAAATAGAGATGACCAAACGCCGCCGCCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGGTATAAGTCGCCAGCGCAGGGTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGAGCCCACGGGTTATCGATGGATTTTCCATGGCCTAATaact                            <  2:239809/250‑4 (MQ=255)                                                                                                                                                                                                                          ggcatacgaggAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCGCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTc                           <  2:135159/242‑1 (MQ=255)                                                                                                                                                                                                                                  aggAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCg                    >  1:126430/2‑250 (MQ=255)                                                                                                                                                                                                                                  aggAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGGCCATTTCGTCCATTCGATCCCCAGGGTTATCGATAGATTTTTCTTCGCCTAATaaatcctgtgcg                    >  2:197687/2‑239 (MQ=255)                                                                                                                                                                                                                                       aaTAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGtt              <  2:126430/251‑1 (MQ=255)                                                                                                                                                                                                                                          agagATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTc            <  1:46837/250‑1 (MQ=255)                                                                                                                                                                                                                                           gagaTGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTGCGCCTAATAACTCCTGAGCGACGGTTTCCa          <  2:337307/251‑1 (MQ=255)                                                                                                                                                                                                                                            agaTGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAg         <  1:2767/251‑1 (MQ=255)                                                                                                                                                                                                                                             gaTGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGc        <  1:197687/251‑1 (MQ=255)                                                                                                                                                                                                                                                   cAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAAt  >  1:307096/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       AAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAAT  >  NC_000913/64419‑64911

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: ‑