Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	288,310	A→C	18.6%	T508P (ACC→CCC)	yagH →	CP4‑6 prophage, putative xylosidase/arabinosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	288,310	0	A	C	18.6%	110.5 / 2.1	67	T508P (ACC→CCC)	yagH	CP4‑6 prophage, putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base A (8/39);  new base C (4/11);  total (14/53)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAGGGTAGAACCATCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTT  >  NC_000913/288074‑288544                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            gACAGGGTAGAACCATCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑c                                                                                                                                                                                                                                           <  2:10562/238‑1 (MQ=255)    aGGGTAGAACCATCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGTAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑cccg                                                                                                                                                                                                                                        >  1:929036/1‑236 (MQ=255)       gTAGAACCATCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGgcgc                                                                                                                                                                                                                                     <  2:1441442/238‑1 (MQ=255)        tAGAACCATCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGGGC‑GGCG‑CG‑TT‑CT‑CCGTcgcgg                                                                                                                                                                                                                                   >  2:595903/1‑238 (MQ=255)        tAGAACCATCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCC‑TCTTT‑CA‑CCgggcgg                                                                                                                                                                                                                                    >  1:577518/1‑235 (MQ=255)        tAGAACCATCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCC‑TC‑TT‑TA‑CCcgcgcgg                                                                                                                                                                                                                                   >  1:1259441/1‑238 (MQ=255)            aCCATCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTtg                                                                                                                                                                                                                                <  2:502229/238‑1 (MQ=255)            aCCATCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGGCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTtg                                                                                                                                                                                                                                <  2:709926/238‑1 (MQ=255)                tCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGGGCGGGCT‑CC‑TT‑CC‑CCGCCGCCTTgggggg                                                                                                                                                                                                                            >  1:293862/1‑239 (MQ=255)                tCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGG‑GGCG‑TC‑TT‑CC‑CCCGCCCGCGTGTgggg                                                                                                                                                                                                                           >  2:414412/1‑238 (MQ=255)                tCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGGC‑TC‑TT‑TA‑CCGgcggggtttttgg                                                                                                                                                                                                                            >  1:598911/1‑228 (MQ=255)                   aaGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTg                                                                                                                                                                                                                        <  2:1308644/239‑1 (MQ=255)                    aGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGc                                                                                                                                                                                                                       <  1:5388/239‑1 (MQ=255)                       tATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCGCGAGCAGCGCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCa                                                                                                                                                                                                                      <  2:405453/237‑1 (MQ=255)                              cTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGGGG‑GGGC‑TC‑TT‑TACCCGCCGGGTTTTTGGGCC‑GCCCTGCggggg                                                                                                                                                                                                             >  1:779530/1‑238 (MQ=255)                               tCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAg                                                                                                                                                                                                               <  1:521641/236‑1 (MQ=255)                                        aCGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTGCTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAgcggc                                                                                                                                                                                                   <  1:174480/239‑1 (MQ=255)                                        aCGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAgcggc                                                                                                                                                                                                   <  1:927746/239‑1 (MQ=255)                                             ccGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGGT‑TC‑TT‑TC‑CCGGCGCGGTTTTTGGGCTGCAcctcgcgggcatgcacggcggag                                                                                                                                                                                               >  1:911048/1‑216 (MQ=255)                                               gtcgtGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGgct                                                                                                                                                                                            <  1:828324/239‑1 (MQ=255)                                                         cTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCTCCC‑TT‑CC‑CCGCGCCGTTTTTGGGCC‑GCCCTGGGGGGGCATCCGCGGGGGCgggcggcagcgcgc                                                                                                                                                                                  >  1:796641/1‑227 (MQ=255)                                                          tGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGGG‑TC‑TT‑CA‑CCGGCGCGCTTGTGGGGCTGCCCCGCGAGGGCACCAACTGGGGCGGGCGCTACGCgg                                                                                                                                                                                   >  1:630546/1‑237 (MQ=255)                                                              cGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTTCCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCg                                                                                                                                                                             <  2:458498/239‑1 (MQ=255)                                                                cgcagcCCATGCCGGTGCCGGAGCGTTCGGTGGGGGTTCGGTAGCGGGTGGACGTGGATGCGCTGGTCTCCCGCTCAAGCTAGTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGGGGCGTGGAAGCTGTCGGACGGCTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGGTTGTGAGGCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCg                                                                                                                                                                             <  2:309551/236‑1 (MQ=255)                                                                 gcagcCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGTAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGac                                                                                                                                                                           <  2:929036/238‑1 (MQ=255)                                                                 gcagcCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGac                                                                                                                                                                           <  1:414412/238‑1 (MQ=255)                                                                 gcagcCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGac                                                                                                                                                                           <  1:664147/238‑1 (MQ=255)                                                                 gcagcCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGac                                                                                                                                                                           <  1:595903/238‑1 (MQ=255)                                                                   agcCTATTCCGGAGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGacta                                                                                                                                                                         <  1:216603/238‑1 (MQ=255)                                                                   agcCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGacta                                                                                                                                                                         <  1:1492778/238‑1 (MQ=255)                                                                      ccATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACtt                                                                                                                                                                      <  1:744393/238‑1 (MQ=255)                                                                                gCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAg                                                                                                                                                            <  1:1442041/238‑1 (MQ=255)                                                                                gCCGGAACAGGGGGAGAGCGTCTGGCGGCGGGTGGACGTGGATAGGCTGGTCTCCCGTTACAGCTACTCGTTTGATGGGGAGACGGGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCGGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAg                                                                                                                                                            <  2:680625/238‑1 (MQ=255)                                                                                 ccGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGGT‑TC‑TT‑CA‑CCCGCGCGCTTGTGGGGCTGCACTGCGGGGGCATCAGCGGGGGCGGCTGCTACGCGGGCTTCGACTACTTCACCTACGAGcc                                                                                                                                                          >  1:1429627/1‑239 (MQ=255)                                                                                  cGGAACATGCGGGGAGCGTCTGGCTGCGGGTGGAGGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGcc                                                                                                                                                          <  1:1309361/238‑1 (MQ=255)                                                                                       cgtGCGGGGCGGGTCGGGCTGGGGGGGGAGGGGGATACGCTGGTCTCCCGTTACAGCTACTCGTTTGAGGGCGAGACGGGGCACCCCGTGCCGGTGACGTAGGAGGCGTGGAAGCTGTCGGACGACTACACCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTCCTTCACCTACGAGCCGGTCTa                                                                                                                                                    <  2:224702/237‑1 (MQ=255)                                                                                       cATGGGGAGGGCGTTGGGCGGCGGGTGGACGGGGTTGCGCGGGTCTCCCGCCACAGCTACTCGTTTGAGGGCGAGACGGGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CC‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTa                                                                                                                                                    <  2:779530/239‑1 (MQ=255)                                                                                             gaggGCGTCGGGCTGCGGGTGGAGGTGGATACGCGGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGG‑GGCT‑TC‑TT‑CA‑CGGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGc                                                                                                                                               <  2:458324/238‑1 (MQ=255)                                                                                                 gggTCTGGCGGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGCCGTGGCCCACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAg                                                                                                                                           <  1:537098/236‑1 (MQ=255)                                                                                                  cgttgggcggcGGGGGGCTGTGGATCCGCTGGTTTACCGCTAAGGCTACGCGTTGGATGGCGGGACTGGGCCCACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTGGGACGATTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CC‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGAAATCAGGGGCGCCGGCTGCTACGCGGACTTCGACTACTTCCCCTCCGAGCCGGTCTAACGGCTCAgg                                                                                                                                          <  2:577518/233‑1 (MQ=255)                                                                                                  cGTCTGGCTGGGGGTGGAGTGGGAGACGCGGTTCTACCGCTACAGCTACTCGTTTGATGGGGAGACGTGGCACACCGTGCCGGGGACGTAGGGGGCGTGGAGGCTGGCGGACGACTCCATCGGCGGGCGC‑GGCT‑TC‑TT‑CC‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAgg                                                                                                                                          <  2:861050/238‑1 (MQ=255)                                                                                                  cGTCTGGCTGCGGGGGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAgg                                                                                                                                          <  1:444030/238‑1 (MQ=255)                                                                                                  cGTCGGGCTGGGGGGGGACGGGGAGCCGCTGGTCGCCCGCTACAGCTGCTCGTTTGAGGGCGGGACGTGGCCCCCCGGGCCGGTGACGTATGGGGCGGGGAAGCTGTCGGACGACTACAGCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTCCGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAgg                                                                                                                                          <  2:598911/238‑1 (MQ=255)                                                                                                  cGTCGGGCTGCGGGTGGAGGTGGATACGTTGGTCTACCGCTACAGCTACTCGTTTGAGGGCGAGACGTGGCACACGGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAgg                                                                                                                                          <  1:1482137/238‑1 (MQ=255)                                                                                                  cGTCGGGCGGCGGGTGGGCGTGGATACGGTGGTCTACCGCTACGGCTCCTCGTTTGAGGGCGAGCCGTGCCCCCCCGTGCCGTGGACGTATGAGGCGTGGAAGCTGTCGGACGACTACTGCGGCGGGCGC‑GTCT‑TC‑TT‑CC‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTCCGCGGACTTCGCCTACTTCCCCTACGAGCCGGTCTAACGGCTCAgg                                                                                                                                          <  2:49718/238‑1 (MQ=255)                                                                                                      tGGCGGGGGGGGGACGGGTGTACGCCGGTCTGCCGCTACGGCTACTCGTTGGAGGGCGGGACGGGGCCCCCCGCGCCGGTGCCGTATGGGGCGTGGAAGCTGGCGGACGACTACATCGGCGGGCGG‑GGCT‑TC‑TC‑CA‑CCGGCGCGTTTGGGGGCCGGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGGGGACTTCGACTCCTCCACCTCCGAGCCGGTCTAACGGCTCAGGCCgg                                                                                                                                      <  2:824030/238‑1 (MQ=255)                                                                                                                       ggttACGCGGGTTTCCCGCTCCCGCTACTCGTTTGAGGGGGAGCCGTGGCGCCCCGGGCGGGTGCCGTAGGGGGCGTGGAAGCTGTGGGACGACTACATCGGGGGGCGC‑GGCT‑TC‑TT‑CC‑CCGGCGCGTTTGTGGGCCTGCACGGCGAGGACATCAGCGGCGACGGCTGTTACGCGGACTCCGACTCCTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCga                                                                                                                    <  2:1240140/236‑1 (MQ=255)                                                                                                                            ggCTGGTCGCCCGGTACGGCTACTGGTTTGAGGGCGGGACGGGGCCCCCCGTGCGGGTGACGTATGAGGCGGGGAGGCTGGCGGCCGATTACATCGGGGGGCGG‑GGCT‑TC‑TT‑CG‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTCCGACTCCTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAgcg                                                                                                               <  2:1093243/238‑1 (MQ=255)                                                                                                                                  tCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGc                                                                                                          <  1:1267961/238‑1 (MQ=255)                                                                                                                                   cTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGcc                                                                                                         <  1:745092/238‑1 (MQ=255)                                                                                                                                     cccGCTACGGCGACTCGTTTGATGGCGAGACGTGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACACCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCCCTGCGAGGACACCAGCGGCGACGGCTGCTCCGCGGACTTCGACTCCTCCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCgg                                                                                                      <  2:783959/238‑1 (MQ=255)                                                                                                                                     aCCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTGGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTCCTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCgg                                                                                                      <  1:1264040/239‑1 (MQ=255)                                                                                                                                     aCCGCTACAGCTACGTGTTGGAGGGCGAGGCGGGGCACACCGTGCCGGGGAGGTATGAGGCGGTGAAGCTGTCGGCCGACTACATCGGCGGGGGC‑GGCT‑TC‑TC‑CG‑CCGGCGCGTTTGTGGGCGTGCACTGCGAGGACATCGGCGGCGACGGCTGCTACGCGGACTTCGACTCCTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCgg                                                                                                      <  2:39733/239‑1 (MQ=255)                                                                                                                                                    cGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGGGC‑GGCG‑TC‑TT‑CT‑CCCGCGCGGTTGTGGGGCTCCCCCGCGCGGGGGTCCAGGGCGCCGCGCGCCCCGCCGGCCTCCGCCCCTTCACCTACCAGCCGGGCCAACGGCTCCGGCCGGGGCGCCCCCCCCGCGCGACCGCGGGAGGCCGCCCCGCTGGAgctgc                                                                                        >  2:857464/1‑238 (MQ=255)                                                                                                                                                                  gcGTGGCCCGCCGTGCCGGTGACTTATGAGGCGTGGAAGCTGTCGGACGGCTCGATCGGCGGGCGG‑GGCT‑TC‑TT‑CC‑CCGGGGGGTTTGGGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCGCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGc                                                                          <  2:1437571/237‑1 (MQ=255)                                                                                                                                                                  cgcgggccccgcgggcgggTGACGTGGGGGGGGTGGAGGCTGTCGGACGCCTCCACCGGCGGGCGC‑GGCT‑TC‑TT‑CC‑CCGGCGCGTTTGTGGGCTTGCACTGCGAGGACATCGGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGc                                                                          <  2:1497976/221‑1 (MQ=255)                                                                                                                                                                  aCGTGGCACGCGGTGCCGGTGCCGAAGGGGGGGTGAAGGTTGTGGGACGCCTCGAACGGGGGGCGC‑GTTT‑TT‑TT‑CC‑CCGGGGCGTTTGTGGGCCTGCACTGCGAGAACACCAGCGGCGGCGGCTGCTCCGCGGACTTCGACTCCTCCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGc                                                                          <  2:630546/238‑1 (MQ=255)                                                                                                                                                                     tGGCACCCCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAg                                                                                                                                                            <  1:721983/153‑1 (MQ=255)                                                                                                                                                                             ggTGCGGGGGCCGTAGGGGGGGGGGAACCGGGCGGACGACTCCATCGGCGGGCGG‑GCCT‑TC‑TT‑CC‑CGGGCGCGTTTGTGGGCCTGCACTGGGAGGACATCAGCGGCGCCGGCTGCTACGCGGACTTCGACTCCTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGcgtc                                                              <  2:468368/238‑1 (MQ=255)                                                                                                                                                                                                           gTCGGACGACTACAGCGGGGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCt                                <  1:815828/239‑1 (MQ=255)                                                                                                                                                                                                            gtcggagacctcggtggggggggC‑GGGC‑TC‑TTCCA‑CGGGCGGGTTGGTGGGCCGGCCCTGCGAGGACACGAGCGGGGCCGGCGGCTACGGGGACTGCGACTACTTCACCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGc                                 <  2:1429627/220‑1 (MQ=255)                                                                                                                                                                                                                              gggcgg‑ggcg‑TC‑TT‑CA‑CGGGGGCGTTGGGGGGCCTGCACTGGGAGGACAGCAGCGGCGACGGCGGCGACGCGGGCTGCGACTACTTCCCCGACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGAcgcg              <  2:68245/238‑1 (MQ=255)                                                                                                                                                                                                                               ggcgC‑GGCT‑TC‑TT‑CC‑CCGGCGCGTTTGGGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCgg             <  1:637873/238‑1 (MQ=255)                                                                                                                                                                                                                                gcgc‑GGCT‑TC‑TT‑CA‑CCGGCGGGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGCCTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCggg            <  1:444205/238‑1 (MQ=255)                                                                                                                                                                                                                                      gCT‑TC‑TT‑CG‑CCGGCGCGTTTGGGGGCCTGCCCGGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCg       <  1:545590/238‑1 (MQ=255)                                                                                                                                                                                                                                      ctt‑tc‑tt‑cC‑CCGGCGCGTTTGGGGGCCTGCACTGCGAGGACATCAGCGGCGCGGGCGTATACCCGGCCTCCGACTACTCCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCg       <  1:161417/236‑1 (MQ=255)                                                                                                                                                                                                                                             tt‑cA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGtt  >  1:79980/1‑238 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            GACAGGGTAGAACCATCAAAGTTATCCAGCTCGACCACAACGTGCCGTCGTGGCCGCTGCACGAGCAGCCCATTCCGGTGCCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGCGC‑GGCT‑TC‑TT‑CA‑CCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTT  >  NC_000913/288074‑288544

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