Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,956	G→C	L263L (CTG→CTC)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,956	0	G	C	93.8%	205.4 / 6.2	65	L263L (CTG→CTC)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base C (26/35);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTG  >  NC_000913/14709‑15199                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    agACCTGTCCACACTGTCAGGGCCGCGTTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTGGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGCGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGTAACAACCTCTAt                                                                                                                                                                                                                                                  <  2:326285/251‑1 (MQ=255)    ccattCCACACTGTCAGGGCCGCCGTACGCTGTCCAAAGATGCGTGCAACAAAATTCATGGTCATGGTCTTGTTGATCGCAGCAAAACGCTGTCCGTTAAATGCCCGGCAGGGGTGGACATTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCTGGTGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTCCAGGTTAAACAGCCC‑CCGATTTTCGAGCGTGAAGTCAACAACCTCTATTGc                                                                                                                                                                                                                                               <  2:265105/247‑1 (MQ=255)             acTGTCAGGGCCGCGGTACGTGGATGAAAGACCCGTGCAACAAATGTCAGGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAAGCCCGGCAGGGGTGGGCACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCg                                                                                                                                                                                                                                      <  2:56594/251‑1 (MQ=255)                gTCAGGGCCGCGGAACGCTGATACAAGCTCCGTGCAACAAATGTCATGGTCATGGTCGCGTTGCGCGCAGCAAAACGCTGTCCGTTACACTCCCGGCAGGGGTGGACACTGGAGACAGCCTCCGTCTTGCGGGAGAAGGTGAAGCGGGCGCGCATGGACCACCGGCAGGCGATCTGACCGTTCCGGTCCAGGTTACACAG‑ACACCGATTTTCGAGCGTGAAGGCAACAACCTCTAGTGCGAGGTCCCGATc                                                                                                                                                                                                                                   >  2:58281/1‑251 (MQ=255)                                   gATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGGAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATATTCGAGCGTGAATGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGgcggcg                                                                                                                                                                                                                 <  1:204580/250‑1 (MQ=255)                                      cAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTg                                                                                                                                                                                                              <  1:232400/250‑1 (MQ=255)                                               gTGCAACAAATGTCATGGTCATGGGCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGaaa                                                                                                                                                                                                    <  1:103257/251‑1 (MQ=255)                                                 gCAACAAATTTCATGGTCATGGTCGGGTTGAGCGCAGCAAAACGCTGTCCGGTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCGTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGTGCTGGGTGGCGAAATc                                                                                                                                                                                                  <  2:174681/251‑1 (MQ=255)                                                 gCAACAAATGTCATGTTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACGGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCCC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATGGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATc                                                                                                                                                                                                  <  2:474215/251‑1 (MQ=255)                                                          ggcatggtcatggtCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGTCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATGTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGAGTAACTTCGCGATGGCGGCGCTGGGTGGCGAAATCGAAGTACCg                                                                                                                                                                                         <  1:277735/249‑1 (MQ=255)                                                                gtcatggtcGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCtt                                                                                                                                                                                   >  1:85164/1‑251 (MQ=255)                                                                  catggtcGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATggcg                                                                                                                                                                                                                    <  1:102215/216‑1 (MQ=255)                                                                     ggtcGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGAGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCGGGGGGGCGGCATCGAGGTACCGACCCTTGATgg                                                                                                                                                                              >  1:153326/1‑251 (MQ=255)                                                                         gtgtTGAGCTCAGCAAACCGCTGTCCGTTAAAATCACGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGCTCAGCCTCAGGTTAAACAGCAA‑CCGGTTTTCGAGCGAGACGGCAACACCATCTCTTGCGAAGTCCCGCTCACCTTCGCTATGGCGGCGCCGGGAGGCGAAAACGAAGTAcagccccacgatgcaccc                                                                                                                                                                          >  2:245482/1‑234 (MQ=255)                                                                             ggAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTc                                                                                                                                                                       <  1:94074/249‑1 (MQ=255)                                                                                       aaaaCGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGt                                                                                                                                                            >  1:86112/1‑251 (MQ=255)                                                                                         aaCGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCtt                                                                                                                                                                                   <  2:209649/226‑1 (MQ=255)                                                                                         aaCGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCtt                                                                                                                                                                                   >  1:209649/1‑226 (MQ=255)                                                                                         aaCGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGc                                                                                                                                                          >  1:86519/1‑251 (MQ=255)                                                                                           cGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTCTGGCGGCACTGGGGGGCGAAATCGAAGGACAGCCCCTTGATGGTCGCCTCAAACTGAAAGAGccg                                                                                                                                                        >  1:14086/1‑250 (MQ=255)                                                                                            gCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTg                                                                                                                                                       <  1:457535/251‑1 (MQ=255)                                                                                                   gTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAAATGCAAATGCCTGGCGAAAc                                                                                                                                                >  1:209186/1‑251 (MQ=255)                                                                                                      aaaaTCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAAcccc                                                                                                                                             >  1:405601/1‑250 (MQ=255)                                                                                                             cGGCAGGGGTGGACACTGGAGACCGCATCCGTCGTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAAGCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGt                                                                                                                                      <  2:24568/251‑1 (MQ=255)                                                                                                              acgacgcgcgtccgatctagaCCGCATCCGTCTTGCGGGCGAAGGGGAAGCTGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGAGTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGGCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTa                                                                                                                                     <  2:230641/233‑1 (MQ=255)                                                                                                                cAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAg                                                                                                                                   <  1:2218/251‑1 (MQ=255)                                                                                                                cAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTa                                                                                                                                     >  1:104783/1‑249 (MQ=255)                                                                                                                cAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAg                                                                                                                                   <  1:232665/251‑1 (MQ=255)                                                                                                                  ggggTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGGGGGGAAATCGAAGTAACGACCCTTGTTGGTCGCGTCAAACTGAAAGGGCCTGGCGTAACCCAAACCGGTAAGCt                                                                                                                                 >  1:387589/1‑251 (MQ=255)                                                                                                                    ggTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCttt                                                                                                                               <  1:12737/251‑3 (MQ=255)                                                                                                                      gtgACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGACGCACCGGCAGGCGATCTGTACGGTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCGCTATTGCTAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTc                                                                                                                             <  1:34804/249‑1 (MQ=255)                                                                                                                        gACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCg                                                                                                                           <  2:29874/251‑1 (MQ=255)                                                                                                                                agaCCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGACACCCAGACCGGTa                                                                                                                                     >  1:230641/1‑233 (MQ=255)                                                                                                                                   ccGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAAc                                                                                                                                                >  1:143261/1‑219 (MQ=255)                                                                                                                                   ccGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAAc                                                                                                                                                <  2:143261/219‑1 (MQ=255)                                                                                                                                      cATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCAACGGCAGGCAATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAg                                                                                                             <  1:350378/251‑1 (MQ=255)                                                                                                                                        tCCGTCTTGCAGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGCTTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGCAATCTACGTACCGACCCTTGATGGTCGCGTCAAACTTCAAGGGCCTGGCGAAACCCCGACCGGTAAGCTTTTCCGTCTGCGCGGTCAAGGc                                                                                                           >  2:192750/1‑251 (MQ=255)                                                                                                                                         ccggcTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCGGTACGTTCAGGTTAAGGTGAAACAGCAC‑CTGATGTTCGAGCGTGAAGGCAACAACATCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGATGGGTGGCGAAATCGAAGTACCGACCCTGGATGGTCGCGTCAAACTGAAAGTGCATGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCg                                                                                                          <  2:38957/247‑1 (MQ=255)                                                                                                                                               tGCGGGGGAAGGTGAAGCGGGCGGGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCa                                                                                                       <  2:103674/248‑1 (MQ=255)                                                                                                                                                   ggCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCa                                                                                                                                             >  1:145253/1‑206 (MQ=255)                                                                                                                                                   ggCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCa                                                                                                                                             <  2:145253/206‑1 (MQ=255)                                                                                                                                                      gAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTGTGCGCGGTAAAGGCGTCAAGTCTGTCcg                                                                                             <  1:375088/251‑1 (MQ=255)                                                                                                                                                                cGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGgcgcc                                                                                   >  1:177521/1‑250 (MQ=255)                                                                                                                                                                      atcagggagcACCGGCATGCGATCTGGTCGTTTAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGATGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGt                                                                             <  2:453689/243‑1 (MQ=255)                                                                                                                                                                                   ggCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGcc                                                                <  1:134947/251‑1 (MQ=255)                                                                                                                                                                                      aGGCGATCTGTACGTTCAGGTTCCGGTTAAGCAGCTC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTAGTGCGAAGTCCCGATCAACGTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCt                                                                     <  2:397597/243‑1 (MQ=255)                                                                                                                                                                                      aGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCcgcg                                                             <  2:85164/251‑1 (MQ=255)                                                                                                                                                                                      aGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGGGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGAGCCTGGCGAAACCCAGACCGGTAAGCTTTTACGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGGGGCGCACAGGGTGATTTGCt                                                                     >  1:397597/1‑243 (MQ=255)                                                                                                                                                                                       ggCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGATTTGCTGTGCCGCGt                                                            >  1:419806/1‑251 (MQ=255)                                                                                                                                                                                            tCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGGGGCGCCCAGGGTGATTTGCTGTGCCGCGTTttcg                                                       >  2:177084/1‑251 (MQ=255)                                                                                                                                                                                               gTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTgtcgtc                                                     <  1:13654/250‑1 (MQ=255)                                                                                                                                                                                                      caggttcaggttAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAacac                                              <  1:149157/250‑1 (MQ=255)                                                                                                                                                                                                      caggttcaggttAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACAcc                                             >  1:397741/1‑251 (MQ=255)                                                                                                                                                                                                       atgtgcaggttAAACTGCTC‑CAGATTTTCGAGCTTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACGTTGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAAACCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACAcc                                             <  2:339571/245‑1 (MQ=255)                                                                                                                                                                                                                  aaaCAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGa                                  <  1:177084/250‑1 (MQ=255)                                                                                                                                                                                                                   aaCAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGGATGCGCGGTAAAAGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGATTTGCTGTGCCGAGGTGTCGTCGAAACCCCGGGAGGCCTGAAc                                >  2:413903/1‑251 (MQ=255)                                                                                                                                                                                                                      aGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGag                             <  1:15297/251‑2 (MQ=255)                                                                                                                                                                                                                        ccc‑ccGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGagaa                           <  2:405601/249‑4 (MQ=255)                                                                                                                                                                                                                          c‑ccGATTTTCGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACgggcagc                         >  2:30496/1‑245 (MQ=255)                                                                                                                                                                                                                                ttttAGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCCCTGGGGGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGAATACGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGCACCACCGGTAGGCCTGAACGAGAAGCAGaa                     >  1:36252/1‑250 (MQ=255)                                                                                                                                                                                                                                 tttCGAGCGTGAAGGCAACAACCTCTCTTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAAc                   >  1:184584/1‑251 (MQ=255)                                                                                                                                                                                                                                    cGAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGGAAAGGCGTCAAGTCTGTCCGCGGTGGCGCCCAGGGGGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAAccac                >  2:325378/1‑248 (MQ=255)                                                                                                                                                                                                                                         gTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGctgctg           >  1:327153/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aacaacCTCTATTGCGAAGTCCCGATCAACTGCGCTAGGGCGGCGCTGGGTGGCGAAAGCGAAGTACCGACCCTTGATTGTCGCGTCAAACTGAAAGTGCCTGGCTAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGGGTTGTCGTCGAAACACCGGTAGGCCTGAACGAGAAGCAGAAACAGCTgctgcaggagct   <  2:160403/251‑1 (MQ=255)                                                                                                                                                                                                                                                  acaacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGGGATTTGCTGTGCCGCGTTGTCGTCGAAACCCCGGTAGGCCCGAACGAGAAGCAGAAACAGCTgctgcaagagctg  >  2:128436/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCAC‑CCGATTTTCGAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGATTTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACGAAAGGCAGAAACAGCTGCTGCAAGAGCTG  >  NC_000913/14709‑15199

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: ‑