Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,686	C→T	R173R (CGC→CGT)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,686	0	C	T	97.1%	217.6 / NA	70	R173R (CGC→CGT)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/36);  minor base G (0/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

ATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGAT  >  NC_000913/14444‑14931                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      atTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGagt                                                                                                                                                                                                                                               <  2:98101/251‑3 (MQ=255)    tttGGTGACGGTTGCGGCGATATGTTTGGCGGCGGACGGGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGGGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACGCAGTCGCAGACCGGGCCGACCTGGCATGGTTCTGGTGAGGGGCAGAGGCGTCAGGGAttct                                                                                                                                                                                                                                             <  2:224323/250‑1 (MQ=255)      tGGTGACGGTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGAttcttt                                                                                                                                                                                                                                           <  1:348856/250‑2 (MQ=255)          gACGTTTGCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGAttctttgccg                                                                                                                                                                                                                                       <  1:400238/250‑6 (MQ=255)             gttggcGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGAttctttgccgtgca                                                                                                                                                                                                                                   <  2:488105/246‑10 (MQ=255)                   ggCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGAc                                                                                                                                                                                                                             >  1:466678/1‑251 (MQ=255)                               ggcggcGGACTTGGTCGTCAACGTGCTGCGCGCGGTGCTGAATTACTCTAGATCAGGGAGTGCACCCTCGAAGAAGCTGGACGTGGCGTGACCAAAGAGATCCGCATTCCGTCTCTGGAAGAGTGTTACGTTGGTCACGGTAGCGGTGCAAAACAAGGTACACAGCCGCAGACCTGTCCGCCCTGTCATGGTTCTGGTCAGGTGCAGATGCTTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTg                                                                                                                                                                                                                 <  2:302240/251‑1 (MQ=255)                                        cGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAACCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCTGTGCAGCAGACCTGTCCACACTGTCAGGGCcg                                                                                                                                                                                                        >  2:139643/1‑251 (MQ=255)                                           ggTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCgg                                                                                                                                                                                                     >  1:426028/1‑251 (MQ=255)                                               gtcAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACg                                                                                                                                                                                                 <  2:346059/251‑1 (MQ=255)                                                 gaaCGTGCGGCGCGCGGTGATGATTTACGCTATAACATGGAGCACACCCTCGACGCAGCTGCGCGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGATGAGTGTGACGTTTGCCACGGTAGCGGCGGAAAACCAGGCAAACAGCCGCAGACCTGTCCGAACTGACATGCGTCTGGTCAGGTGCAGGTGCGTCAGGGCTTCTTGGCCGTGCAGCGAACCTGTCCACACTATCAGGCCCGGGGTACg                                                                                                                                                                                                 >  2:24310/2‑249 (MQ=255)                                                        cggcgCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCaa                                                                                                                                                                                         >  1:365702/1‑250 (MQ=255)                                                                ggTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTg                                                                                                                                                                                >  1:284245/1‑251 (MQ=255)                                                                ggTGCTGATTTACGCTAGAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGAGCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGGACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTg                                                                                                                                                                                <  2:216385/251‑1 (MQ=255)                                                                    cagattgacggtagaaCACGGAGCTCACCCCCGAAGAAGCTGTACGTGGCTGGACCAAAGAGAGCCGCATTCCGACTCTGGAAGGGCGGGACGTTTGCCACGGTTGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGCCCGGCCTGTCATGGTTCTGGTCAGGTGCAGATGCGGCAGGGATTCTTTGCCGTGCAGGAGACCTGGCCACACTGGCAGGGCCGCGGTACGCTGATCAAATATCCGTGcatc                                                                                                                                                                            <  2:164501/237‑3 (MQ=255)                                                                      gATTTACGCTATAACATGGAGCTCACTCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATGCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGcaacaa                                                                                                                                                                          <  2:47375/251‑1 (MQ=255)                                                                      gATTTACGCTATAACATGGAGCTCAACCTCGAAGAAGCTGTAAGTGGCGTGACCAAAGAGATCCGCCTTCCGACTCTGGAAGAGTGTGACGGTTGCCACGGTAGCGGTGCAAAACTAGGTACACAGCCGCAGGCCTGTCCGACCTGTCATGGGTCCGGTCATGTGCAGATGCGGCAGGGATTCTTGGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGTTACGCTGATCAAATATCCGTGcaacaa                                                                                                                                                                          <  2:73375/251‑1 (MQ=255)                                                                         ttACGCTATAACATGGAGCTCACACTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAAt                                                                                                                                                                        <  2:352310/250‑1 (MQ=255)                                                                           aCGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGATCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtc                                                                                                                                                                     >  1:440539/1‑251 (MQ=255)                                                                            cgctagaacagtgAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGATATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCATATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCGGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtca                                                                                                                                                                    <  2:460792/239‑1 (MQ=255)                                                                             gCTATAACAGGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGGGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtcat                                                                                                                                                                   <  2:365702/251‑1 (MQ=255)                                                                                   aCATGGAGCCCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtcatggtcat                                                                                                                                                             <  2:197633/251‑1 (MQ=255)                                                                                      tGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTcatggtcatggt                                                                                                                                                          >  1:476917/1‑251 (MQ=255)                                                                                        gAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCAGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCTGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCg                                                                                                                                                        <  1:128996/251‑1 (MQ=255)                                                                                          gCTCACCCTCGAAGAAGCTGTACGTGGCGTGCCCAAAGCGATCCGCATTCCGACTCTGGCAGCGTGTGACGTTTGCCACGGTAGCGGTGCAACACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCCGCCCTGACCACACTGTCAGGGCCGCGGGCCGCTGATCACAGATCCGAGCAACCAATgccatggtcaaggcggtg                                                                                                                                                      >  2:261103/1‑243 (MQ=255)                                                                                           cTCACCCTCGAAGAGGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGCCCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGGACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCgtgt                                                                                                                                                     >  2:393782/1‑251 (MQ=255)                                                                                            gcACCCTCGAAGAAGCTGTACGTGGCGTGGCCAAGCAGATCGGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAACACCCGGTACACAGCCGCAGACCTGTCAGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTTCGATGATAACAGAGCCGTGCACCGAATgtcatgtgcatggcagtgtt                                                                                                                                                    >  2:41993/2‑237 (MQ=255)                                                                                                ccTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAgc                                                                                                                                                <  1:425599/251‑1 (MQ=255)                                                                                                   cGAAGAAGCTGTACGGGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAgcg                                                                                                                                               <  2:138341/249‑1 (MQ=255)                                                                                                     aagaagCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCgcagc                                                                                                                                           >  2:184060/1‑251 (MQ=255)                                                                                                      agaagCGGTACGTGGCGTGACCAAAGAGATCCGCATGCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCgcagca                                                                                                                                          <  2:466678/251‑1 (MQ=255)                                                                                                        aagCTGTACGTGGCGTGCCCAAAGAGATCCGCATTCCGACTCTGGCAGAGTGTGACGCTTGCCACGGTAGCGGTGAAAAACCAGGTACACAGCAGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGAGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGCACGCTGATCAAAGATCCGCGCAACAATTGTCATGGGCATGGTCGTGTTGTGCGCAGCaac                                                                                                                                        >  2:467426/1‑250 (MQ=255)                                                                                                        aagCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCaaa                                                                                                                                        <  1:498596/251‑1 (MQ=255)                                                                                                          gCTGTACGTGGCGTGACCAAAGAGATCCGCTTTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTGCACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATTGTCATGGTCGTGTTGAGCGCAGCAAAAc                                                                                                                                      <  1:97247/251‑1 (MQ=255)                                                                                                          gCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAAc                                                                                                                                      <  2:237346/251‑1 (MQ=255)                                                                                                             gTACGGGGCGTTACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCt                                                                                                                                   <  2:354522/251‑1 (MQ=255)                                                                                                                  tggggcgACCAAAGAGAGCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAACCCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTGCGCTGATCAAAGATCCGTGCAACAAATGTCATGGACATGGTCGTGTTGAGCGCAGCAAGAAGCTGTCCg                                                                                                                              >  2:333074/7‑251 (MQ=255)                                                                                                                     cGTGACCAAAGAGATCCGCATTCCGACTCTGGCAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAACCGCTGCCCGtt                                                                                                                            >  2:482783/1‑250 (MQ=255)                                                                                                                        gACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCGCGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTaaa                                                                                                                         <  1:460188/250‑1 (MQ=255)                                                                                                                                 gaTCCGCATTCCGACTCTGGAAGAGTGTGACGGTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAg                                                                                                               <  1:450241/251‑1 (MQ=255)                                                                                                                                            cGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACt                                                                                                    <  1:184060/251‑1 (MQ=255)                                                                                                                                            cGACTCTGGAAGAGTGTGACGTTTACCACGGTAGCGTTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCCTGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGAGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCTTGTTGGGCGCAGCAAAACGCTGTCCGTAAAAATACCGGCGGGGGTGGacacg                                                                                                    >  2:401246/1‑250 (MQ=255)                                                                                                                                             gACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTCCGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGGCACTg                                                                                                   >  2:15134/1‑251 (MQ=255)                                                                                                                                                tctGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtcatggtcatg                                                                                                                                                            >  1:105879/1‑191 (MQ=255)                                                                                                                                                tctGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtcatggtcatg                                                                                                                                                            <  2:105879/191‑1 (MQ=255)                                                                                                                                                      agagTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGc                                                                                           >  1:363770/1‑250 (MQ=255)                                                                                                                                                            tgACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCGGCAAAACGCTGTCCGTTAAAATCCCGGCCGGGGGGGACACTGGAGACCGCaccagtc                                                                                    >  2:104682/1‑245 (MQ=255)                                                                                                                                                                   tGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCggg                                                                             >  1:481109/1‑251 (MQ=255)                                                                                                                                                                   tGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCgg                                                                              <  1:502216/250‑1 (MQ=255)                                                                                                                                                                        cGGTAGCGGTGCAAAAACAGGTACACAGCAGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTAATGGTCGTGTTGAGCGCGACAACACGCTGTCCGTTACAATCCCGGCATGGGTGGGCACTGGGGGCCGCATCCGTCTTGCGGGCGagg                                                                        >  2:48435/1‑249 (MQ=255)                                                                                                                                                                           tAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTg                                                                     >  1:330622/1‑251 (MQ=255)                                                                                                                                                                            aGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATccc                                                                                                                    >  1:197913/1‑203 (MQ=255)                                                                                                                                                                            aGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGGCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATccc                                                                                                                    <  2:197913/203‑1 (MQ=255)                                                                                                                                                                                 tGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCgg                                                               >  1:465633/1‑251 (MQ=255)                                                                                                                                                                                      aaCCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAg                                                          <  2:426028/251‑1 (MQ=255)                                                                                                                                                                                          aGGTACACAGCCGCATACCTGTCCGCCCTGTCATGGTTCTGTTGAGGTGCAGATGTGTCAGGGATTCTTTGCCGTGTAGCTGACCGGTCCACACTGTCAGGGCCGCGGTACGCTGATCAACGATCCGTGCAACAAATGTCATGGTCATGGGTGTTGTGAGTGCAGCAAAACGCTGTCCGGTAAAATCCTGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCgcgaatg                                                      <  2:416153/251‑7 (MQ=255)                                                                                                                                                                                                  aGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGGGGACACTGGAGACCGAATCCGTCTTGCGGGCGGAGGGGAAGCGGGCGGGCATGGCGCCCCg                                              >  2:160874/1‑251 (MQ=255)                                                                                                                                                                                                            ctGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGACGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATc                                    >  1:433006/2‑251 (MQ=255)                                                                                                                                                                                                                ccGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTa                                <  1:458677/251‑1 (MQ=255)                                                                                                                                                                                                                  gACCTGTCATGGTTCTGGTCAGGTGCATATGCGTCAGGGATTCATGGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCTGTACGCGGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCTCACCGGCAGGCGATCTGTACg                              <  2:464669/251‑1 (MQ=255)                                                                                                                                                                                                                  gACCTGTCATGGTTCTGGACAGATGAAGATGCGTCAGGGCTTCTTTGCCGAGCAGCAGACCTGTCCAAACTGTCAGGGCCGCGGTACGCTGATCACAGATCCGTGCAACCACTGTCATGGTCATGGTCGTGTTTAGCGCAGCACAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGCGACCGCATCCGTCTTGCGGGCGACGGCGCCGCGGGCGGGCATGGCGACCCGGCATGCgctccgtccg                              >  2:59544/1‑242 (MQ=255)                                                                                                                                                                                                                   aCCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgt                             >  1:193724/1‑251 (MQ=255)                                                                                                                                                                                                                     cTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgttc                           >  1:119606/1‑251 (MQ=255)                                                                                                                                                                                                                       gTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgttc                           >  1:259855/1‑249 (MQ=255)                                                                                                                                                                                                                         cATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgttcaggt                       >  1:159008/1‑251 (MQ=255)                                                                                                                                                                                                                            ggTTCTGGTCAGGTGCAGATGGGTCAGGGATTCTTTGCCGTGCAGCAGCCCTGTCCACTCTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGGGCAACAAAGGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgttcaggttca                    <  2:300110/251‑1 (MQ=255)                                                                                                                                                                                                                                     cAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAAc            <  1:15134/250‑1 (MQ=255)                                                                                                                                                                                                                                      aGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAg          <  1:478831/251‑1 (MQ=255)                                                                                                                                                                                                                                      aGGTGCAGAGGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACTCTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAg          <  2:396825/251‑1 (MQ=255)                                                                                                                                                                                                                                              aTGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGAt  <  1:441951/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      ATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGAT  >  NC_000913/14444‑14931

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