Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,264,918	A→T	6.5%	D44V (GAC→GTC)	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,918	0	A	T	6.5%	177.9 / 2.2	61	D44V (GAC→GTC)	vanB	vanillate O‑demethylase oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base T (2/2);  total (32/29)
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.21e-01

CTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCAGTG  >  NC_002947/4264772‑4265064                                                                                                                                                    |                                                                                                                                                   ctgAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACg                                                                                                                                                  >  6:123480/1‑149 (MQ=255)  tgAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGCGCACCAGCCCGTCgg                                                                                                                                                 >  4:80231/1‑149 (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                                                                                                                           >  2:251573/1‑149 (MQ=255)                        aTGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCtg                                                                                                                           <  3:150344/149‑1 (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                                                     <  5:167536/141‑1 (MQ=255)                                                                                                      tgcCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGtt                                                     >  6:167536/1‑141 (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)                                                                                                                                                  cGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCAGTg  >  2:29113/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CTGAGGAACAGCCGACATGATCGATGCCGTAGTGGTATCCCGTAACGATGAAGCGCAGGGTATCTGCAGCTTCGAGCTGGCCGCGGCAGATGGCAGCCTGCTGCCGGCGTTCAGCGCCGGCGCCCATATCGA‑CGTGCACCTGCCCGACGGGCTGGTGCGCCAGTATTCGCTGTGCAACCACCCCGAAGAACGCCATCGCTATCTGATTGGCGTACTCAACGACCCGGCTTCGCGGGGCGGTTCTCGTAGCCTGCACGAACAGGTGCAAGCCGGTGCCCGGCTGCGTATCAGTG  >  NC_002947/4264772‑4265064

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