Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,316	T→A	I155I (ATA→ATT)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,316	0	T	A	97.1%	93.3 / NA	36	I155I (ATA→ATT)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (11/22);  minor base C (1/0);  total (12/24)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCG  >  NC_000913/65071‑65558                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   gCATGTTGTTGCAGCATTCGCAGAACGAACTGCACCGCGATCCAACCGATTATCACGTCAGGAGCGTTGTTGGCAAACCCAGCGTTGAGTGTTTCCAGTAGCTGCGTGCGGCTGGCGACTTAGTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGTCCCAGCATAGAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGTTGGTGGTTTCAAGATCaata                                                                                                                                                                                                                                                 <  2:330391/249‑1 (MQ=255)   aTGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGaaa                                                                                                                                                                                                                                             <  1:449900/251‑1 (MQ=255)              gcaTTCGCAGATCGCACTGCACCACGTTCCAACCGATCATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGCTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCAGGGGGGTTCGATATCAATAAAAACCCACTTggac                                                                                                                                                                                                                                  >  1:149321/1‑248 (MQ=255)                  tCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTTACTATCACTAGAAAcaccctttcgcggcg                                                                                                                                                                                                                              >  2:423131/1‑237 (MQ=255)                                    cacGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCg                                                                                                                                                                                                            <  1:297404/251‑1 (MQ=255)                                     ccGTTCCAACCGATGATCCAATCAGGATCGTAGTTGGAAAACCAGGCGTTGAGTTTTTCCAGTAGCGGCGGGCGGCTTGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCGCACCGTGGCGGGTGGTTTCAATAGCAATATAAACCCACTTGAGCGGCGGACGATAGTCGTGATGCgg                                                                                                                                                                                                           <  2:242779/250‑1 (MQ=255)                                          ccAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTc                                                                                                                                                                                                       <  1:86426/250‑1 (MQ=255)                                               cggtgatCACCTCAGGTTCGTAGTTGGCAAACCAGGGGTGGAGTGTTTCCAGTAGCTGCGGGCGGTGGGCGATGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCTTTCCAGTCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACgg                                                                                                                                                                                                 <  2:362278/248‑1 (MQ=255)                                                       acacaaGGATCGTAGTTGGCAAACCAGGCGTTCAGTTTTTCCAGAAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAAc                                                                                                                                                                                         <  2:124193/246‑1 (MQ=255)                                                             ggATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCCATATCAATAGAAACCCCCTTGAGCGGCGGACGGTAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAgg                                                                                                                                                                                   >  1:107414/1‑250 (MQ=255)                                                              gATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAgcg                                                                                                                                                                                  <  1:395723/251‑3 (MQ=255)                                                                  gTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCAt                                                                                                                                                                              <  1:171719/251‑1 (MQ=255)                                                                       gtgCAAACCAGGGGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATAGCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGc                                                                                                                                                                         <  2:229418/249‑1 (MQ=255)                                                                          cAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGTCATTAACGATAGCGCCATTGTGCat                                                                                                                                                                       >  1:355551/1‑250 (MQ=255)                                                                           aaaCCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCatat                                                                                                                                                                     <  1:81881/251‑1 (MQ=255)                                                                           aaaCCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCatat                                                                                                                                                                     <  1:241127/251‑1 (MQ=255)                                                                           aaaCCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGGTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGAGGCGGGTTCCGACGGGACTTAAAGATAACGCAATTGTGCatat                                                                                                                                                                     >  1:64069/1‑251 (MQ=255)                                                                                aGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGGAGCGGTTTCCGACGGGCATTAACGATATCGCCATTGTGCctattcccc                                                                                                                                                                >  1:461315/1‑246 (MQ=255)                                                                                                cagtagcTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCAATATCAATAGAAAACCACTTGAGCGGCGGACGAAAGGCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCCTTGTGCCTATCACCCTCGAccccacccggtg                                                                                                                                                >  1:4662/7‑242 (MQ=255)                                                                                                cagtagcTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGGCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGTCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCggtg                                                                                                                                                <  2:355551/245‑1 (MQ=255)                                                                                                     gcTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAAGCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCggtgaggtg                                                                                                                                           <  1:283422/250‑1 (MQ=255)                                                                                                                              cattacaaaatcgcgcgcggtggCGTCGCCATTCTCTGGACCCAGCATATAAACGCTGCGCTGCCCGCAGCCTTCCAGGCCTATGCAGTCCAGCTCACCGTGGCGGGTGGTTGCAATAGCGATAGAAGCCCAGTTGAGCGGCGGACGATAGTCGGGATGCTGTTTCAGACGGGCATTAACGAGAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGtt                                                                                                                  <  2:132493/230‑1 (MQ=255)                                                                                                                                 ttcaaaaTCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGGGGTTTCAATATCAATAGAAACCCACCTTAGCGGCGGACGATAGGCGGGATGCGGTTTTAGACGGGCATTAAAGCAAGCGCCATTTTGCCTATCAACCTCTACCCACACCGGGGGGGGGAATAACCGCTCCCTCAGATAACGTTCTg                                                                                                               >  1:83784/5‑251 (MQ=255)                                                                                                                                         cgaGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCgggcg                                                                                                        <  1:32868/248‑4 (MQ=255)                                                                                                                                              gcggaggCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATc                                                                                                  <  1:92475/249‑1 (MQ=255)                                                                                                                                                  aggCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGcc                                                                                              <  1:366417/251‑1 (MQ=255)                                                                                                                                                        cgccttacgttggCCCCACAATATAAACGATGCGCTGCCCGCAGCCTTACAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCTCATCGGCCTCGTAg                                                                                        <  2:330722/240‑1 (MQ=255)                                                                                                                                                                              ataATCGATGCGCTGCCCGCAGCCGTCCAGGCCGAGGCAGTACAGCTCCCCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGca                                                                                                      <  2:117940/215‑1 (MQ=255)                                                                                                                                                                                 aaCGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCCCACCgggg                                                                                                                                                >  1:194902/1‑168 (MQ=255)                                                                                                                                                                                 aaCGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCggtg                                                                                                                                                <  2:194902/170‑1 (MQ=255)                                                                                                                                                                                        cgcTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTc                                                        <  2:345117/251‑1 (MQ=255)                                                                                                                                                                                           tGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGGGGACGATAGTCGGGCTGCGGTTTTAGGCGGGCATTAACGGTAGGGCCATTTTGCATATCACCCTCGACCCCCACCGGGGAGGTGATAAACCGATCCCTCAGATAGCGTTCTTGCGAGGGCACATCGGGCTCGCTGACGGGAAAGCCACCCTCAcaccggccgtttttgtt                                                     >  1:93031/1‑235 (MQ=255)                                                                                                                                                                                                       tCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAAt                                           <  1:375269/249‑1 (MQ=255)                                                                                                                                                                                                         cAGGCCGATGCAGTACAGCTCACCGTGGCGGGAGGTTTTAATATCCATAGAAACCCAATTGAGCGGCGGACGATACTGGGGATGAGGATTCACACAGGCATTTACGATAGCGCAATTGTGCATATCAACCTCGACCCACCCCGGGGAGGGGATAAAACGCCCCATCAGAAAACGTTCTGGCGGGCGCACATCGGCCCCGAAAACCGGAAAGCCCCACTCACGCCAGCGCTTTTCGTAATTTCTCAAtttgg                                       >  1:30244/1‑248 (MQ=255)                                                                                                                                                                                                                          gCTCACCGTGGCGGGTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTaa                      <  1:146973/251‑1 (MQ=255)                                                                                                                                                                                                                                  tGGCGGGTGGTTTCAATCTCAATAGAACCCCACTTGAACGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCCGCCTCGGAGACGGGAAAGCCACCTTCAAGCAGGCGGTTTTCGTAATTCCTCCATTGGCGACGGGCGCGGCAGTAAAAGccca               >  1:397458/1‑248 (MQ=255)                                                                                                                                                                                                                                      gggTGGTTTCAATATCAATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACAcc          <  1:361300/251‑1 (MQ=255)                                                                                                                                                                                                                                              tcaatatcaataGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGCGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGGCGCACATCGGCCTCGTAGACGGTAAAGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACCGTAAAAGCCATACACCGGCTGGCg  >  1:221559/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   GCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGATGCAGTACAGCTCACCGTGGCGGGTGGTTTCAATATCTATAGAAACCCACTTGAGCGGCGGACGATAGTCGGGATGCGGTTTCAGACGGGCATTAACGATAGTGCCATTGTGCATATCACCCTCGACCCACACCGGTGAGGTGATAAACCGCTCCATCAGATAGCGTTCTGGCGGACGCACATCGGCCTCGTAGACGGTAACGCCACCTTCACGCAGGCGCTTTTCGTAATTCATCAATTGGCGATGGGCGCGACAGTAAAGGCCATACACCGGCTGGCG  >  NC_000913/65071‑65558

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