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.1%	190.2 / NA	53	A48A (GCG→GCA)	intR	putative integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/27);  minor base A (1/0);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCGTTTGAAGT  >  NC_000913/1412848‑1413334                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    cGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGctct                                                                                                                                                                                                                                                <  2:468458/250‑1 (MQ=255)       aaCGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAAttt                                                                                                                                                                                                                                         >  1:144526/1‑251 (MQ=255)       aaCGAATATGTGATAACGGGGTACTTACACAAATTATTATTATGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTATATTTAATTTTTGTCCCGATATCACAAAGATCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAACACGTCTTACATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGAGCTATTAAAttt                                                                                                                                                                                                                                         >  2:695251/1‑251 (MQ=255)          gAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTANTTCATGAACTATCAGTGCTCTTAAATTTCCg                                                                                                                                                                                                                                      >  2:339319/1‑251 (MQ=255)           aaTATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCgg                                                                                                                                                                                                                                     <  2:150422/251‑1 (MQ=255)                gtgATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCttt                                                                                                                                                                                                                                <  2:144526/251‑1 (MQ=255)                                      aaTTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAAtct                                                                                                                                                                                                          <  1:742232/251‑1 (MQ=255)                                             attCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGc                                                                                                                                                                                                    <  1:64028/250‑1 (MQ=255)                                                    gtgtTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAAt                                                                                                                                                                                            <  1:755054/251‑1 (MQ=255)                                                           tAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTccc                                                                                                                                                                                     >  2:214455/1‑251 (MQ=255)                                                                      tagTTTACTTCATTGTGTCTGTTGATAACTCTGTATCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGATTCTCGTTGTCACCATTTTTGCCACCGTGCTGGACTGGCGAAAACGTCTTCAATACTCAAATTCACCGGCTTTTCTTTCCTGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTACTATCCCAGCCTTTTAATGTTTCCAGACATCGattt                                                                                                                                                                          >  2:279389/3‑251 (MQ=255)                                                                           ttcttcATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTc                                                                                                                                                                      <  1:241141/250‑1 (MQ=255)                                                                               tcATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGaa                                                                                                                                                                   <  1:636621/249‑1 (MQ=255)                                                                                cATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATCTCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAAc                                                                                                                                                                 >  2:331187/1‑250 (MQ=255)                                                                                cATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTTTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTTTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTTTTAAATTTTCGCCTTTTTTATTATTACTGTTTGTATTTTCCCAGCTTTTTATTTTTTCCGGACATCGTTTTCCTCGAAAct                                                                                                                                                                >  2:487660/1‑250 (MQ=255)                                                                                       tttgttGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAAcc                                                                                                                                                          >  1:322290/1‑250 (MQ=255)                                                                                          gttggtAACTCTGACTCTTTAATTTATCTCCAGATCTCAAAAAGCTATCCGAACGTTTTTATGACTATCCTTGACACCATTATTGCCCCAGTGCTGGACTGGGGAAAACGTCTTTAATACTCCAATTCACCGGAGTTTCTTCAATGCACTATCAGTGTTCTTAAATTTCCGGCCTTTTTAATATTAATGTTTGTAATCTCCCAGCCTTTTAATGTTTACCGACATCGTTTTCCTCGAACCATGaccc                                                                                                                                                          >  2:694517/1‑244 (MQ=255)                                                                                                      gTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCt                                                                                                                                          >  1:235285/1‑251 (MQ=255)                                                                                                               aTTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑t                                                                                                                                <  1:339319/251‑1 (MQ=255)                                                                                                               aTTTTTGTCCAGAAATCACAAAGCTATCCGTACGTGTTTATGACTCTCGTTGTCCTCATCTTAGCCCCCGTGCTGGACTGGCGTACACGTATTAAATACTCATATTCACCGGAGTTTATTTCATGAACTATCAGTGCTATTACATTTCCGGCGTTTTTAATATCACTGTTTGTAATATCCAGGCCTTTTAACGTTTCCCGCCATAGTTTTCCTCGAAACATGACCCAAATGAGAATTTATATACCTATAA‑t                                                                                                                                >  2:686638/1‑251 (MQ=255)                                                                                                                 ttttGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑tct                                                                                                                              >  2:483575/1‑251 (MQ=255)                                                                                                                     gTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGac                                                                                                                          <  1:331187/251‑1 (MQ=255)                                                                                                                             atCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGtt                                                                                                                   <  1:6731/250‑1 (MQ=255)                                                                                                                              tCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTgg                                                                                                                 >  1:406777/1‑251 (MQ=255)                                                                                                                                acaAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGt                                                                                                                >  1:287243/1‑250 (MQ=255)                                                                                                                                       tctcCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACAcc                                                                                                                       >  1:442710/1‑236 (MQ=255)                                                                                                                                       tctcCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACAcc                                                                                                                       <  2:442710/236‑1 (MQ=255)                                                                                                                                          ccggaCGTTTGTAGGAATCTCGTTGTCACCATTTTTGCCCCAGTGCTGGAATGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGACCTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACGTTTTGTAAACTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCACTAC‑TCTCGACACCTGTTGGTAATTTATACag                                                                                                     <  2:76280/247‑2 (MQ=255)                                                                                                                                          ccGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTACATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCTGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCC‑ACATTGTTTTCATCGAAACCTGAACCAGATGCTAATGTCACTACATCTAAATC‑CGACACCTGTTTGTAATTTAGCCat                                                                                                     >  1:444727/1‑250 (MQ=255)                                                                                                                                          ccGAACGTTTTTATGACTCCCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACat                                                                                                     >  1:348150/1‑251 (MQ=255)                                                                                                                                            gAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGAGCTGGACAGGGGAAAACGTATTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCACAGACTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTATAA‑TCTCGACACCTGATGGTAATTTAGACatat                                                                                                   >  1:460866/1‑251 (MQ=255)                                                                                                                                                   tatATGACTCTCGTTGTCACGATTTTTGCCCCGGTGCTGGACTGGGGAAACCGTCTTAAATACTCAAATTCACCGGAGTTGATTTCATGATCTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATAGAGTTCCTCGAAACATTAACCAGATGCGAATGTATCCACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGt                                                                                            >  2:568833/3‑251 (MQ=255)                                                                                                                                                    ttATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGt                                                                                            <  1:243541/250‑1 (MQ=255)                                                                                                                                                      aTGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCtt                                                                                         >  1:431285/1‑251 (MQ=255)                                                                                                                                                        gACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTg                                                                                       <  1:487660/251‑1 (MQ=255)                                                                                                                                                        gACTCTAGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCCAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTACTATTACTGTTTGTACTCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCATCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGAAATATCATGAGTCTTTg                                                                                       >  1:736752/1‑251 (MQ=255)                                                                                                                                                               gTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACt                                                                                >  2:206286/1‑251 (MQ=255)                                                                                                                                                                     aCCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGAttt                                                                           >  1:457850/1‑250 (MQ=255)                                                                                                                                                                      ccATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGACCTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAACCTGATTTAt                                                                         >  2:116746/1‑251 (MQ=255)                                                                                                                                                                      ccATTTTGGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTTCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTGGTAATCTCCCAGCCTTTTAATTTTTCCC‑ACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCACGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGAATTATc                                                                        <  2:444727/251‑1 (MQ=255)                                                                                                                                                                       cATTTTTGCCCCAGTGCTGGACTGGGGAAAANGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATc                                                                        <  2:348150/251‑1 (MQ=255)                                                                                                                                                                            ttGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTg                                                                    <  1:279874/250‑1 (MQ=255)                                                                                                                                                                               ccccAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATa                                                                <  1:182305/251‑1 (MQ=255)                                                                                                                                                                                               gggAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTc                                                                                           >  1:183810/1‑208 (MQ=255)                                                                                                                                                                                               gggAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTc                                                                                           <  2:183810/208‑1 (MQ=255)                                                                                                                                                                                                ggAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATcc                                               <  1:461168/251‑1 (MQ=255)                                                                                                                                                                                                       gtcgtAAATACTCAAATTCACCGGAGTTTTTTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGttt                                        <  1:568833/247‑1 (MQ=255)                                                                                                                                                                                                                aCTCAAATTCACCGGAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTTTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCGTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTg                                <  2:692663/250‑1 (MQ=255)                                                                                                                                                                                                                              gAGTTTATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGgag                 <  2:287243/251‑1 (MQ=255)                                                                                                                                                                                                                                  ttATTTCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATa              <  1:116746/250‑1 (MQ=255)                                                                                                                                                                                                                                      ttCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCttt                                                                     >  1:361787/1‑191 (MQ=255)                                                                                                                                                                                                                                      ttCATGAACTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCttt                                                                     <  2:361787/191‑1 (MQ=255)                                                                                                                                                                                                                                             aCTATCAGTGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCGTTTGAAGt  <  2:406777/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    CGCTATAACGAATATGTGATATCGGGGTACTTTCACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAACATGAACCAGATGCGAATGTATCTACCTCTAA‑TCTCGACACCTGTTGGTAATTTAGACATATCATGAGTCTTTGATAAACTGATTTATCTTTGGATAGTTGTACCAGATAATCCCTCGTTTGCTGTCTGGCTTACCTAAAGGAGATACTCGTTTGAAGT  >  NC_000913/1412848‑1413334

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