Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	64,749	G→A	H344H (CAC→CAT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	64,749	0	G	A	97.3%	112.1 / ‑6.0	37	H344H (CAC→CAT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/24);  minor base T (1/0);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAG  >  NC_000913/64504‑64989                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  taCGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCTAGGAAAAATGGCATGATTTTAGTTTTATggt                                                                                                                                                                                                                                               >  2:35623/2‑249 (MQ=255)                ttATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCa                                                                                                                                                                                                                                <  1:254165/250‑1 (MQ=255)                   taGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGGAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTc                                                                                                                                                                                                                            >  2:90426/1‑251 (MQ=255)                       ccagcaccagtaCCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCCCTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCa                                                                                                                                                                                                                         <  2:212021/250‑1 (MQ=255)                        cagcacaagtacCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCAGACAGGGCTGACGTGTGGCGGCACTTCGCCGAGCTTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGCTGACCAAACGCCGCCACCGATCACCCGTGGCGGTATACCGGCAGGCATGTCACAGTCGCTCGTTCGCGTTAATATGGCATGTTTTTAGTTTTTTGGAAAATCCACGTCCCCAGcgcacaa                                                                                                                                                                                                                        >  1:249418/1‑244 (MQ=255)                           caccagtacCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCtt                                                                                                                                                                                                                    <  1:39194/251‑1 (MQ=255)                              cagtaccgAATCATAAAGCCCTGGCCGTGAATCCAGCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCa                                                                                                                                                                                                                 <  1:26429/244‑1 (MQ=255)                                   ccgAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCTCCAGGGCTGTCGTGTGTCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGATAATACCAAACGCCGCCACTGATCCCCCTTGGCTGTCTACCGTCAGGCCGTTCACCGTCGCTCTTTCGAGTAAAAATGGCATGTTTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTATCAGGTta                                                                                                                                                                                                            <  1:172065/249‑1 (MQ=255)                                        tCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCAACAGCTCGCAATCTTTCAGGTTATAAGt                                                                                                                                                                                                       >  2:231923/1‑251 (MQ=255)                                         aagaaagaccgggCCGTGAATCCGTCACGTAGCCGCCAGAGCTGGCGTGTGGAGGCACTTCGCCGAGATTAGGCGCGACATAGCCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCAACGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCAACGCCGCTCGTTCGAGTAAAAAGGGCATGATTGCAGTTTTATGGAAAAGCTGCGTCACCAGCTCGCAATCCTTCAGGTTATAAGt                                                                                                                                                                                                       <  2:73763/239‑1 (MQ=255)                                              aGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAg                                                                                                                                                                                                 <  1:129805/251‑1 (MQ=255)                                                                         ccgccAGGGCTGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGATCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAAcggcg                                                                                                                                                                      >  2:120790/1‑251 (MQ=255)                                                                                  atGGCGTGTGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCCGTCACCGTCGCGCGTTCGAGGAAAAAAGGCATGATTTCAGGTTTATGGGAAATCTGCTTCACCAGCTCGCAATCTTTCAGGTTCTAAGTCGCCAGCGCCGGTGTATCTTCGGAGAAACGGCGGTACAAttt                                                                                                                                                             >  2:208897/2‑251 (MQ=255)                                                                                            gcggcACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGa                                                                                                                                                  <  1:231923/251‑1 (MQ=255)                                                                                                           gaTTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCg                                                                                                                                   <  1:203293/251‑1 (MQ=255)                                                                                                               aGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAg                                                                                                                               <  1:246336/251‑1 (MQ=255)                                                                                                                   gcgACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAAATTCGGCCATTCGATCCCACGGGGTTTCGATaacatt                                                                                                                           >  2:220305/1‑246 (MQ=255)                                                                                                                     gACATAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTc                                                                                                                         <  2:52046/251‑1 (MQ=255)                                                                                                                      aCATAACCAGCGCGGTGCATACGAGGGAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCGCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTGTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTGCGTCCATTCGATCCCACGGGTTATCGATAGATTTTcc                                                                                                                        <  2:94827/251‑1 (MQ=255)                                                                                                                      aCATAACAAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGACAATCTGCGTCACCAGCTCGCACTCTTTCAGGTTATCAGTAGCCCGCGCAGGTTTATCTTCGGCGAAACGGCGGG‑CCATTTCGTCCATTTTATCCCACGGGGTATCTATAGATTTTcc                                                                                                                        >  2:53938/1‑251 (MQ=255)                                                                                                                         tAACCAGCGCGGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCCTTCGATCCCACGGGTTATCGATAGATTTTCCTTc                                                                                                                     >  1:151641/1‑251 (MQ=255)                                                                                                                                cgcgGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAAt                                                                                                              <  1:214615/251‑1 (MQ=255)                                                                                                                                cgcgGTGCATACGAGGAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGGTATCGATAGATTTTCCTTCTCCTAAt                                                                                                              >  1:82465/1‑251 (MQ=255)                                                                                                                                            gaggAAAATAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑TAATTTCGTCCATTCGATCCCACGGGTTATCGATAGAATTTCCTTCGCCTACTAACTCCTGAGCg                                                                                                  >  2:90344/3‑251 (MQ=255)                                                                                                                                                   aTAGAGATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGAcgggtt                                                                                           >  1:226213/1‑246 (MQ=255)                                                                                                                                                     agagATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTc                                                                                          <  1:148857/250‑1 (MQ=255)                                                                                                                                                      gagaTGACCAAACGCCGACACCGATTCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTAGATCGGTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCGTCGGCGAAACGGCGGT‑CGATTGCGTCCATTCGATCCCACGTGTGATCGATAGATTTTCCTTCGCCTAATAACGCCTGAGCgcagggttcta                                                                                        <  2:178985/251‑11 (MQ=255)                                                                                                                                                                         aCCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAAtt                                                                     <  1:166707/251‑1 (MQ=255)                                                                                                                                                                          ccgatcccccgTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTc                                                                    <  1:70926/243‑1 (MQ=255)                                                                                                                                                                                                                          cGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACTGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCATAATAACTCCTTAGCGTCGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGTATTTCAGCGCCTCGATACCTTCGATAATTAGCCGAcc                     <  1:180493/250‑1 (MQ=255)                                                                                                                                                                                                                          cGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGAcc                     <  1:180482/250‑1 (MQ=255)                                                                                                                                                                                                                               aaaaaTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGc               >  2:242759/1‑251 (MQ=255)                                                                                                                                                                                                                                aaaaTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGcc              <  1:236712/251‑1 (MQ=255)                                                                                                                                                                                                                                  aaTGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTg            <  2:209172/251‑1 (MQ=255)                                                                                                                                                                                                                                       cATGATTTCAGTTGTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTGAGCCTGGGCaa       <  2:25253/251‑1 (MQ=255)                                                                                                                                                                                                                                          gATTTCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCaaaaa    <  1:242759/251‑1 (MQ=255)                                                                                                                                                                                                                                             ttCAGTTTTATGGAAAATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAg  <  1:26468/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  GACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGT‑CAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAAG  >  NC_000913/64504‑64989

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