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	98.6%	211.5 / ‑4.9	71	R173R (CGC→CGT)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/42);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

ATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATT  >  NC_000913/14442‑14932                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        atatTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGa                                                                                                                                                                                                                                                   <  1:668959/251‑1 (MQ=255) atatTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGa                                                                                                                                                                                                                                                   <  1:173501/251‑1 (MQ=255)    tttttGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGAttc                                                                                                                                                                                                                                                >  1:662648/1‑251 (MQ=255)      tttGGTGACGGTTTCGGCGATAGTGTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTTGACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGAttct                                                                                                                                                                                                                                               <  1:398192/250‑1 (MQ=255)       ttGGTGACGTTTTCGGCGATATTTTTGGCTGCGGACGTGGGCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCTCCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCATACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGAttctt                                                                                                                                                                                                                                              <  2:663216/250‑1 (MQ=255)                    cGGCGATATTTGTGGGGGCGGACGTGGTCGTCACCGTGCGGCGCGCGGTGCTGATTTCCGCTAGAACATGTAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGGCCAAAGAGAGCCGCATTCCGACTCTGGAAGAGTGTG‑ACGGTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGGCATGGTTCGGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTTcagc‑ag                                                                                                                                                                                                                                <  2:445905/250‑1 (MQ=255)                       cGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCt                                                                                                                                                                                                                            <  1:353713/251‑1 (MQ=255)                                tGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCacac                                                                                                                                                                                                                    <  1:89264/250‑1 (MQ=255)                                      gggACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGTATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAgg                                                                                                                                                                                                             <  1:206518/250‑1 (MQ=255)                                         gcGTGGGCGGCAACGGGGGGCGCGTGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGTCCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACTCTGTCAGGGcc                                                                                                                                                                                                          <  2:477125/250‑1 (MQ=255)                                          cGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCcg                                                                                                                                                                                                         <  1:277735/251‑1 (MQ=255)                                           gTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCcgc                                                                                                                                                                                                        >  1:43044/1‑251 (MQ=255)                                                cgtcAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTa                                                                                                                                                                                                    <  1:379560/250‑1 (MQ=255)                                                    aaCGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTg                                                                                                                                                                                               >  1:137006/1‑251 (MQ=255)                                                            gcgcgcGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGa                                                                                                                                                                                       <  1:119651/251‑1 (MQ=255)                                                               cgcgGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGCAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACACTGATCAACGATcc                                                                                                                                                                                    >  2:229076/1‑251 (MQ=255)                                                                 cgGTGCTGATTTACGCTATAAAATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCTTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCGGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGt                                                                                                                                                                                  <  2:279058/251‑1 (MQ=255)                                                                     gCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGca                                                                                                                                                                               <  1:437208/250‑1 (MQ=255)                                                                          tttACACTATAACCTGGAGCTAACCCTCGAAATCGCTGTCCGTGGCGCGACCACAGAGATCCGCATTCCGACTCTGGGAGAGTGTG‑ACGTTTGCCACGGCAGCGGTGCACAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGACGTGCAGCCAG‑CCTGTCAACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAAt                                                                                                                                                                         >  2:668427/1‑251 (MQ=255)                                                                             aCGCTATAACATGGAGCTTACCCTCGAAGAAGCTGTACGTGGCGGTACCAAAGAGATCCGCATGCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgt                                                                                                                                                                       <  2:678971/250‑1 (MQ=255)                                                                              cGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtca                                                                                                                                                                     <  1:439230/251‑1 (MQ=255)                                                                               gCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCTGTGCAACAAATgtcat                                                                                                                                                                    <  1:514616/251‑1 (MQ=255)                                                                                cattaACATGGAGCTCATCCTCGCATAAGCGTTACGTGGCGTGAACAAATAGATCCGCATTCCGACCCTGTAAGTGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCTACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCAGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGGGCAACAAATgtcatg                                                                                                                                                                   <  2:150789/248‑1 (MQ=255)                                                                                       aTGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGtcatggtcatgg                                                                                                                                                            >  1:629722/1‑251 (MQ=255)                                                                                         ggAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGAGCCGCATTCCGACTCTGGAAGAGTGTG‑ACGGTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTGTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTcatggtcatggt                                                                                                                                                           <  1:93515/250‑1 (MQ=255)                                                                                              tCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGTTCACAGAt                                                                                                                                                                                      >  2:188067/1‑218 (MQ=255)                                                                                              tCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGAt                                                                                                                                                                                      <  1:188067/218‑1 (MQ=255)                                                                                              ccACCCTCGAATAAGCTGTACGTGGCGTGACCAAAGAGAGCCGCATTCCGCATCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGTTGCAAAACCAGGTACACAGCCGCAGACCGGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCGGGGCCGCGGTACGCTGATCAAAGATCCGAGCAACAAATGGCATGGTCATGGTCGTGtt                                                                                                                                                     <  2:22600/250‑1 (MQ=255)                                                                                                     cGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAgcg                                                                                                                                                <  1:594888/249‑1 (MQ=255)                                                                                                     cGAAGAAGCTGGACGGGGCGTGACCAAAGAGATCCGCATTCCGACTCTTGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGAACTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTGGAGCgca                                                                                                                                              <  2:286495/251‑1 (MQ=255)                                                                                                      gacgaagCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGAttctttgc                                                                                                                                                                                                                                           >  1:566580/4‑155 (MQ=255)                                                                                                      gaagaaGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGAttctttgc                                                                                                                                                                                                                                           <  2:566580/158‑4 (MQ=255)                                                                                                        agaagCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCgcagc                                                                                                                                            >  1:94133/1‑250 (MQ=255)                                                                                                        agaagCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGAACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGAgcagg                                                                                                                                            >  2:474798/1‑249 (MQ=255)                                                                                                                 aCGTGGCGTGACCAAAGAGATCCGTATTCCGACTCTTGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTTTCAGGGCCGCGGTACGATGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCATCAAAACGCTGt                                                                                                                                  <  2:393207/251‑1 (MQ=255)                                                                                                                     ggcggtaCCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTATCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCGGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGt                                                                                                                              <  2:648939/245‑1 (MQ=255)                                                                                                                          gACCAAAGAGTTCCGCATTCCGACTCTGGAAGGGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTTGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTaaa                                                                                                                          <  2:14059/250‑1 (MQ=255)                                                                                                                            ccAAAGAGATCCTCATTACGATTCTGGCAGAG‑GTGTACGTTTGCCTCGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGGCATGGGTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGGCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATc                                                                                                                       <  2:449685/251‑1 (MQ=255)                                                                                                                                   gaTCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGtt                                                                                                                                                     >  2:594568/1‑214 (MQ=255)                                                                                                                                   gaTCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGtt                                                                                                                                                     <  1:594568/214‑1 (MQ=255)                                                                                                                                       cGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAggggg                                                                                                            >  1:36208/1‑250 (MQ=255)                                                                                                                                           ttCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACCAATGTCATGGTCATGGTCGTGTTGAGCGAAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGac                                                                                                        >  1:434790/1‑251 (MQ=255)                                                                                                                                           ttCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCTTGTTGAGCGCAGCAAAACGCGGTCCGTTAAAATCCCGGCAGGGGTGGac                                                                                                        >  1:434745/1‑251 (MQ=255)                                                                                                                                             ccGACTCTGGAAGAGTGTG‑ACGTTTACCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTCAAATCCCGGCAGGGGTGGacac                                                                                                      >  1:107125/1‑251 (MQ=255)                                                                                                                                                aCTCTGGAAGAGGGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTAAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGACGCGGTACGCTGATCAAAGATCCGTGCAACAACTGTCATGGGCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGGGGGCACTgg                                                                                                   >  2:494426/1‑251 (MQ=255)                                                                                                                                                        agagTGTG‑ACGTTTGCCACGGAAGCGGAGCAAGACCAGGTACACAGCCGCAGACTTGGCCGACCTGTAATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGAACAGC‑AAACCTGTCCACACTGTCAGAGCCGCGGTACGCTGATCAAAGACCCGTGCAACAAATGTCATGCTCATGGTCGTGTTGCGCGCACAAAAACGCTATCAGGTAAAATCCCGGCGGGGGTGGACCCAGGCGACCGcc                                                                                           >  2:351149/1‑250 (MQ=255)                                                                                                                                                           gtgtG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGGCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGTACACTGGAGACTGCATcc                                                                                        <  2:94133/251‑1 (MQ=255)                                                                                                                                                             gtg‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGt                                                                                      >  1:581874/1‑251 (MQ=255)                                                                                                                                                               g‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCt                                                                                    >  1:367476/1‑251 (MQ=255)                                                                                                                                                                 aCGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACCAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACCATGGGGACCGCATCCGTCtt                                                                                   >  1:532170/1‑251 (MQ=255)                                                                                                                                                                  cGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGACGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTg                                                                                  >  1:549574/1‑251 (MQ=255)                                                                                                                                                                      tGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGAGGGAAACTGGAGACCGCATCCGTCTTTCggg                                                                              >  2:642351/1‑251 (MQ=255)                                                                                                                                                                         cACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCg                                                                            >  1:512015/1‑250 (MQ=255)                                                                                                                                                                            ggTAGCGGTGCAAAACAAGGCACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCAGCGGTACGCTGATCAACGATCCGTGCAACACATGTCATGGTCATGGTCGTGTTGAGCGCAGCTAAACGCTGTCCGTTAAAATCCCGGCAGGGGGGGACCCTGGAGAcgacatgcggattgaggtcgacgg                                                                        >  2:310865/1‑228 (MQ=255)                                                                                                                                                                                gCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAAGCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGaa                                                                    <  1:406188/251‑1 (MQ=255)                                                                                                                                                                                     gCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCggg                                                               <  1:132665/251‑1 (MQ=255)                                                                                                                                                                                     gCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCgg                                                                <  1:561286/250‑1 (MQ=255)                                                                                                                                                                                         aaCCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTACTGGGCAGGTGCAGATGCGTCAGGGATTCGTTGCCGTGCAGC‑AGACCGGTCCACACTGTCAGGGCCGCGGTACGCTGAACAAAGATCCCTGCAACAAATGTCATGGTCATGGGCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTCGACACTGGAGACCGCATCCGTCTTGCGGGc                                                                             <  1:257726/233‑1 (MQ=255)                                                                                                                                                                                                gacacAGCCGCAGACCTGGCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGgcg                                                    <  1:429653/250‑1 (MQ=255)                                                                                                                                                                                                 acacAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGAGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AACCCTGTCCACACTGTCAGAGCCGCGGTACGCTGATAAAAGATCCGTGCAACAAATGTCATGGACATGGTCATGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATACCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGGAGGAGAAGCGGGCGTGCATGgcgc                                                   >  1:630821/1‑251 (MQ=255)                                                                                                                                                                                                              cctGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGAt                                      <  1:370282/249‑1 (MQ=255)                                                                                                                                                                                                              cctGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGAt                                      >  1:44961/3‑251 (MQ=255)                                                                                                                                                                                                                tGTCCGACCTGTCCTGGCTCTGGACCGGTGCAGCTGCGTCAGGGCTTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGCTCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGCCCGCCTCCGTCTTGCGGGCGAAGGTACAGCGGGCGCGCATGGCGAACCGGCAgggcatcg                                    >  2:44897/1‑245 (MQ=255)                                                                                                                                                                                                                   ccGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTa                                 <  1:306945/251‑1 (MQ=255)                                                                                                                                                                                                                    cGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTAc                                <  1:552910/251‑1 (MQ=255)                                                                                                                                                                                                                            cATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgttcaggt                        <  1:376501/251‑1 (MQ=255)                                                                                                                                                                                                                            cATGGTTCTGGTCAGGTGCAGAGGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgttcaggt                        <  1:484155/251‑1 (MQ=255)                                                                                                                                                                                                                                     ggTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGCCACTGGAAACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTcaggatccggttt                >  2:323391/1‑249 (MQ=255)                                                                                                                                                                                                                                     ggTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTaa               <  1:197023/251‑1 (MQ=255)                                                                                                                                                                                                                                           gTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCa         <  1:637219/251‑1 (MQ=255)                                                                                                                                                                                                                                                  tGCGTCAGGGATTCTTTGCCGTGCAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAACGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCCTCCGTCTTGCGGGCGAAGGGGAAGAGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAACACCCGAtt  >  2:163833/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        ATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTG‑ACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGC‑AGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATT  >  NC_000913/14442‑14932

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