Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,264,912	T→A	6.3%	L42Q (CTG→CAG)	vanB →	vanillate O‑demethylase oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,264,912	0	T	A	6.3%	194.7 / 2.8	63	L42Q (CTG→CAG)	vanB	vanillate O‑demethylase oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base A (2/2);  total (33/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 = 9.48e-01

CCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTAT  >  NC_002947/4264767‑4265059                                                                                                                                                   |                                                                                                                                                    cctgcctgAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGcc                                                                                                                                                  >  1:306545/1‑149 (MQ=255)      ctgAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACg                                                                                                                                             >  6:123480/1‑149 (MQ=255)       tgAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGCGCACCAGCCCGTCgg                                                                                                                                            >  4:80231/1‑149 (MQ=255)        gAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGc                                                                                                                                                   <  4:220098/141‑1 (MQ=255)        gAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGc                                                                                                                                                   >  3:220098/1‑141 (MQ=255)         aGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGc                                                                                                                                          <  1:30306/149‑1 (MQ=255)           gAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTg                                                                                                                                        >  6:80464/1‑149 (MQ=255)             aCAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGCGCACCAGCCCGTCGGGCAGGt                                                                                                                                      <  3:80231/149‑1 (MQ=255)               aGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTgc                                                                                                                                    >  8:316513/1‑149 (MQ=255)                 ccGACATGATGGATGGCGTAGTGGTATCCCGTATGGATGAAGCGCAGGGTATCTGCAGCTTCGAGGTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTgcgc                                                                                                                                  <  2:34161/149‑1 (MQ=255)                 ccGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTgcgc                                                                                                                                  <  7:316513/149‑1 (MQ=255)                 ccGACATGATCGATGCCGTAGTGGCATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCAGCCCGTCGGGCAGGTgcac                                                                                                                                  <  3:129653/149‑3 (MQ=255)                          tCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCg                                                                                                                         >  1:251661/1‑149 (MQ=255)                            gATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGc                                                                                                                        <  2:289500/148‑1 (MQ=255)                             aTGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCtg                                                                                                                      <  3:150344/149‑1 (MQ=255)                             aTGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCtg                                                                                                                      >  2:251573/1‑149 (MQ=255)                                 cGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGACCGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGc                                                                                                                   <  1:251573/149‑1 (MQ=255)                                            cccGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCga                                                                                                        <  7:192631/148‑1 (MQ=255)                                                tAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAAc                                                                                                   >  1:253683/1‑149 (MQ=255)                                                      tGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATc                                                                                             <  5:92210/149‑1 (MQ=255)                                                          gcgcAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCCCCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCt                                                                                          >  8:286820/1‑148 (MQ=255)                                                            gcAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATc                                                                                       <  5:123480/149‑1 (MQ=255)                                                             cAGGGTATCTGCAGCTTCGAGCTGGCCGCGGTAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCt                                                                                      >  3:99438/1‑149 (MQ=255)                                                                    ccgtgcgCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGc                                                                               <  4:247381/143‑1 (MQ=255)                                                                        cAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCg                                                                              <  5:80464/146‑1 (MQ=255)                                                                        cAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTAc                                                                           <  5:30056/149‑1 (MQ=255)                                                                            ttCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCaa                                                                       >  2:86624/1‑149 (MQ=255)                                                                              cGAGCTGGCCGCGGTAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCaa                                                                       <  4:99438/147‑1 (MQ=255)                                                                                  cTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGAccc                                                                 <  5:322207/149‑1 (MQ=255)                                                                                    ggCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCgg                                                               >  8:64369/1‑149 (MQ=255)                                                                                          ggCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTcgc                                                         >  6:135681/1‑149 (MQ=255)                                                                                             aGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCggg                                                      >  4:347553/1‑149 (MQ=255)                                                                                                 ggCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCgg                                                  >  3:77778/1‑149 (MQ=255)                                                                                                   cAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGtt                                                >  7:101567/1‑149 (MQ=255)                                                                                                   cAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGt                                                 >  6:29999/1‑148 (MQ=255)                                                                                                    aGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTtc                                               <  2:270740/149‑1 (MQ=255)                                                                                                       ctgctgCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCg                                            <  7:286820/149‑1 (MQ=255)                                                                                                        tgctgcCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGt                                           >  2:183793/1‑149 (MQ=255)                                                                                                          gcgcCGGCGCTGAACGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGAccc                                                                 <  6:321581/123‑1 (MQ=255)                                                                                                          gcgcCGGCGCTGAACGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGAccc                                                                 >  5:321581/3‑125 (MQ=255)                                                                                                          ctgcCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAg                                         <  7:64369/149‑1 (MQ=255)                                                                                                           tgcCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGtt                                                >  6:167536/1‑141 (MQ=255)                                                                                                           tgcCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGtt                                                <  5:167536/141‑1 (MQ=255)                                                                                                           tgcCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGc                                        >  5:280376/1‑149 (MQ=255)                                                                                                             ccGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGc                                        <  2:253683/147‑1 (MQ=255)                                                                                                             ccGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGc                                        <  2:306545/147‑1 (MQ=255)                                                                                                              cGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTg                                     >  6:6499/1‑149 (MQ=255)                                                                                                               ggCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTg                                     >  4:212166/1‑148 (MQ=255)                                                                                                                 cGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCAc                                  <  5:29999/149‑1 (MQ=255)                                                                                                                   ttCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACg                                 >  2:338271/1‑148 (MQ=255)                                                                                                                      aGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAAc                              <  5:135681/148‑1 (MQ=255)                                                                                                                                 ccATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCg                  >  1:330992/1‑149 (MQ=255)                                                                                                                                 ccATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCg                  >  8:63086/1‑149 (MQ=255)                                                                                                                                 ccATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCg                  >  5:150200/1‑149 (MQ=255)                                                                                                                                   atatCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGt                >  6:97049/1‑149 (MQ=255)                                                                                                                                    tatCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTg               <  5:6499/149‑1 (MQ=255)                                                                                                                                        gA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCg           <  7:63086/149‑1 (MQ=255)                                                                                                                                         a‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCgg          >  5:341047/1‑149 (MQ=255)                                                                                                                                          tggcgCACCAGCCCGTCGGGCAGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCgg          >  8:240431/5‑149 (MQ=255)                                                                                                                                               cACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCgg          >  5:146858/1‑144 (MQ=255)                                                                                                                                               cACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCgg          <  6:146858/144‑1 (MQ=255)                                                                                                                                                aCCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCg     >  7:72729/1‑148 (MQ=255)                                                                                                                                                  cTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTAt  <  8:72729/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCTGCCTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTAT  >  NC_002947/4264767‑4265059

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