Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,032,219	G→A	Q353Q (CAG→CAA)	pepQ →	Xaa‑Pro dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,032,219	0	G	A	98.2%	199.7 / NA	57	Q353Q (CAG→CAA)	pepQ	Xaa‑Pro dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/31);  minor base C (1/0);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAGGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAA  >  NC_000913/4031976‑4032462                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    gCTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGAGCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATa                                                                                                                                                                                                                                              <  1:85602/251‑1 (MQ=255)   tGGTGAATGACGTTAATGATGAACAACTGGCGCTCATCGCCACCATGAAAGCAGGCGTCAACTATGTGTATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACAACAAATCATTACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGTGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGGCTGCAGGTGCATGACGTCGCTGTTTTTATGCAAGATGATAg                                                                                                                                                                                                                                             <  2:253026/250‑1 (MQ=255)        aaGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTAc                                                                                                                                                                                                                                       <  1:181041/251‑1 (MQ=255)             cGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGNAcc                                                                                                                                                                                                                                  >  2:358148/1‑251 (MQ=255)               ttAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTc                                                                                                                                                                                                                                >  2:329694/1‑251 (MQ=255)                tAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTc                                                                                                                                                                                                                                <  1:122299/250‑1 (MQ=255)                          caacTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGc                                                                                                                                                                                                                     >  2:148483/1‑251 (MQ=255)                               ggCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTAGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTTAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCGGGTTTTATGCAAGATGATAGCGGTACGCACCGCGCGTCACCGGCAAAat                                                                                                                                                                                                                <  2:169693/251‑1 (MQ=255)                                  gcTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATc                                                                                                                                                                                                             >  1:214452/1‑251 (MQ=255)                                          ccaccaTGAAAGCAGACGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGGCGCTGGTTTTATGCAAGATGGTAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTg                                                                                                                                                                                                     >  1:74820/1‑251 (MQ=255)                                                tGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGc                                                                                                                                                                                               >  1:144175/1‑251 (MQ=255)                                                       aGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGt                                                                                                                                                                                         >  2:228538/1‑250 (MQ=255)                                                                cTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGc                                                                                                                                                                               >  1:33414/1‑251 (MQ=255)                                                                  aTGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGcc                                                                                                                                                                              <  2:332409/250‑1 (MQ=255)                                                                  aTGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAAGCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCTCCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGc                                                                                                                                                                               <  1:220688/249‑1 (MQ=255)                                                                    gtgGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCGGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTAGCTGCGCTGCACCCGTATTCTCCAGCCgg                                                                                                                                                                            <  1:375355/250‑1 (MQ=255)                                                                      ggATTACCACATCCAGTTCCATCAGCGAATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGc                                                                                                                                                                          <  1:38944/250‑1 (MQ=255)                                                                           aCCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGgt                                                                                                                                                                     <  1:388684/250‑1 (MQ=255)                                                                               caTCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGCTATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTCTTATGCAAGATGATAGCGGTACGCAACTCGCGGCAACGGCAAAATATCCGTACCTGGGCTGCACCCGTACTCTCCAGCCGGGCATGGTGTTa                                                                                                                                                                 >  2:84813/1‑250 (MQ=255)                                                                                aTCCAGTGCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCGCCGATAGGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAAc                                                                                                                                                               <  1:233581/251‑1 (MQ=255)                                                                                   cAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCAt                                                                                                                                                            <  1:123749/251‑1 (MQ=255)                                                                                        ccATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAAc                                                                                                                                                       <  1:235836/251‑1 (MQ=255)                                                                                                   tCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGCGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTAc                                                                                                                                            >  1:352881/1‑251 (MQ=255)                                                                                                    cGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCg                                                                                                                                                          >  2:311324/1‑236 (MQ=255)                                                                                                    cGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCg                                                                                                                                                          <  1:311324/236‑1 (MQ=255)                                                                                                    cGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACt                                                                                                                                           >  2:117864/1‑251 (MQ=255)                                                                                                           ctgctgCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCg                                                                                                                                     >  1:148167/1‑250 (MQ=255)                                                                                                             gctgcGGAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTAGCTACTTCATCGAAt                                                                                                                                  <  1:118228/251‑1 (MQ=255)                                                                                                                   tAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTAc                                                                                                                            >  1:344132/1‑251 (MQ=255)                                                                                                                      aCATCAAATCATCACCGATATGAGTGAAGAAGTGATGGTCGAAGACGATCTTACCGAGGCGTTTATGCCGCATGGTATCGGCCATCCGCTGGACCTGCGGGTGCATGACGTCGCTGGTTTTATGCACGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGTATCGCTACTgg                                                                                                                         >  1:133355/1‑251 (MQ=255)                                                                                                                              tcatcaCCGATGTGAGTGAAGAAGCGATGGTCGCCAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTTCGCCTGCAGGTTCATGACGTCGATGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCCTGGCATTAACCAACGAACCGGGTATCTACATCATCGAATCACTACTGGCACCGTgg                                                                                                                 >  2:16713/1‑251 (MQ=255)                                                                                                                               catcaCCGATATGAGTTAAGAAGCGATGGTCGAAACCGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCTGTACGCACCTCGCGGCACCGGCAAAATATCCTTACCTGCGCTGCACCCGTATTCTCCAGCCGGGTATGGTGTTAACCAGCGAACCGGGTATCGACTTCATCGAATCGCTACTGGCACCGTGGc                                                                                                                <  2:233215/251‑1 (MQ=255)                                                                                                                                atcaCCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTNGCg                                                                                                               >  2:358916/1‑251 (MQ=255)                                                                                                                                         gtGAGTGAAGAATCGATTGTGGAAAGCGATCTTACCGGGCCGGTTATGCCGCATGGTATCGGCCATCCGCTTGGCCTGCAGGTGCATGACGTCGCTGGTTGTATGCAAGATGATAGCGGTGCGCACCTCGCGGCACCGGCAAAATATCCGTGCCTGCGCTGCACCCGGATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAATGGCa                                                                                                      <  2:7175/250‑1 (MQ=255)                                                                                                                                                 agaagCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCagca                                                                                              <  2:192936/251‑1 (MQ=255)                                                                                                                                                    agCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCagcaagc                                                                                           <  1:262832/251‑1 (MQ=255)                                                                                                                                                    agCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGAGCATGGTGTTAACCATCGAACCGGGCATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCagcaagc                                                                                           >  1:394362/1‑251 (MQ=255)                                                                                                                                                     gCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCagcaagca                                                                                          <  1:82115/251‑1 (MQ=255)                                                                                                                                                          ggTCGAAAACTATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAAGCGCTGCTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCa                                                                                     <  1:223309/251‑1 (MQ=255)                                                                                                                                                                aaaCGATCTTACCGGGCCGTTTATGCCGGATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGc                                                                               <  2:302088/251‑1 (MQ=255)                                                                                                                                                                        ttACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAACATATCCGTAcc                                                                                                                                                                                                       >  1:6/1‑123 (MQ=255)                                                                                                                                                                        ttACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAAtt                                                                       >  2:216391/1‑251 (MQ=255)                                                                                                                                                                        ttACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTAAGCACCTCGCGGCACCGGCAAAATATCCGTAcc                                                                                                                                                                                                       <  2:6/123‑1 (MQ=255)                                                                                                                                                                                 cGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTg                                                              >  1:299208/1‑251 (MQ=255)                                                                                                                                                                                     tATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAAc                                                          >  1:340312/1‑251 (MQ=255)                                                                                                                                                                                     tATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGAACTTCAACTGGCAGAAAATTGAAGCACTGaaa                                                           >  2:233549/1‑250 (MQ=255)                                                                                                                                                                                        gCCGCATGGTATCGGCAATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGAACCGGCAAAAGATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACGGGCAGAAAATTGAAGAACTGAAAcagt                                                       <  2:74820/251‑4 (MQ=255)                                                                                                                                                                                               ggTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTcggcg                                                 <  1:331439/250‑1 (MQ=255)                                                                                                                                                                                                           ccGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCg                                     <  1:358148/250‑1 (MQ=255)                                                                                                                                                                                                            cGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGa                                    >  2:297680/1‑250 (MQ=255)                                                                                                                                                                                                                   ccTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACg                            <  2:214452/251‑1 (MQ=255)                                                                                                                                                                                                                   ccTGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACg                            <  2:11590/251‑1 (MQ=255)                                                                                                                                                                                                                     tGCAGGTGCATGACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACgtg                          >  1:132117/1‑251 (MQ=255)                                                                                                                                                                                                                                gACGTCGCTGGTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACg                <  1:216391/250‑1 (MQ=255)                                                                                                                                                                                                                                         ggTTTTATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACg       <  2:269527/250‑1 (MQ=255)                                                                                                                                                                                                                                          gTTTTATGCAAGATGATAGCGGTACGCACATCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACTGTGGCGTGAAGGGCAGTTCAGCGAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTg     <  2:352881/251‑1 (MQ=255)                                                                                                                                                                                                                                              tATGCAAGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGaa  >  2:215359/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    GCTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCATGGTATCGGCCATCCGCTGGGCCTGCAGGTGCATGACGTCGCTGGTTTTATGCAGGATGATAGCGGTACGCACCTCGCGGCACCGGCAAAATATCCGTACCTGCGCTGCACCCGTATTCTCCAGCCGGGCATGGTGTTAACCATCGAACCGGGTATCTACTTCATCGAATCGCTACTGGCACCGTGGCGTGAAGGGCAGTTCAGCAAGCACTTCAACTGGCAGAAAATTGAAGCACTGAAACCGTTCGGCGGCATTCGTATCGAAGACAACGTGGTGATCCACGAAAACAACGTGGAA  >  NC_000913/4031976‑4032462

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