Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,053,627	T→C	8.3%	intergenic (‑27/+39)	PP_1828 ← / ← PP_1829	hypothetical protein/alpha/beta family hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,053,627	0	T	C	8.3%	170.9 / 5.5	60	intergenic (‑27/+39)	PP_1828/PP_1829	hypothetical protein/alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base C (3/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

TCGGTCGGCAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTT  >  NC_002947/2053479‑2053773                                                                                                                                                     |                                                                                                                                                    tcggtctgCAGTGCTGACAGGGTTCGCAGCTTGAAATGTTCGGCGGTCAGCTAGGCCCAGGGCAGTGCGGGGGTTTCGTGGTTTTCTTCGGCCGGCTCGAGTGCGGTTGTTGAGGTGGTCATCAGTGCTTGTCCTCAGCTATGCctgct                                                                                                                                                     <  8:91917/149‑1 (MQ=255)       ggCAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCg                                                                                                                                               >  8:351193/1‑149 (MQ=255)         cAGTGCTGCCAGGCGTAGCAGCTGGAGATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCgg                                                                                                                                              >  7:213137/1‑148 (MQ=255)             gctgcCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCggcagg                                                                                                                                          >  7:62623/1‑148 (MQ=255)             gctgcCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCggcagg                                                                                                                                          >  3:241133/1‑148 (MQ=255)                gcCAGGCGTCGCAGCTGGAGATGTGCGGCGGTCCGCGCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGCCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATg                                                                                                                                      <  8:213137/149‑1 (MQ=255)                 ccAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATgg                                                                                                                                     >  7:272800/1‑149 (MQ=255)                  ctggCGTCGCAGCGGGAAAGGTGCGCCGGCCCCCGCGCCCCAGGGCCGGGGGGGGGTTTCCGGGGGGTCGGCGGCCGGCGCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGc                                                                                                                                    <  3:350957/147‑1 (MQ=255)                   aGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCt                                                                                                                                   >  5:375916/1‑149 (MQ=255)                     gCGTAGCTGCTGGAAATGTGCGGCGGTCAGCTCGGCCCAGGGCTGGGCGGGGGTTTCCTGGGGGTCTTCGGCCGGCGCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTga                                                                                                                                 <  2:140905/149‑1 (MQ=255)                     gCGTAGCAGCTGGAAATGTGCGGCGGTCAGCGCGGCCCCGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCCGATACAGCAGGCATAGCTga                                                                                                                                 <  3:111285/149‑1 (MQ=255)                            aggtGGAAATGTTGGGGGTGCAGCTGGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGCCGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGa                                                                                                                          <  6:299966/146‑1 (MQ=255)                             gcTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAg                                                                                                                         >  1:19669/1‑149 (MQ=255)                                   aaggtccgggggTCTCCTCGGCCCTGGGACGGGCGGGGGTTGGCGGGGGTTCTGCGGCCGGCTGGCCGGGGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGt                                                                                                                   <  8:220065/142‑1 (MQ=255)                                       ttCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGa                                                                                                               <  7:351193/148‑1 (MQ=255)                                       ttCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGa                                                                                                               <  4:241133/148‑1 (MQ=255)                                          ggcggGCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTa                                                                                                           >  7:331980/1‑149 (MQ=255)                                              tctaGCGCGGCCCCGGGCCGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACg                                                                                                         <  4:10187/144‑1 (MQ=255)                                                  gCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACg                                                                                                         >  3:10187/1‑143 (MQ=255)                                                  gCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAg                                                                                                   <  8:331980/149‑1 (MQ=255)                                                        cccAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAAt                                                                                             >  4:141702/1‑149 (MQ=255)                                                           agggcagggcGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCt                                                                                          >  8:320053/1‑149 (MQ=255)                                                             ggcagggcGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTg                                                                                         >  2:161850/1‑148 (MQ=255)                                                                agggcGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGAt                                                                                      >  3:320875/1‑148 (MQ=255)                                                                   gcGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAAc                                                                                  >  8:286443/1‑149 (MQ=255)                                                                   gcGGGGGTTTCCTGGGGTGCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGTCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAAc                                                                                  <  3:44423/149‑1 (MQ=255)                                                                       gggTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCggc                                                                                                                                             >  5:43978/1‑87 (MQ=255)                                                                       gggTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCggc                                                                                                                                             <  6:43978/87‑1 (MQ=255)                                                                       ctgTTTCCTTGGGTTCTTCGGCCGGCTCTACGGCGGTTGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACGTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTc                                                                               <  4:120453/146‑1 (MQ=255)                                                                        ggTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTcac                                                                             <  8:272800/149‑1 (MQ=255)                                                                         gTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTcac                                                                             >  4:320071/1‑148 (MQ=255)                                                                          tttCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTc                                                                               >  3:120453/1‑145 (MQ=255)                                                                                  ggTTCTTCGGCCGGCTCGACGGCGGGGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATtcgctcg                                                                   <  6:375916/149‑1 (MQ=255)                                                                                      cttcGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCa                                                               >  2:283426/1‑149 (MQ=255)                                                                                         cggccggcTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAAt                                                                          <  5:332265/135‑1 (MQ=255)                                                                                         cggccggcTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAAt                                                                          >  6:332265/1‑135 (MQ=255)                                                                                           gccggcTCGACGGAGGTGGGGGAGGTCGTCAACAGGCAATCTACTCACCTATGCCTGCTGTATCGGCACGCATGGCTGAGAACTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGCATGAACCTCACAATACGCTCGCCCAGCtaa                                                          >  7:14748/1‑147 (MQ=255)                                                                                            ccggcTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCAc                                                         <  2:19669/149‑1 (MQ=255)                                                                                            ccggcTCGACGGCGGGGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTTAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCAc                                                         <  8:62623/149‑1 (MQ=255)                                                                                                cTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACgcgc                                                     <  3:141702/149‑1 (MQ=255)                                                                                                            gggggAGGTGGTCATCAGGGCCTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTc                                         >  8:42282/1‑149 (MQ=255)                                                                                                               ggAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCt                                          >  5:361079/1‑145 (MQ=255)                                                                                                               ggAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCt                                          <  6:361079/145‑1 (MQ=255)                                                                                                               ggAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCggg                                      <  1:52902/149‑1 (MQ=255)                                                                                                                  ggtggtCATCAGGGCGTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGtt                                   <  8:312438/149‑1 (MQ=255)                                                                                                                   gtggtCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAg                                            >  7:296261/1‑139 (MQ=255)                                                                                                                   gtggtCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAg                                            <  8:296261/139‑1 (MQ=255)                                                                                                                    tggtCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGa                                 <  7:320053/149‑1 (MQ=255)                                                                                                                      gtCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAgg                               >  5:248253/1‑149 (MQ=255)                                                                                                                       tcatcaGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAgg                               >  1:52685/1‑148 (MQ=255)                                                                                                                              ggCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTc                       <  3:320071/149‑1 (MQ=255)                                                                                                                                    tCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctg                 >  4:244646/1‑149 (MQ=255)                                                                                                                                     tcTCAGCCATGCCTGCCGATACAGCAGGCATAGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGCGTTGAGGTAGCTGTCGCgctgc                >  4:40849/2‑149 (MQ=38)                                                                                                                                      cTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctgct               >  8:46762/1‑149 (MQ=255)                                                                                                                                       tCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCGTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctgct               >  1:76238/1‑149 (MQ=255)                                                                                                                                       tCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGctgctg              >  5:191903/1‑149 (MQ=255)                                                                                                                                       tCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGctgctg              >  2:355189/1‑149 (MQ=255)                                                                                                                                       tCAGCCATGCCTGCCGATACAGCAGGCATAGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTcgc                     <  4:91988/142‑1 (MQ=38)                                                                                                                                       tCAGCCATGCCTGCCGATACAGCAGGCATAGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTcgc                     >  3:91988/1‑142 (MQ=38)                                                                                                                                         aGCCATGCCTGCCGATACAGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGc            >  7:34000/1‑149 (MQ=255)                                                                                                                                           cTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTc          <  7:286443/149‑1 (MQ=255)                                                                                                                                                  tgctgTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAgg                               <  7:94063/121‑1 (MQ=255)                                                                                                                                                  tgctgTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAgg                               >  8:94063/1‑121 (MQ=255)                                                                                                                                                  tgctgTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCt   >  4:358220/1‑149 (MQ=255)                                                                                                                                                    ctgTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCtt  <  7:46762/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TCGGTCGGCAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTT  >  NC_002947/2053479‑2053773

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: ‑