Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,336,353	C→T	P175S (CCA→TCA)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,336,353	0	C	T	92.5%	147.3 / 1.1	54	P175S (CCA→TCA)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (26/23);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTTGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTA  >  NC_000913/1336119‑1336593                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 gttgttGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTg                                                                                                                                                                                                                                  <  1:9009/251‑1 (MQ=255)            ggcggcAACGCGCGAAGCGGTTAAACGCCCCGGCTGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGAGCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAg                                                                                                                                                                                                                           <  2:293049/247‑1 (MQ=255)              cggcAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAAc                                                                                                                                                                                                                     <  1:23624/251‑1 (MQ=255)                 cAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAACGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTc                                                                                                                                                                                                                  >  1:513860/1‑251 (MQ=255)                  aaTGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTTTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCg                                                                                                                                                                                                                 >  1:275929/1‑251 (MQ=255)                         gAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATAtct                                                                                                                                                                                                          >  1:489268/1‑251 (MQ=255)                          aaGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATAtctc                                                                                                                                                                                                         <  1:205652/251‑1 (MQ=255)                           aGCGGTTAAACGCCTCGGCGGCGATACGGCAATGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGAGCGGAGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCATGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTATCCTCGAATATCTCg                                                                                                                                                                                                        <  2:553784/251‑1 (MQ=255)                                   aaCGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCAACGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCa                                                                                                                                                                                                <  1:255699/251‑1 (MQ=255)                                           ggcggcGATACTGCAAAGGTTAACCCGCTCTCACCGGACGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCATGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGCATATCTCGGCAAAGCAGTGTgg                                                                                                                                                                                          >  1:18683/1‑249 (MQ=255)                                             cggcgATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGACTATCTAGGCAAAGAAGGGGGGAGTg                                                                                                                                                                                      >  2:107641/1‑251 (MQ=255)                                              ggcgATACTGCAAAGGTTAACCCGCTCTCACCTGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGCTGATGAGGCCTTTGCAGAGAACGTACGCCTGGAAAGGGAGCGCAACCACGAACGTTATGTGTACCTGAAATGGGGAGAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACCGGCATTGGCCTTCAGGTACACCTCGAACATCTCGGCAAAGAAGTGTGGGGTg                                                                                                                                                                                      >  2:56318/1‑250 (MQ=255)                                                 gATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAAc                                                                                                                                                                                                                     >  1:241464/1‑216 (MQ=255)                                                 gATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAAc                                                                                                                                                                                                                     <  2:241464/216‑1 (MQ=255)                                                           aGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGt                                                                                                                                                                        <  1:1392/251‑1 (MQ=255)                                                             gTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTg                                                                                                                                                                       >  1:518738/1‑250 (MQ=255)                                                                     gCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGAACGTCGATCGTTTTGGTGATGGTGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATAGCAGGACGGTGAATGGATTg                                                                                                                                                              >  1:509585/1‑251 (MQ=255)                                                                                     cTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACt                                                                                                                                              >  1:437343/1‑251 (MQ=255)                                                                                       ggTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCg                                                                                                                                            >  1:517439/1‑251 (MQ=255)                                                                                          cATTGACCACTCGGTGAACGTCGATCGTTTTGGCGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGCCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTg                                                                                                                                         >  1:461945/1‑251 (MQ=255)                                                                                           aTTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCATGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTgg                                                                                                                                        >  1:115789/1‑251 (MQ=255)                                                                                             tGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGt                                                                                                                                       >  1:473580/1‑250 (MQ=255)                                                                                              gACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATTGATTGCTTATCCGGATACACTCGTTGGTa                                                                                                                                      <  1:193940/250‑1 (MQ=255)                                                                                               aCCACTCGGTGACCGACGATCGTTTTGGTGATGATGAGGCATTTGACGAAAACGTACGCCTGGAAATGGAGAGCAACCACGAACGTTATGGGTTCCTGAAATGGGGAAAGCAAGCGTTCAGCCGGTTTAGCGTCGTGCCGTCAGGGACAGGCATTTGCCATAAGGTTAACCTCGAATATCCCGGCAAAACAGTGTGGAGTGAATTGCCGGACGGTGAAGGGATTGCTTATCCggtaacacacggtggta                                                                                                                                      >  2:374034/1‑234 (MQ=255)                                                                                                  ccTCGGTGTCCGGCGATCGTTTTGGGGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGAc                                                                                                                                 <  2:39050/250‑1 (MQ=255)                                                                                                     cGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGAATCg                                                                                                                              >  2:541089/1‑251 (MQ=255)                                                                                                      ggTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGACAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCGCGGAAAAGCAGTGGGGAGTGAATTGCAGGACGGGGAATGGATTGCTTATCCGGATACACTCGTTGGTATTGCCTCGc                                                                                                                             >  2:526923/1‑251 (MQ=255)                                                                                                          aCCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCacac                                                                                                                         <  2:518738/251‑1 (MQ=255)                                                                                                              tCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCATGAAATAGGGAAAGCAAGCGTTCAGTAGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACaccacc                                                                                                                     >  1:29209/1‑251 (MQ=255)                                                                                                               cGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACaccacc                                                                                                                     <  1:358746/250‑1 (MQ=255)                                                                                                                gATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCACGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTCCTGAATAGAACaccaccct                                                                                                                   >  2:196578/1‑249 (MQ=255)                                                                                                                gATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCAt                                                                                                                   <  1:546435/251‑1 (MQ=255)                                                                                                                  tCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATg                                                                                                                  >  1:407555/1‑250 (MQ=255)                                                                                                                      tttGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCa                                                                                                              >  1:434675/1‑250 (MQ=255)                                                                                                                                              gaaAACGTCCGCCTGGAAATGGAGAGCAACCACGAACGTTATGTGTTAATGAAATGGGGAAAGCAAGCGTTCAGACGGTTTAGCGTCGAGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGACTATCTCGGCAAAGCAGTGTGGAATGAATTGCAGACCGGTGAATGGATTGCTTATCCGGAAACACGCATTGGAACTGACTCGCACACCACCAGTAACAACGGCCTTAGCGTGAaggggtgggg                                                                                     >  2:434579/1‑251 (MQ=255)                                                                                                                                              gaaAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCtggggtgggg                                                                                     <  1:56318/251‑1 (MQ=255)                                                                                                                                                                  ggAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGc                                                                   >  1:121197/1‑249 (MQ=255)                                                                                                                                                                  ggAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCa                                                                  <  2:283138/250‑1 (MQ=255)                                                                                                                                                                   gAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGa                                                                <  1:541089/251‑1 (MQ=255)                                                                                                                                                                       gcAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGCTCGAAGCAGAAGcc                                                            >  1:562293/1‑251 (MQ=255)                                                                                                                                                                         aaCCACGAACGTGATGGGATACTGAAATGGGGAAAGCAAGCGTGCAGTCGGTTTAGCGGCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAATCAGGGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGTATACACGCGTTGGTACTGTCGCGCACACCACCATGACCAACGGCCTTGGCGTGCGGGGGTGGGGCGTTGGTGGGATCGAAGCATAAGCCGc                                                          <  2:513860/251‑1 (MQ=255)                                                                                                                                                                          aCCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCa                                                         <  1:158040/251‑1 (MQ=255)                                                                                                                                                                              cGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGAGTGGAGTGAATTGCAGGCCGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGACGTGCTGGGGTGGGGCGTTGGTGGGATCGGAGCAGAAGCCGCCATGt                                                     >  2:295955/1‑251 (MQ=255)                                                                                                                                                                                             ccTGAAATGGGGAAAGCAAGTGTTCATTCGGTTTAGCGTCGGGCCGTCAGGCACAGGCATTTGCCAGCAGGTTAACCGCGAATATCTCGGCAAAGCAGTTTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCTGATACACTCGTTGGTACTGACTCGCACCCCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGTATCGAAGCAGAAGCCGCAATTTTAGGCCAGCCGGttt                                      <  2:560185/251‑1 (MQ=255)                                                                                                                                                                                               tGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTcc                                    >  1:543075/1‑251 (MQ=255)                                                                                                                                                                                                           agcaagtGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCAGCAGGTTATCCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACTGCCTTGGCGTGCTGGGGTGGGGCGTTGGAGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCg                        <  2:489268/251‑1 (MQ=255)                                                                                                                                                                                                             caagcGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGa                                                                                                                                                                            >  2:47164/1‑101 (MQ=255)                                                                                                                                                                                                             caagcGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGa                                                                                                                                                                            <  1:47164/101‑1 (MQ=255)                                                                                                                                                                                                                 cGTTCAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTTCCATCAGGTTAACATCGACTATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTCTCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCACCGGCCTTGGAGTGCTGGGGTGGGGCGATGGTGGGATCGAAGCAGACGCCGCAATGTTAGGCCGGCCGGTTTCCATGCTTATCCCGGATGtg                  >  2:540729/1‑250 (MQ=255)                                                                                                                                                                                                                     cAGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTgg              >  1:341756/1‑251 (MQ=255)                                                                                                                                                                                                                      aGTCGGTTTAGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTggg             <  1:196578/251‑1 (MQ=255)                                                                                                                                                                                                                               aGCGTCGTGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACt    <  1:295955/251‑1 (MQ=255)                                                                                                                                                                                                                                 cgtcgtGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCtggg                                                                                           <  1:165191/162‑1 (MQ=255)                                                                                                                                                                                                                                 cgtcgtGCCGTCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTa  <  1:24930/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 GTTGTTGATCTGGCGGCAATGCGCGAAGCGGTTAAACGCCTCGGCGGCGATACTGCAAAGGTTAACCCGCTCTCACCGGTCGACCTGGTCATTGACCACTCGGTGACCGTCGATCGTTTTGGTGATGATGAGGCATTTGAAGAAAACGTACGCCTGGAAATGGAGCGCAACCACGAACGTTATGTGTTCCTGAAATGGGGAAAGCAAGCGTTCAGTCGGTTTAGCGTCGTGCCGCCAGGCACAGGCATTTGCCATCAGGTTAACCTCGAATATCTCGGCAAAGCAGTGTGGAGTGAATTGCAGGACGGTGAATGGATTGCTTATCCGGATACACTCGTTGGTACTGACTCGCACACCACCATGATCAACGGCCTTGGCGTGCTGGGGTGGGGCGTTGGTGGGATCGAAGCAGAAGCCGCAATGTTAGGCCAGCCGGTTTCCATGCTTATCCCGGATGTAGTGGGCTTCAAACTTA  >  NC_000913/1336119‑1336593

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