Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	17,893	T→C	A135A (GCT→GCC)	nhaA →	Na(+):H(+) antiporter NhaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	17,893	0	T	C	98.0%	150.6 / ‑4.4	51	A135A (GCT→GCC)	nhaA	Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (26/24);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

GAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGA  >  NC_000913/17648‑18137                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     gAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTATATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGc                                                                                                                                                                                                                                                 >  1:73513/1‑251 (MQ=255)                       tATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGtt                                                                                                                                                                                                                            <  1:205449/250‑1 (MQ=255)                            aTAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTATGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCCGTTATCGCCGCTATTGGTGGGATGATTGTTCAGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGCTTACCCGCGAAAGGTGGGGGATCCCGGCGGCTACTGCCATTGCCTTTGCACTTTGTGTGCTGGCGCTGTTgggga                                                                                                                                                                                                                       >  2:423580/1‑248 (MQ=255)                               aaTGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACCTTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCg                                                                                                                                                                                                                   >  1:626040/1‑251 (MQ=255)                               aaTGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCg                                                                                                                                                                                                                   >  1:258334/1‑251 (MQ=255)                                               cGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTg                                                                                                                                                                                                   >  1:135105/1‑251 (MQ=255)                                                        tCCTGTTAGTCGGTCGTGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCtt                                                                                                                                                                                           <  1:25090/250‑1 (MQ=255)                                                            gTTAGGCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTg                                                                                                                                                                                       <  1:605847/250‑1 (MQ=255)                                                             ttAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACACGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCCCTTGGGGGGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAACGCTGTAGCTCTTTTTGAt                                                                                                                                                                                     >  2:225756/1‑251 (MQ=255)                                                               agtacgTCTGAAAGTGAAACTTGAATTCATGCAAGGAGTGCTAGCCTGCTTACTCCAGGCCGCATTTCCAGTTAGCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATgg                                                                                                                                                                                   <  2:472212/246‑1 (MQ=255)                                                                  cGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGctc                                                                                                                                                                                <  1:8056/251‑1 (MQ=255)                                                                        ggAAGTTAAACGTGAACTGATGCAAGGTTCGCTAGCCAGCTTGCGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGACATTTCTCTATCTGGGCTTTAACTATGCCGCGCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGTCTTTGCACTTGGTGTCCTGGCGCTGGTGGGAAGGCGTGTTTCGTTAGCGCTGAAGATCTTTTTAATGGCCCTGGCta                                                                                                                                                                          <  2:238315/251‑1 (MQ=255)                                                                           aGTTAAACGTGAACTGTTGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTTTCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGTGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTTTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCtatta                                                                                                                                                                       <  1:562080/251‑1 (MQ=255)                                                                            gTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCtattat                                                                                                                                                                      >  1:9131/1‑251 (MQ=255)                                                                             ttAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATc                                                                                                                                                                     >  1:416352/1‑251 (MQ=255)                                                                              gtaccgtgaaatgtatgaaggttcccTAGCCAGCTTACTCCAGGCCACATTTCCATTGATCGTCGCTATTGGTGTGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGAGTCGATTACCCGCGTAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCCCGTGGTGTGCTGTCGCTTTTGTTAAGTCGTGTTCCGTTAGCTCTGAAGATCTTTTTGATGGCGCTGGCTAGTATcg                                                                                                                                                                    <  2:631843/226‑1 (MQ=255)                                                                                    tGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATc                                                                                                                                                              >  1:427146/1‑251 (MQ=255)                                                                                             gCAAGGATCGCTAGCCAGCTTACGCCTAGCCGCATTTCCAGTTTCCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTGGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGGGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCca                                                                                                                                                     <  2:469232/251‑1 (MQ=255)                                                                                              cAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCcat                                                                                                                                                    <  1:647145/251‑1 (MQ=255)                                                                                                    tCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGCTCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGGGTTCCGTTAGCGCTGAAGGTCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCCTTat                                                                                                                                              >  2:292247/1‑251 (MQ=255)                                                                                                       cTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCCTTatcat                                                                                                                                           >  1:201437/1‑251 (MQ=255)                                                                                                           ccAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGc                                                                                                                                        >  1:215764/1‑250 (MQ=255)                                                                                                               cTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGt                                                                                                                                   >  1:333330/1‑251 (MQ=255)                                                                                                                  aCGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCt                                                                                                                                <  1:423580/251‑1 (MQ=255)                                                                                                                         gCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTaa                                                                                                                         <  1:85317/251‑1 (MQ=255)                                                                                                                                     gTTATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGTCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGAt                                                                                                             <  2:385949/251‑1 (MQ=255)                                                                                                                                     gTTATCGCCGCTAGTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACGTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGAt                                                                                                             <  1:611336/251‑1 (MQ=255)                                                                                                                                       tATCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATgg                                                                                                           >  1:93461/1‑251 (MQ=255)                                                                                                                                        aTCGCCGCTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATgg                                                                                                           <  1:603276/250‑1 (MQ=255)                                                                                                                                             cgcTATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTTATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCctctc                                                                                                     <  1:290283/251‑1 (MQ=255)                                                                                                                                                tATTGGTGGGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTg                                                                                                  >  1:624202/1‑251 (MQ=255)                                                                                                                                                    ggttggATGATTGTTCCTGCATTTCTCTGTCCGGCCTTTAAGTATGCCGATCCGATTACCCGTTGAGGGTGGGCTATCCCGGCTGCTACTGAGATTGCCTTTGCACTTGGTCTGCTGGCTCTGTTGGGAAGTCGTGTTCCGTTATCGCTGAAGATCTTTTTGATGGCTCTGGCTATTACCTACGATCTTGGGGCCATCATGATCATCGAATTGTTCTACAATAATGACTTATCGATGGCCTCTCTTGGCg                                                                                               <  2:176102/246‑1 (MQ=255)                                                                                                                                                      tggGATGATTGTTCCGGCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTcg                                                                                            >  1:199864/1‑251 (MQ=255)                                                                                                                                                                      gCATTACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCTCTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGCTCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGTCGGCCTCTCTTTGCGTCCCGGCTGTAGCCATTg                                                                             >  2:66732/1‑250 (MQ=255)                                                                                                                                                                          tACTCTATCTGGCCTTTAACTATGCCGATCCGATTGCCCGCGAAGGGTGGGCGAGCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTg                                                                             <  2:198983/246‑1 (MQ=255)                                                                                                                                                                          tACTCTATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTg                                                                             >  1:198983/1‑246 (MQ=255)                                                                                                                                                                            ctctATCTGGCCTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACt                                                                      >  1:324926/1‑251 (MQ=255)                                                                                                                                                                                atctgtcctttAACTAGGCCGATCCGATTTCCCGCGAAGGGGGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAGTTGCGGTACTcg                                                                    <  2:9131/241‑1 (MQ=255)                                                                                                                                                                                         ttAACTATGCCGATCCGGTTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACAGTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATGGCGGTACTCGCGGTATTGa                                                           <  2:93461/249‑1 (MQ=255)                                                                                                                                                                                         ttAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCTCTTGGTGTGCTGGCGCTGTTTGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGa                                                           <  1:61488/249‑1 (MQ=255)                                                                                                                                                                                          tAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGCCGATCTTTTTGATGGCTCTGGCTATTATCGCCGATCTTGGGGCCATCCTTTTCATCGCATTGTTCTCCACTACTGACTTATCGTTGGCATCTCTTGGCGTCGCGGCTGTATCAATTGCGGTAATCGCGGTTTTGAATc                                                        >  1:624084/1‑251 (MQ=255)                                                                                                                                                                                            aCTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATGGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCtg                                                      >  1:49227/1‑251 (MQ=255)                                                                                                                                                                                             cTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCtgt                                                     >  1:461422/1‑251 (MQ=255)                                                                                                                                                                                                        cgatTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGCCTTATCGATGGCCTCTCTTGGCGTCGCGGCTGGAGCAATTGCGGTACTGGCGGTATTGATTCTGTGTGGTGTAcgc                                          >  2:28521/1‑251 (MQ=255)                                                                                                                                                                                                              cccGCGAAGGGTGGGCGATCCCGGCGGCTACTGACCTTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTCTTGAATCTGTGTGGTGTACGCCGCACg                                    >  1:162933/1‑251 (MQ=255)                                                                                                                                                                                                                           ggcgaatccGGCTGCTTCTTACATTGCCTTTGCACTTGGTGTGCTGGCGCTGGTGGGAATTCTGGTGCCGTTAGCGCTGATGATCTTTTTGATGGCTCTGGCTATTATCGTCGATCTTGGTGCCATCATTATCATCGCATTGTTCTACACTAATGTCTTATCGATGGCCTCTCTTGGCGTCGCGGCTGGAGCAATTGCGGTACTCGCGGTATTGAATCTGTGCGGTGTTCGCCGCACGGTCGTCtatacta                       <  2:202592/244‑4 (MQ=255)                                                                                                                                                                                                                               aTCCCGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATAGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCtggt                   <  1:416918/251‑4 (MQ=255)                                                                                                                                                                                                                                  ccGGCGGCTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGACGCGGCTGTAGCAATTGCGGTACTCGGGGGATTGGATTTGTGTGGTGTACGCCGCCCGGGCGGCTTTTTTCTGGTTgg                >  2:605916/1‑251 (MQ=255)                                                                                                                                                                                                                                         cTACTGACATTGCCTTTGCACTTGGTGTGCTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCgtagtg         >  1:297171/1‑251 (MQ=255)                                                                                                                                                                                                                                             tGACATTGCCTATGCACTTGGCGTGCTGGCGCTGTTGGGACGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATAATTATCATCGCCTTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTAACCGCGGTATTGAATCTTTGTGGGGTACGCCGCACGGGCGTATATATTCTGGTTGGGGTTGtgttgg     >  2:206444/1‑250 (MQ=255)                                                                                                                                                                                                                                                  ttGCCTTTGCACTTGGTGTGCTGTCGCTGTGTGGAAGTCGTGTTCCTTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTTTTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACTCTAATGACTTAGCGTTGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTTCGGTACTCTCGGTATTGAATCTGTGTGGTGTACGCCGCACGGGCGTCTATATTCTGGTTGGCGTAGTGTTGTGGa  <  2:487534/249‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     GAAATCAACAAAAACATGCTGTTATGGATAAATGACGCGCTGATGGCGGTATTTTTCCTGTTAGTCGGTCTGGAAGTTAAACGTGAACTGATGCAAGGATCGCTAGCCAGCTTACGCCAGGCCGCATTTCCAGTTATCGCCGCTATTGGTGGGATGATTGTGCCGGCATTACTCTATCTGGCTTTTAACTATGCCGATCCGATTACCCGCGAAGGGTGGGCGATCCCGGCGGCTACTGACATTGCTTTTGCACTTGGTGTACTGGCGCTGTTGGGAAGTCGTGTTCCGTTAGCGCTGAAGATCTTTTTGATGGCTCTGGCTATTATCGACGATCTTGGGGCCATCATTATCATCGCATTGTTCTACACTAATGACTTATCGATGGCCTCTCTTGGCGTCGCGGCTGTAGCAATTGCGGTACTCGCGGTATTGAATCTGTGTGGTGCACGCCGCACGGGCGTCTATATTCTTGTTGGCGTGGTGTTGTGGA  >  NC_000913/17648‑18137

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