Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,054,830	A→C	16.5%	T396T (ACA→ACC)	shiA →	shikimate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,054,830	0	A	C	16.5%	149.8 / 4.9	72	T396T (ACA→ACC)	shiA	shikimate transporter
Reads supporting (aligned to +/- strand):  ref base A (38/19);  new base C (2/12);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 5.89e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTACGTTAATCGGAACGTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACA  >  NC_000913/2054604‑2055057                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     caGGTACGTTAATCGGAACGTTGAGCGCATCTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATACTGGCAAACATTGCGCATGCCATGGGGGTGTGTGTGCAACAACAGATGTTTACCGAAATGTTTGGTGCCAGTTATGGCTATAGTGGCGCTGGAGTCGGTTATAAAGTTGCCAGGTTGGTTGGGCGGGGGTTTACACATtttttttca                                                                                                                                                                                                                          >  2:27222/1‑236 (MQ=255)  aGGTACGTTAATCGGAACGTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATATTGGCGCTGGAGTCGGTTTTCAGGTTGCCAGTTTGGTTGGCGTGGGATTTACACCCTTTTTTGc‑c                                                                                                                                                                                                                         >  1:288709/1‑238 (MQ=255)     tACGTTAATCGGAACGTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑Cg                                                                                                                                                                                                                        <  1:80835/236‑1 (MQ=255)      tcGTTAATCGGAACGTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGc                                                                                                                                                                                                                    <  1:721184/238‑1 (MQ=255)              cGGAACGTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGGGTGGTTGGCGGGGGATTTACACCTTTTATTGC‑AGCTGCACtcttc                                                                                                                                                                                                             >  1:381723/1‑238 (MQ=255)              cGGAACGTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGGGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCGGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCCCtcatc                                                                                                                                                                                                             >  1:675694/1‑238 (MQ=255)               ggAACGTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGGTTTACCCCTTTTATTTC‑CGCTGCACtcata                                                                                                                                                                                                             >  2:169887/1‑236 (MQ=255)                  aCGTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGATGCCAGTTATCGCTATAGGCGCGCTGGAGTCGGGTATCAGGTTGCAATTGTTGTTGGCGGTGGTTTTAACCCTTTTTTTGC‑CGCTGCCCtcatacact                                                                                                                                                                                                         >  2:249481/1‑233 (MQ=255)                                         tATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCATGGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATacg                                                                                                                                                                                  >  1:662297/1‑236 (MQ=255)                                                           aTCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACCTGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTCTAGTGACGCAGGAGTCGGTTATCAGTTTGCCCGTGTGGTTGGCGGTGGGTTTACACCCTTTATTGC‑CGCTGCAACCATCACTTACTTTGCCGGGACCTGCCATAGCGTTCCCatttatttcccgg                                                                                                                                                               >  2:129634/1‑234 (MQ=255)                                                                    cTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATg                                                                                                                                                        >  1:689859/1‑237 (MQ=255)                                                                       gATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCGGa                                                                                                                                                   >  1:10556/1‑239 (MQ=255)                                                                        aTAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGAt                                                                                                                                                  >  1:688354/1‑239 (MQ=255)                                                                                 ttctCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAACTGTTTGGTGCAAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTGTTGC‑TGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGAGTGGCCATTTCTTTGCTGGCTGTATGCCTGATTTCCGCAAt                                                                                                                                         >  2:116972/1‑239 (MQ=255)                                                                                                gCACACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCCCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTTCACCTTTTATTGC‑CCCTGCACTCATCCCTTACTTTGCCGGGAACTGGCTAAGCGTCGCAATTTATTTATTGGGTGGATGCCTGAATTCCCAAATGACCCCTttggtgga                                                                                                                          >  2:302177/1‑237 (MQ=255)                                                                                                gCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTttgttgtt                                                                                                                          >  1:849586/1‑237 (MQ=255)                                                                                                gCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATACGGTCGCCATTTATTTGCTTGCTGGGTGCCTTATTTCCCCAATGACCCCtttttgttt                                                                                                                          >  2:580989/1‑233 (MQ=255)                                                                                                    aCATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGaaa                                                                                                                      <  1:309261/239‑1 (MQ=255)                                                                                                       ttGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCCTCACTTACTTTGCCGGGAACTGGCATAACGTCGCCGTTTATTTGCTGCCTGGATGCCTGCTTTCCGCCATGACCGCTTTGTTGATGAAGGAc                                                                                                                   >  1:143490/1‑239 (MQ=255)                                                                                                        tGCGCATGACATGGTGGTGTGGGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCCTCACTTACTTTGCCGGGAATTGGCATAGCCTCGCCCTTTTTTTTCTGGCTGTATGCCTGCTTTCCCCATTGACCGCTtttttgcgtaaagaaa                                                                                                                  >  2:623177/1‑229 (MQ=255)                                                                                                               gACATGGTTGTGTGTGTGCAATAACCGATGTTTACCGAAATGTGTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTTCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAAcg                                                                                                           <  1:464358/239‑1 (MQ=255)                                                                                                                atgtggtggTGTGGTGGCAACAACCGATGTTTACCGAAATTTTTGGGGCCAGTTTTCGGTATAGTGGCGCTGGGGTTGGTTATCAGGTTGCCAGTGTGTTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAAcgc                                                                                                          <  2:683974/236‑1 (MQ=255)                                                                                                                    ggtggtGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGTTGAAAGACAGTCAACGCGCtt                                                                                                      >  1:649054/1‑239 (MQ=255)                                                                                                                      tggtgTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGCTGCCTGATTTCCGCAAGGACCGCTTTGTTGTTGAAAACAAGTCAACAGGCTTGa                                                                                                    >  2:287788/1‑239 (MQ=255)                                                                                                                       gttgtttgtgTAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATTGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGa                                                                                                    <  2:689859/236‑1 (MQ=255)                                                                                                                         tgtgtgtgCCACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCCTCACTTACTTTGCCGGGAACATGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCCCAATGACCGCTTTGTTGTTGAAAGAAAGGCaaacggcttaaaa                                                                                                  >  2:284815/1‑227 (MQ=255)                                                                                                                         tgtgtgtgCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGA‑TTACACCTTTTATTGC‑CGCTGGACTCATCACTTACTTTGCCGGGAACTGGCCTAGCGTCGCACTTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGCCCCCTTTGTTGTTGAAAAAAAGTCAACGGGCtttaaag                                                                                                 >  2:412280/1‑233 (MQ=255)                                                                                                                           tgtgtgCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTTGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCTTCGCCATTTATTTGTCGGCCTGATGCCTGTTTTCGCCACTGACCGCTTTGTTGATGAAAGAAAGTCAAAGCCCTTgcaagca                                                                                               >  2:220416/1‑233 (MQ=255)                                                                                                                           tgtgtgCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGGTAg                                                                                                 >  1:496977/1‑237 (MQ=255)                                                                                                                           tgtgtgCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAg                                                                                                 >  1:139419/1‑237 (MQ=255)                                                                                                                           tgtgtgCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGcc                                                                                               >  2:98780/1‑239 (MQ=255)                                                                                                                           tgtgtgCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGAGGATTTAAACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGACCTGGCATAGCGTCCCCATTTATTTGCTGGCTGGATGCCGGATTTCCGCATTGACCGCTTTGTTGATGAAAACAAGCAAACGCGCtggtaaacc                                                                                               >  2:628368/1‑231 (MQ=255)                                                                                                                           tgtgtgCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGATGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGCCAGTCAACGAGCTTGATAccc                                                                                               >  1:817446/1‑239 (MQ=255)                                                                                                                               tgCAATAACCGATGTTTACCGAAATGTTTGGTTCCAGTTTTCGCTATAGTGGCGCTGGCGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTTGATTTACACTTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGc                                                                                           <  1:803351/239‑1 (MQ=255)                                                                                                                                     aacCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTTCTGGCTGGATGCCTGATTTCCCAAATGACCCCTTTGTTGTTGAAAGCCAGTCAACGCGCTTGATAACCTgggaaaaaa                                                                                     >  2:142544/1‑232 (MQ=255)                                                                                                                                        cGATGTTGACCGAAATGGTTGGGTCCAGTTATCGCTATAGGGGCGCTTGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACCCCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGatcat                                                                                  <  1:249481/239‑1 (MQ=255)                                                                                                                                            gTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGCTGAAAGAAAGTCAACGCGCTTGGTAGCCTGTGGAAAATAAGCCGa                                                                              >  2:777013/1‑239 (MQ=255)                                                                                                                                               ttccGAAATTTTTGGTGCCAGTTTTTGCTATGGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCt                                                                           <  2:849586/237‑1 (MQ=255)                                                                                                                                                cctggaaTGTTTGGTGCCAGTTATCGTTCTAGTGGGCTTGGAGTCGTTTATCAGGTTGCCCTTTTGGTTGGCGGTGGATTTACCCCTTTTATTGC‑CGCGTCACTCATCACTTACTTTGCCGGGAACTGGCATAGCTGCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCtt                                                                          <  2:472481/234‑1 (MQ=255)                                                                                                                                                         tttGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCa                                                                                                              >  1:829764/1‑194 (MQ=255)                                                                                                                                                         tttGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCa                                                                                                              <  2:829764/194‑1 (MQ=255)                                                                                                                                                                     tATCGCTATGGTGGCGCTGGAGTCGGTTACCAGGTTGCCAGTGTGGTTGGCTGTGGATGTACACCTTTTATTGC‑CTCTGCACTTATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTcac                                                     <  1:27222/239‑1 (MQ=255)                                                                                                                                                                         gcttttgTGGCCGTGGAGTCGGTTATCAGGTTGCCAGTTGGGTTGGCGGTGGATTTACCCCTTTTATGGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACAttt                                                 <  2:649054/233‑1 (MQ=255)                                                                                                                                                                               gTGGCGCTGGAGTCGGTTAGCAGGTTGCCAGTGTGGTTGGCGGTGAATTTACGCCTTTTATTGC‑CGCTGCACTCATCGCTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGa                                           <  1:284815/239‑1 (MQ=255)                                                                                                                                                                                 ggCGGTGGATTCGTTTATCAGGTTGCCAGTGTGGTTGGGTGTGGATTTGCACCTTTTAATGC‑CGCTGCACTCATCACTTACGTTGCCGGGAACTGGCATAGCGGCGCCATTTATTTGATGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCTATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACa                                         <  1:410835/239‑1 (MQ=255)                                                                                                                                                                                 ggCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGAc                                          <  1:206163/238‑1 (MQ=255)                                                                                                                                                                                 ggCGCTGGAGTCGGTTATCAGGGTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGAc                                          <  1:577584/238‑1 (MQ=255)                                                                                                                                                                                  gcgcTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTTGATGCCTGATTTCCGCCATGACCGCTTTGTTTATGGAAGAAAGTCAACGCGCTTGATAGCCTGGCGCAGATACACCGATCTTCACATTAAACCTTTTTTTTACATTtgtggtaaaa                                        >  2:124776/1‑233 (MQ=255)                                                                                                                                                                                  gcgcTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGCTGAAAGACAGTCAACGCGCTTGCTAGCCTGGCGAAAATCCTCCGATCTTCGCCTTACACTTTTGTTTTACATTTCTGTGCCAt                                        >  2:225025/1‑239 (MQ=255)                                                                                                                                                                                  gcgcTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACAt                                        >  1:570310/1‑239 (MQ=255)                                                                                                                                                                                   cgcTGGAGGCGGTTATCAGGTTGCCAGTGGGGTTGGCGGTGGATTGACACCTTTTAGTGC‑CGCTGCCCTCAGCACTTACTTGGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACAt                                        <  1:113631/238‑1 (MQ=255)                                                                                                                                                                                    gcTGGAGACGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCCTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCCTACCCATTTGTTTTAAATTTTTGTGACaaaa                                      >  2:155024/1‑236 (MQ=255)                                                                                                                                                                                     cTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTTTTGATGAAAGACAGTCAACCCGCTTGATAGCCTGGCGAGAATCCTCCGCTCTTCGCCTTACACTTTTTTTTCCCATTTCGTTGAcctccc                                     >  2:113631/1‑234 (MQ=255)                                                                                                                                                                                     cTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATCCt                                     >  1:102279/1‑239 (MQ=255)                                                                                                                                                                                     cTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTAGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCACCGCGCATGATAGCCTGGCGGAGATCATCCGATCTTCGCCTTcca                                                                 >  1:522638/1‑211 (MQ=255)                                                                                                                                                                                          gTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGGTCTTCGCCTTACACTTTTGTTTCACATTTCTTTGAAATAATATCgg                                >  2:206163/1‑239 (MQ=255)                                                                                                                                                                                             ttttttCAGGTTGCCCGTTTGGTTGGCGGTGGGTTTACCCCTTTTATGTC‑CGCTTCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGAtg                             <  2:381723/234‑1 (MQ=255)                                                                                                                                                                                             tgttttttggTTGGCCGTGTGGTTGGGTGTGGAATTCCACCTTTT‑TTGCGCGCTGCACTCATCACTTACTGTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGAtg                             <  2:484175/231‑1 (MQ=255)                                                                                                                                                                                             ggTTTTCAGGTTGCCAGTGTTGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACGGGCATAGCGTCGCCATTTAATTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGAtg                             <  1:658508/239‑1 (MQ=255)                                                                                                                                                                                                tttCAGGTTGCCGGTTGGGTTGGCGGTGGATTTACCCCTTTTATGGC‑CGCTGCACTCATCCCTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGc                          <  2:496977/237‑1 (MQ=255)                                                                                                                                                                                                tatacggTTTCCAGTGGTGTTGGCCGTGGATTTACCCCTTTTAATGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCAGCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGc                          <  1:129634/234‑1 (MQ=255)                                                                                                                                                                                                         gCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGGAAGGCAGTCAACGCGCTTGATACCCTGGCGAGAGTCATCCGTTCTTCGCCTTACCCTTTTTTTTTACATTTTTTTGACATACTACCGGATGTTCGGCAATTga                 >  2:69178/1‑238 (MQ=255)                                                                                                                                                                                                         gCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTACCGGATGTGCGGTAATTGt                 >  1:880877/1‑239 (MQ=255)                                                                                                                                                                                                          ccAGTGTGGTGGGCGGTGGATTTTCCCCTTTTATGGC‑CGCTGCAGTCATCAATTGCTTTGCCGGGAACTGGCATAGCGTCGTCATTTATGTGCTGGCTGGATGCCTGATTTCCGCAATGACCGGTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTa                <  2:662297/239‑1 (MQ=255)                                                                                                                                                                                                           cAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTa                <  1:98780/238‑1 (MQ=255)                                                                                                                                                                                                              tttggtggGGGGTGTATTTTCCCTTTTTATTGC‑CGCTTCACTCACCACTTACTTGGCCGGGAACTGGCATCGCCTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGCCCGCTTTGTTGATGAAAGAAAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGa            <  2:160995/237‑1 (MQ=255)                                                                                                                                                                                                               gtgGTTGGCGGTGGATTTACCCCTTTTATTGC‑CGCTGCACTCATCCCTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATg              <  1:623177/236‑1 (MQ=255)                                                                                                                                                                                                                 ggTTGGCGGTGGATTTACCCCTTTTATTTC‑CGCTGCAATCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAAc          <  1:567307/238‑1 (MQ=255)                                                                                                                                                                                                                      gCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTTATCCCTTACTTGGCCGTGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGgaga    <  1:368220/239‑1 (MQ=255)                                                                                                                                                                                                                       cGGTGGATTTACCACTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCAACCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGgaga    <  1:287788/238‑1 (MQ=255)                                                                                                                                                                                                                        ggTGGA‑TTACACCTTTTATTGC‑CGCTTCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTGGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGaca  <  1:412280/238‑1 (MQ=255)                                                                                                                                                                                                                         gTGGATTTACCCTTTTTTTTGC‑CGCTGCATTCATCACTTTTTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGaca  <  2:288709/238‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     CAGGTACGTTAATCGGAACGTTGAGCGCATTTCCTTTCTTTATGGCGCTTGAAGCACAATCTATTTTCTGGATAGTTTTCTTCTCCATAATGCTGGCAAACATTGCGCATGACATGGTGGTGTGTGTGCAACAACCGATGTTTACCGAAATGTTTGGTGCCAGTTATCGCTATAGTGGCGCTGGAGTCGGTTATCAGGTTGCCAGTGTGGTTGGCGGTGGATTTACACCTTTTATTGC‑CGCTGCACTCATCACTTACTTTGCCGGGAACTGGCATAGCGTCGCCATTTATTTGCTGGCTGGATGCCTGATTTCCGCAATGACCGCTTTGTTGATGAAAGACAGTCAACGCGCTTGATAGCCTGGCGAAGATCATCCGATCTTCGCCTTACACTTTTGTTTCACATTTCTGTGACATACTATCGGATGTGCGGTAATTGTATGGAACAGGAGACA  >  NC_000913/2054604‑2055057

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