Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,866,221	C→T	8.7%	V186I (GTC→ATC)	glllA ←	gallate dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,866,221	0	C	T	8.7%	161.4 / 4.0	57	V186I (GTC→ATC)	glllA	gallate dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base T (3/2);  total (30/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.31e-01

CGCCGAGCGTGGCGTACTCGGCCAGTGTCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGCTTTCATCGCAG  >  NC_002947/2866074‑2866369                                                                                                                                                     |                                                                                                                                                     cgccgAGCGTGGCGTACTCGGCCAGTGTCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTGTTTTAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGcgac                                                                                                                                                      <  1:107624/149‑1 (MQ=255) cgccgAGCGTGGCGTACTCGGCCAGTGTCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGcgac                                                                                                                                                      >  8:125240/1‑148 (MQ=255)        cGTGGCGTACTCGGCCAGTGTCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCa                                                                                                                                              >  3:42214/1‑149 (MQ=255)               tACTCGGCCAGTGTCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAg                                                                                                                                       <  4:86415/149‑1 (MQ=255)                   cGGCCAGTGTCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGCTTG‑TTGAAGCCGTAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCt                                                                                                                                   >  4:265860/1‑149 (MQ=255)                       cAGTGTCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGACAAGGAACTGCGCGTCCCACTCTGGGATG‑TTGAAGCCGCAACGCTCCCCGTGCACCAGCTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGttgatct                                                                                                                               >  6:104616/1‑143 (MQ=255)                        aGTGTCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCggg                                                                                                                              <  5:400992/149‑1 (MQ=255)                        aGTGTCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCggg                                                                                                                              >  8:297689/1‑149 (MQ=255)                          tgtCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCgggg                                                                                                                             >  5:109524/1‑148 (MQ=255)                            tCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAg                                                                                                                          >  2:26358/1‑149 (MQ=255)                              aTCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGct                                                                                                                        >  2:226680/1‑149 (MQ=255)                                     tAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATg                                                                                                                 <  6:177525/149‑1 (MQ=255)                                     tAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATg                                                                                                                 <  2:431408/149‑1 (MQ=255)                                       aGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATggc                                                                                                               <  4:403853/149‑1 (MQ=255)                                       aGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATgg                                                                                                                <  7:278020/148‑1 (MQ=255)                                       aGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACGCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATgg                                                                                                                <  8:393119/148‑1 (MQ=255)                                         gcgcTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATggcgc                                                                                                             >  7:26624/1‑149 (MQ=255)                                                 ggTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCTGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAAt                                                                                                     <  2:54639/149‑1 (MQ=255)                                                     gTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCt                                                                                                 <  5:359616/149‑1 (MQ=255)                                                      ttCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTg                                                                                                >  6:63983/1‑149 (MQ=255)                                                        cACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCaa                                                                                                      <  2:15441/141‑1 (MQ=255)                                                        cACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCaa                                                                                                      >  1:15441/1‑141 (MQ=255)                                                        cACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTgg                                                                                               <  4:42214/148‑1 (MQ=255)                                                        cACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGc                                                                                              >  7:4021/1‑149 (MQ=255)                                                         aCCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGATCGTCGCCAccg                                                                                                                                           <  2:330293/103‑2 (MQ=255)                                                         aCCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGATCGTCGCCAccg                                                                                                                                           >  1:330293/1‑102 (MQ=255)                                                         aCCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGcc                                                                                             >  4:316790/1‑149 (MQ=255)                                                         aCCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGcc                                                                                             >  5:231892/1‑149 (MQ=255)                                                               agGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGttt                                                                                       <  6:109524/149‑1 (MQ=255)                                                                  tCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTa                                                                                    <  6:153207/149‑1 (MQ=255)                                                                   cAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTag                                                                                   >  8:198682/1‑149 (MQ=255)                                                                     aGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTagca                                                                                 <  2:371573/149‑1 (MQ=255)                                                                             cgcgTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAgcgccgc                                                                          <  3:412321/148‑1 (MQ=255)                                                                                 tCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGt                                                                     >  1:9652/1‑149 (MQ=255)                                                                                   ccACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTgg                                                                   <  8:4021/149‑1 (MQ=255)                                                                                   ccACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTgg                                                                   <  2:9652/149‑1 (MQ=255)                                                                                    cACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTg                                                                                                >  3:78805/1‑119 (MQ=255)                                                                                    cACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTg                                                                                                <  4:78805/119‑1 (MQ=255)                                                                                       tctGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATg                                                                                                    <  4:356518/112‑1 (MQ=255)                                                                                       tctGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATg                                                                                                    >  3:356518/1‑112 (MQ=255)                                                                                             ttg‑ttgAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAg                                                         >  3:202538/1‑149 (MQ=255)                                                                                             ttg‑ttgAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAg                                                         <  7:297689/149‑1 (MQ=255)                                                                                                   gAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTg                                                      >  6:249105/1‑147 (MQ=255)                                                                                                   gAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGc                                                     <  7:125240/148‑1 (MQ=255)                                                                                                    aaGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTggg                                                                  <  7:96979/134‑1 (MQ=255)                                                                                                    aaGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTggg                                                                  >  8:96979/1‑134 (MQ=255)                                                                                                         gCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGc                                              >  2:177773/1‑149 (MQ=255)                                                                                                         gCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGc                                              >  2:356256/1‑149 (MQ=255)                                                                                                              gCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGAcc                                         >  2:412576/1‑149 (MQ=255)                                                                                                               cTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCt                                        >  1:265828/1‑149 (MQ=255)                                                                                                                          ccTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCa                                                    <  1:154090/126‑1 (MQ=255)                                                                                                                                gtgACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTg                        >  6:270321/1‑148 (MQ=255)                                                                                                                                     acCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACggg                  >  5:382476/1‑149 (MQ=255)                                                                                                                                        ctcaaggTGGCGATCGTCGCCACCTTGAGGTCTTAGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCAACATCGGTACGGGCCa               >  8:195079/6‑149 (MQ=255)                                                                                                                                        ccgccgGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCa               <  5:256099/149‑1 (MQ=255)                                                                                                                                          gccgGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGc             >  8:387755/1‑149 (MQ=255)                                                                                                                                             ggTGGCGATCGTCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACg                             <  6:99786/130‑1 (MQ=255)                                                                                                                                             ggTGGCGATCGTCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACg                             >  5:99786/1‑130 (MQ=255)                                                                                                                                             ggTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGCttt          <  1:412576/149‑1 (MQ=255)                                                                                                                                                   gacgaTCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGCTTTCATCg     >  7:112968/1‑148 (MQ=255)                                                                                                                                                   gacgaTCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGCTTTCATCGc    <  5:270321/149‑1 (MQ=255)                                                                                                                                                   gacgaTCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGCTTTCATCGc    >  1:256407/1‑149 (MQ=255)                                                                                                                                                   gacgaTCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGCTTTCATCGc    <  1:226680/149‑1 (MQ=255)                                                                                                                                                    acgaTCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGGAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGCTTTCATCGCa   <  8:26624/149‑1 (MQ=255)                                                                                                                                                    acgaTCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGCTTTCATCGCa   <  3:265860/149‑1 (MQ=255)                                                                                                                                                     cgaTCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCAGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGCTTTCATCGCAg  <  1:356256/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CGCCGAGCGTGGCGTACTCGGCCAGTGTCATCTCGGTAAGGCGCTGCGGGTCGTTCACCAGCAGGTCAAGGAACTGCGCGTCCCACTCTGGGTTG‑TTGAAGCCGCAACGCTCGCCGTGCACCTGGTGTGACACCCCGCCGGTGGCGACGATCGCCACCTTGAGGTCTTCGGGGTAGCTCTCGATGGCGCGGCGCAATGCCTGGCCCAGTTTGTAGCAGCGCCGCGCCGTGGGTATCGGCAGCTGCAGTACGCCGACCTGCAGCGGCACGATCTGTACGGGCCAGCTTTCATCGCAG  >  NC_002947/2866074‑2866369

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