Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,053,630	A→T	7.7%	intergenic (‑30/+36)	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,630	0	A	T	7.7%	190.5 / 5.5	65	intergenic (‑30/+36)	PP_1828/PP_1829	hypothetical protein/alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base A (32/28);  new base T (3/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

GGCAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTGG  >  NC_002947/2053485‑2053775                                                                                                                                                  |                                                                                                                                                   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                                                                                                                 <  4:241133/148‑1 (MQ=255)                                 ttCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGa                                                                                                                 <  7:351193/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                                                                                                                                               <  6:43978/87‑1 (MQ=255)                                                                 gggTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCggc                                                                                                                                               >  5:43978/1‑87 (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                                 >  8:94063/1‑121 (MQ=255)                                                                                                                                            tgctgTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAgg                                 <  7:94063/121‑1 (MQ=255)                                                                                                                                            tgctgTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCt     >  4:358220/1‑149 (MQ=255)                                                                                                                                              ctgTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCtt    <  7:46762/148‑1 (MQ=255)                                                                                                                                                 tATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTgg  <  1:283426/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GGCAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCC‑TGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTGG  >  NC_002947/2053485‑2053775

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