Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	508,771	Δ1 bp	6.9%	coding (240/1890 nt)	estP ←	esterase EstP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	508,771	0	C	.	6.9%	226.3 / 10.3	72	coding (240/1890 nt)	estP	esterase EstP
Reads supporting (aligned to +/- strand):  ref base C (36/30);  new base . (2/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

CCGGCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACAC  >  NC_002947/508625‑508914                                                                                                                                                   |                                                                                                                                                 ccGGCATTGGGGTTTCCTGGAGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGgcgc                                                                                                                                                >  6:336286/1‑149 (MQ=255)  cGGCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGgcgc                                                                                                                                                >  3:353618/1‑148 (MQ=255)  cGGCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGcgcg                                                                                                                                               >  7:12/1‑149 (MQ=255)    gCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGAc                                                                                                                                             <  6:284294/149‑1 (MQ=255)      aTTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCTACGGGTTAAGGTCGCCTGGCGCGACGc                                                                                                                                           <  7:146979/149‑1 (MQ=255)            gTTTCCTGGGGGGATGACGGTCTCGGACGTGGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATGGGGCTGGATACCTACGGATGTCGACGGGTTACGGTCGCCTGGCGCGCGGC‑CCAggc                                                                                                                                    >  4:198927/1‑149 (MQ=255)            gTTTCCTGGAGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAggc                                                                                                                                    <  5:336286/149‑1 (MQ=255)             tttCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGG‑GCGTCGCGCCAggcg                                                                                                                                   <  6:347880/149‑1 (MQ=255)               tCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGCCGCCTGG‑GCGTCGCGCCAggcga                                                                                                                                  <  4:278517/148‑2 (MQ=255)                 cTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCa                                                                                                                               >  7:201030/1‑149 (MQ=255)                      gggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAg                                                                                                                          >  5:377340/1‑149 (MQ=255)                       ggaggaCGGTCTCGGACGTTGTCGTGATCGACTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAgg                                                                                                                         <  3:198927/149‑1 (MQ=255)                       ggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAgg                                                                                                                         <  7:182309/149‑1 (MQ=255)                       ggATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCGGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACGTCCGGATGTCGCCGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAgg                                                                                                                         <  5:308253/149‑1 (MQ=255)                                    ggACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTAGTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTACGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGaca                                                                                                            >  8:323891/1‑149 (MQ=255)                                    ggACGTTGTCGTGATCGAGTACAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGaca                                                                                                            >  2:269902/1‑149 (MQ=255)                                         ttGTCGTGCTCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGGCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCggg                                                                                                       <  2:343546/149‑1 (MQ=255)                                            tCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGc                                                                                                     >  7:263786/1‑148 (MQ=255)                                              gTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCg                                                                                                    >  3:113582/1‑147 (MQ=255)                                                 aTCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGtt                                                                                               >  2:108876/1‑149 (MQ=255)                                                          aGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGc                                                                                      >  5:184874/1‑149 (MQ=255)                                                          aGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGc                                                                                      >  7:26751/1‑149 (MQ=255)                                                          aGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGG‑GCGTCGCGCCAggcgac                                                                                                                                 >  3:123507/1‑104 (MQ=255)                                                          aGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGG‑GCGTCGCGCCAggcgac                                                                                                                                 <  4:123507/106‑3 (MQ=255)                                                                tgctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACg                                                                                >  7:316527/1‑149 (MQ=255)                                                                  ctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCa                                                                                  <  6:377340/145‑1 (MQ=255)                                                                  ctgGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACg                                                                                <  1:365858/147‑1 (MQ=255)                                                                  ctgGGGGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTTGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGt                                                                              <  2:331167/149‑1 (MQ=255)                                                                      gtggtgTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGt                                                                          <  7:340376/149‑1 (MQ=255)                                                                       tggtgTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTg                                                                         <  8:316527/149‑1 (MQ=255)                                                                           gtACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAc                                                                     >  6:77767/1‑149 (MQ=255)                                                                            tACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAc                                                                     <  1:269902/148‑1 (MQ=255)                                                                             aCCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgc                                                                   >  2:211433/1‑149 (MQ=255)                                                                               ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                 >  5:250456/1‑149 (MQ=255)                                                                               ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                 >  5:249296/1‑149 (MQ=255)                                                                               ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                 <  8:37611/149‑1 (MQ=255)                                                                               ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                 >  5:144497/1‑149 (MQ=255)                                                                               ccgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgc                                                                 >  2:122826/1‑149 (MQ=255)                                                                                cgccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAAcgcgcg                                                                <  8:201030/149‑1 (MQ=255)                                                                                 gccgACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACgcgcgc                                                               >  4:88998/1‑149 (MQ=255)                                                                                      cTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTggg                                                          <  8:12/149‑1 (MQ=255)                                                                                       tGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTgggg                                                         >  6:306954/1‑149 (MQ=255)                                                                                          ccAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTAc                                                      >  6:128084/1‑149 (MQ=255)                                                                                             gTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGc                                                   <  4:333515/149‑1 (MQ=255)                                                                                               tattaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGgcg                                                                                                                                                 >  7:158489/1‑54 (MQ=255)                                                                                               tattaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGgcg                                                                                                                                                 <  8:158489/54‑1 (MQ=255)                                                                                                 ttaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaa                                                <  1:108876/148‑1 (MQ=255)                                                                                                 ttaCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaa                                                <  3:88998/148‑1 (MQ=255)                                                                                                   aCCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaaca                                              <  6:144497/148‑1 (MQ=255)                                                                                                    ccATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGt                                                                                          <  3:336975/103‑1 (MQ=255)                                                                                                    ccATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGt                                                                                          >  4:336975/1‑103 (MQ=255)                                                                                                    ccATAGGGCTGGATACATACGGATGTCGACGGGTTAAGGTCGACTGTCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaacaa                                             <  6:272033/148‑1 (MQ=255)                                                                                                     cATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCcaac                                               >  5:234902/1‑145 (MQ=255)                                                                                                     cATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGa                                           <  6:234902/149‑1 (MQ=255)                                                                                                       tCGGGCTGTATACCTACGGATGTCGACGGGTTAAGTTCGCCTGGAGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGTTACCGCCAACAAGATc                                         <  8:263786/149‑1 (MQ=255)                                                                                                         gggCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCgg                                       <  5:77767/149‑1 (MQ=255)                                                                                                          ggCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCgg                                       >  3:250170/1‑148 (MQ=255)                                                                                                             tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                   >  5:357331/1‑149 (MQ=255)                                                                                                             tGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                   >  1:135665/1‑149 (MQ=255)                                                                                                               gATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGaaa                                   <  4:353618/147‑1 (MQ=255)                                                                                                               gATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACt                                 >  4:204832/1‑149 (MQ=255)                                                                                                                    cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                         <  5:207014/136‑1 (MQ=255)                                                                                                                    cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATc                                         >  6:207014/1‑136 (MQ=255)                                                                                                                    cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCg                            >  3:115196/1‑149 (MQ=255)                                                                                                                    cTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCg                            >  2:39922/1‑149 (MQ=255)                                                                                                                    cTACGGATGTCGACGGGATAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCg                            >  5:27305/1‑149 (MQ=255)                                                                                                                     tACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATAGGGAAACTGCCCgg                           >  2:267485/1‑149 (MQ=255)                                                                                                                               gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGa                  <  6:250456/148‑1 (MQ=255)                                                                                                                               gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGa                  <  6:184874/148‑1 (MQ=255)                                                                                                                               gACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGa                  <  5:306954/148‑1 (MQ=255)                                                                                                                                    gTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGt            <  1:211433/149‑1 (MQ=255)                                                                                                                                     ttAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTc           >  1:151622/1‑149 (MQ=255)                                                                                                                                        aGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCCTCGGCGAGGCTGTCGc         >  4:309625/1‑148 (MQ=255)                                                                                                                                         ggccGCCTGG‑GCGTCGCGCCAGGCGACCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGcc        >  2:323080/4‑148 (MQ=255)                                                                                                                                          gTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCa       >  6:75963/1‑148 (MQ=255)                                                                                                                                               cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAacac  >  7:260512/1‑148 (MQ=255)                                                                                                                                               cTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAacac  >  1:255703/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CCGGCATTGGGGTTTCCTGGGGGGATGACGGTCTCGGACGTTGTCGTGATCGAGTCCAGGATCTGCTGGGTGGTGTACCCGCCGACTGCCCAGTTATTACCATCGGGCTGGATACCTACGGATGTCGACGGGTTAAGGTCGCCTGGCGCGACGC‑CCAGGCGGCCACCGAGGATCATCGGCGACACCGGGGCGAAGTTGCCGTCGGCGTCACGGTTGGTGAAACGCGCGCCTGGGGTACCGCCAACAAGATCGGGAAACTGCCCGGCATCGGCGAGGCTGTCGCCAAACAC  >  NC_002947/508625‑508914

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