Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,413,092	C→T	A48A (GCG→GCA)	intR ←	putative integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,092	0	C	T	98.0%	181.1 / ‑6.2	50	A48A (GCG→GCA)	intR	putative integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/21);  minor base G (1/0);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCGTTTGAAGTGGA  >  NC_000913/1412852‑1413337                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      ataACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAAt                                                                                                                                                                                                                                             >  2:268107/1‑251 (MQ=255)     cGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATGGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATGCACCGGAGTTTATTTCATGAACTATCAGTGCGCTTAAATTTcc                                                                                                                                                                                                                                         <  2:768016/251‑1 (MQ=255)          atGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCt                                                                                                                                                                                                                                    >  2:627488/1‑251 (MQ=255)            gtgATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCtt                                                                                                                                                                                                                                   <  1:201012/250‑1 (MQ=255)                      ggTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAAGATTAc                                                                                                                                                                                                                        <  2:272357/251‑1 (MQ=255)                      ggTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGATTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTAc                                                                                                                                                                                                                        >  1:742568/1‑251 (MQ=255)                        tACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTg                                                                                                                                                                                                                      >  1:621315/1‑251 (MQ=255)                         aCTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAAGTTCCGGCCTTTTTAATATTACTGt                                                                                                                                                                                                                     <  1:201640/251‑1 (MQ=255)                                      attattCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCa                                                                                                                                                                                                        <  1:791822/251‑1 (MQ=255)                                         attCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGcc                                                                                                                                                                                                     <  1:228496/251‑1 (MQ=255)                                              gAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCtttt                                                                                                                                                                                                 <  1:627488/250‑1 (MQ=255)                                                   tttttAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTTTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCATTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATATCCCCGCCTTTTAATGt                                                                                                                                                                                            <  2:52354/250‑1 (MQ=255)                                                       tAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTccc                                                                                                                                                                                       >  1:817879/1‑251 (MQ=255)                                                                tttCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGtt                                                                                                                                                                              >  2:601721/1‑251 (MQ=255)                                                                tttCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGtt                                                                                                                                                                              >  1:227401/1‑251 (MQ=255)                                                                    gTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTcc                                                                                                                                                                          >  1:162737/1‑251 (MQ=255)                                                                         cttcATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGa                                                                                                                                                                      >  1:38720/1‑250 (MQ=255)                                                                         cttcATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGa                                                                                                                                                                      >  1:35730/1‑250 (MQ=255)                                                                               tgtgTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGc                                                                                                                                                                                                      <  2:482844/212‑1 (MQ=255)                                                                               tgtgTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGc                                                                                                                                                                                                      >  1:482844/1‑212 (MQ=255)                                                                                                  gTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCt                                                                                                                                            >  2:116950/1‑251 (MQ=255)                                                                                                     tctTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTAcc                                                                                                                                         >  1:87904/1‑251 (MQ=255)                                                                                                                  tCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGac                                                                                                                             >  1:49654/1‑250 (MQ=255)                                                                                                                   ccAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGaca                                                                                                                            <  2:345434/250‑1 (MQ=255)                                                                                                                    cAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAAc                                                                                                                                                             >  1:19347/1‑216 (MQ=255)                                                                                                                    cAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAAc                                                                                                                                                             <  2:19347/216‑1 (MQ=255)                                                                                                                             caAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTa                                                                                                                  >  1:82518/1‑250 (MQ=255)                                                                                                                                 cctctCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTACTCGACTGCGGCAACCGTCTTAAATACTCAAATTCACCGGAGTTTCTTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCCGACGCGAATGTCTCTACCTCTAATCTCGACACCTGTTggcacttt                                                                                                              >  2:708376/2‑244 (MQ=255)                                                                                                                                       cGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAAcc                                                                                                                                                            >  2:206640/1‑198 (MQ=255)                                                                                                                                       cGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAAcc                                                                                                                                                            <  1:206640/198‑1 (MQ=255)                                                                                                                                         aaCGTTTTTATGACTCTCGTTGTCACCATTGTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTACGTTAGACATATc                                                                                                     <  2:75456/251‑1 (MQ=255)                                                                                                                                           cGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCAt                                                                                                   >  2:405786/1‑251 (MQ=255)                                                                                                                                             tttttATGACTCTCGTTGTCACCATTTATGCCCAAGTGCTGGACTGGGGAAAAAGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAAATATCAGTGCTCTTACATTTCCGGCCTTTTTAATATTACAGTTTGTAATCTCCCAGCCTTTTCCTGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTATAATCTCGACACCTGTTGGTACTTTAGACATATCATg                                                                                                  >  2:438546/1‑250 (MQ=255)                                                                                                                                                       tctcGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGAt                                                                                        >  1:646937/1‑250 (MQ=255)                                                                                                                                                               tCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGAtt                                                                               >  1:48929/1‑251 (MQ=255)                                                                                                                                                                             ccAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGt                                                                 <  2:35730/251‑1 (MQ=255)                                                                                                                                                                                       aCTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGAt                                                       <  1:405786/251‑1 (MQ=255)                                                                                                                                                                                          ggggAAAACGTCTTAAATACTCAAATTCACCGGATTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTTATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAAt                                                    <  2:561386/251‑1 (MQ=255)                                                                                                                                                                                            ggAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATcc                                                  >  2:212314/1‑251 (MQ=255)                                                                                                                                                                                                         aaTACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTc                                     >  1:382650/1‑251 (MQ=255)                                                                                                                                                                                                                 aaTTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTAc                             >  1:145643/1‑251 (MQ=255)                                                                                                                                                                                                                     cACCGGAGTTTAATTCATGAACAATCAGTGCTCTGAAATTTCCTGCCTTTTTAATATTACTGTTTGTAATCACCCAGCCTTTTAACGTTTCCCGACATCGTTTTCCCCGAAACATGAACCATATTCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCCGGAGACTTTGATAAACTGAGTTATCTTTGGATAGTTGTACCAGATATTCACTCGTTTGCTGTCTGGCATACCTa                          <  2:829765/250‑1 (MQ=255)                                                                                                                                                                                                                          gAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGgag                    <  2:382650/251‑1 (MQ=255)                                                                                                                                                                                                                            gTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGAt                  <  2:48929/251‑1 (MQ=255)                                                                                                                                                                                                                              ttATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATAc                >  2:261256/1‑251 (MQ=255)                                                                                                                                                                                                                               tATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATAc                <  2:227401/250‑1 (MQ=255)                                                                                                                                                                                                                                 tttCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCCGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCg             >  1:862224/1‑251 (MQ=255)                                                                                                                                                                                                                                   tCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCGt            <  2:646937/250‑1 (MQ=255)                                                                                                                                                                                                                                     aTGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTAGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCGTTTg         >  2:769146/1‑251 (MQ=255)                                                                                                                                                                                                                                             tCAGGGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGAAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGTTAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTAACTAAAGGAGATACTCGTTTGAAGTGGa  >  1:273183/1‑250 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      ATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAATCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCGTTTGAAGTGGA  >  NC_000913/1412852‑1413337

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