Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,054,830	A→C	16.2%	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.2%	150.4 / 4.7	73	T396T (ACA→ACC)	shiA	shikimate transporter
Reads supporting (aligned to +/- strand):  ref base A (39/19);  new base C (2/12);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 5.34e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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