Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,549,922	C→G	8.2%	H418Q (CAC→CAG) ‡	PP_3133 →	oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,549,922	0	C	G	8.2%	182.3 / 4.7	59	H418Q (CAC→CAG) ‡	PP_3133	oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (33/21);  new base G (3/2);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

CTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTGCAA  >  NC_002947/3549776‑3550069                                                                                                                                                   |                                                                                                                                                    cTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATcaca                                                                                                                                                    <  1:13441/148‑1 (MQ=255)         cGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcatgg                                                                                                                                           <  2:355521/149‑4 (MQ=255)         cGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGcacca                                                                                                                                           >  7:307151/1‑149 (MQ=255)         cGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGcacca                                                                                                                                           >  3:121678/1‑149 (MQ=255)                 ggCACACCGTCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATg                                                                                                                                   <  4:74612/149‑1 (MQ=255)                  gCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATgg                                                                                                                                  >  1:345211/1‑149 (MQ=255)                   cacaCCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGgc                                                                                                                                 <  2:265288/149‑1 (MQ=255)                    acacCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCAGATGGGCACCACGCGCATGgcc                                                                                                                                >  4:346244/1‑149 (MQ=255)                        cATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCt                                                                                                                            >  7:207861/1‑149 (MQ=255)                            aaGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTcgc                                                                                                                         <  6:203994/148‑1 (MQ=255)                                         gCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGgt                                                                                                              <  1:33660/146‑1 (MQ=255)                                         gCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGgt                                                                                                              >  2:33660/1‑146 (MQ=255)                                            aaCTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGt                                                                                                         >  6:160825/1‑148 (MQ=255)                                                 gCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                   <  1:288379/149‑1 (MQ=255)                                                 gCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                   >  1:38558/1‑149 (MQ=255)                                                 gCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                   <  7:173306/149‑1 (MQ=255)                                                   cAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGcc                                                                                                  >  3:218571/1‑148 (MQ=255)                                                   cAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                   <  7:202861/147‑1 (MQ=255)                                                   cAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                   >  8:202861/1‑147 (MQ=255)                                                       cAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcatggg                                                                                                                                          >  3:167491/1‑100 (MQ=255)                                                       cAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcatggg                                                                                                                                          <  4:167491/104‑5 (MQ=255)                                                         gTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCa                                                                                           <  2:345211/149‑1 (MQ=255)                                                            cGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGCCACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAgg                                                                                        >  6:21391/1‑149 (MQ=255)                                                            cGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAgg                                                                                        >  3:55159/1‑149 (MQ=255)                                                                 aCATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCaa                                                                                                <  5:260890/136‑1 (MQ=255)                                                                 aCATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGCCACCCCATGACCATCACATCGGCACCACCCCCATCCCCGCCTCGACGCAGTTATGTGTGGACGATGACaa                                                                                                >  6:260890/1‑136 (MQ=255)                                                                         cTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCaa                                                                                          >  1:212388/1‑134 (MQ=255)                                                                         cTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCaa                                                                                          <  2:212388/134‑1 (MQ=255)                                                                           ggggTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGaa                                                                          >  3:100863/1‑148 (MQ=255)                                                                             ggTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAAc                                                                       <  5:361759/149‑1 (MQ=255)                                                                              gTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCACTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAAcc                                                                      >  2:279469/1‑149 (MQ=255)                                                                                 cGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTg                                                                    <  1:280057/148‑1 (MQ=255)                                                                                          cAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCg                                                          >  7:325579/1‑149 (MQ=255)                                                                                           accGATTTCGTCTCCGACACCCCGATACCGCTGAGAATGGCACCCCCTGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGcc                                                           <  8:191105/145‑1 (MQ=255)                                                                                             cGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGTGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCgg                                                         >  1:192273/1‑147 (MQ=255)                                                                                                  tCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGgcagcag                                                   >  6:101818/1‑148 (MQ=255)                                                                                                   cGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGgcagcagca                                                 >  4:102816/1‑149 (MQ=255)                                                                                                    gTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCAt                                                >  4:102541/1‑149 (MQ=255)                                                                                                     tttACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAg                                                                                         >  4:275643/1‑107 (MQ=255)                                                                                                     tttACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAg                                                                                         <  3:275643/107‑1 (MQ=255)                                                                                                     tttACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCAtt                                               >  3:358238/1‑149 (MQ=255)                                                                                                        aCGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCAtttt                                             <  8:207861/148‑1 (MQ=255)                                                                                                         cGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTg                                           >  2:349760/1‑149 (MQ=255)                                                                                                               cTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAgg                                     >  8:145882/1‑149 (MQ=255)                                                                                                               cTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAgg                                     >  3:361464/1‑149 (MQ=255)                                                                                                                 cGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGggc                                   >  6:86701/1‑149 (MQ=255)                                                                                                                       cGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGcc                             <  1:349760/149‑1 (MQ=255)                                                                                                                         cTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTc                           <  8:325579/149‑1 (MQ=255)                                                                                                                             gaATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTTCCAAGGGCGGCGCCTCCAAc                       <  4:55159/149‑1 (MQ=255)                                                                                                                             gaATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGACAAGGGCGGCGCCTCCAAc                       >  5:184234/1‑149 (MQ=255)                                                                                                                                tGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAAccc                     >  3:284430/1‑148 (MQ=255)                                                                                                                                  gCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGa                   <  5:86701/148‑1 (MQ=255)                                                                                                                                        ggtGCCCATGTGATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGACAATTGCAAGGTGTACGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCAAAGGGCGGCGCCTCCAACCCGACCATGcc            >  7:341782/3‑149 (MQ=255)                                                                                                                                        cATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCCATTGCAAGGTGTTCGGTACCGAAAACCTGTATCTCGCAGGCAGCAGCATTTTTGCCAACGGCGGCGCCTCCAACCCGACCATGcc            >  7:176488/1‑149 (MQ=255)                                                                                                                                         aTGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCg           <  2:38558/149‑1 (MQ=255)                                                                                                                                         aTGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGAATTTTTGCCAAGGGCGGTGCCTCCAACCCGACCATGCCg           >  3:22469/1‑149 (MQ=255)                                                                                                                                           gCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGAAGCATTTTTGCCAAGGGCGGCGCCTCCCACCCGACCATGCCg           >  8:147651/1‑147 (MQ=255)                                                                                                                                              cATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGtt      <  4:218571/149‑1 (MQ=255)                                                                                                                                              cATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGtt      >  4:350366/1‑149 (MQ=255)                                                                                                                                                tCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTg     >  2:3017/1‑148 (MQ=255)                                                                                                                                                  acaTGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTGCaa  >  7:94258/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTGCAA  >  NC_002947/3549776‑3550069

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