Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,312,690	A→C	6.5%	intergenic (+31/‑131)	PP_2912 → / → hemB	hypothetical protein/5‑aminolevulinate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,312,690	0	A	C	6.5%	158.4 / 2.3	59	intergenic (+31/‑131)	PP_2912/hemB	hypothetical protein/5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (23/31);  new base C (2/2);  total (26/33)
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.31e-01

CGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAA  >  NC_002947/3312549‑3312832                                                                                                                                              |                                                                                                                                                cGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGa                                                                                                                                         <  4:148775/149‑1 (MQ=255)  gTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGACCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGATTATTCCTGTACCGCCTATTCCCCGGGTTA‑ACCGCGa                                                                                                                                         >  8:297867/1‑148 (MQ=255)     agaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAAC‑CCCGCGAagag                                                                                                                                     <  7:282315/149‑1 (MQ=255)     agaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAagag                                                                                                                                     <  1:329221/149‑1 (MQ=255)          gACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCag                                                                                                                                >  3:290074/1‑148 (MQ=255)           aCGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAac                                                                                                                               <  1:46662/149‑1 (MQ=255)              ccgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAAcaca                                                                                                                            <  2:66409/149‑1 (MQ=255)               cgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAg                                                                                                                           <  5:260668/149‑1 (MQ=255)               cgccCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAg                                                                                                                           <  1:262449/149‑1 (MQ=255)                  ccGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAAC‑CCCGCGAaga                                                                                                                                      >  5:249269/1‑135 (MQ=255)                  ccGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGTAAC‑CCCGCGAaga                                                                                                                                      <  6:249269/135‑1 (MQ=255)                   cGGATGGACGACGCTCGTCGCCATGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGGTC‑CCCGCGAAGAGGCCAACACAGGCa                                                                                                                        >  7:123510/1‑148 (MQ=255)                         gacgacGCTCGTCGCCGTGCAGCACGGGCCGAGGCCGAGCATGGGCCGGCGAAGGCAGAGGAGGTAACGCCGCCCTTGATGCAGTAGTTTATTCCTGTACCGCCTCTTCGGGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACg                                                                                                                 <  7:297867/149‑1 (MQ=255)                         gacgacGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACg                                                                                                                 <  5:127509/149‑1 (MQ=255)                           cgacgCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGc                                                                                                               >  1:52878/1‑149 (MQ=255)                            gacgCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGcc                                                                                                              >  7:79114/1‑149 (MQ=255)                             acgCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCg                                                                                                             <  4:284285/149‑1 (MQ=255)                                 tcgtcgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGc                                                                                                         >  5:85282/1‑149 (MQ=255)                                 tcgtcgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGc                                                                                                         >  8:96314/1‑149 (MQ=255)                                   gtcgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGcc                                                                                                        >  8:223228/1‑148 (MQ=255)                                   gtcgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGc                                                                                                         <  1:301875/147‑1 (MQ=255)                                       ccGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGa                                                                                                   <  3:121077/149‑1 (MQ=255)                                        cGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGAc                                                                                                  <  2:52878/149‑1 (MQ=255)                                        cGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGGTTATTCCTGTACCGCCTCTTCGCGGGGGTA‑CCCGCGAAGAGGGCAACACAGGCATACCCCGGCCCTTGCAATTg                                                                                                    >  1:333470/1‑147 (MQ=255)                                          tGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACg                                                                                                 >  3:71785/1‑148 (MQ=255)                                          tGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGc                                                                                                <  3:145090/149‑1 (MQ=255)                                               caCGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGac                                                                                           <  3:278021/149‑1 (MQ=255)                                                  gcgcCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGAcaca                                                                                        <  7:96314/149‑1 (MQ=255)                                                          gccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCAtt                                                                                                     >  4:151281/1‑128 (MQ=255)                                                          gccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGt                                                                                 >  1:17186/1‑148 (MQ=255)                                                               gCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTa                                                                           <  7:75656/149‑1 (MQ=255)                                                                cATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAg                                                                          <  6:85282/149‑1 (MQ=255)                                                                         gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                   <  4:224338/131‑1 (MQ=255)                                                                         gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                   >  3:224338/1‑131 (MQ=255)                                                                              ggCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGGTA‑GCCGGGAAGAGGCCAACACAGGCAAACCACGGGAGTTGGCATTGTCGGAAGGCACATGGTCGTCAGTTAGGCTTTTCAGTGgcg                                                            >  8:282604/1‑149 (MQ=255)                                                                               gCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAGCACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGgcgc                                                           <  4:290074/149‑1 (MQ=255)                                                                                 agagGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGATGAGGCCAAAACAGGCCAAACACGGCCGTTGGCATAGACGCAAAACACCATAACCGCAATTAGGCTTTTCACTGgccccc                                                         >  2:145963/1‑145 (MQ=255)                                                                                        gTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTAGCCGGGTTT‑CCCCCGAAGAGGCCAACACACGCAAAACACGGCCCTTGGCATTTACGCAAGACACAATCTAGTCAATTATCCATTTCACTGGCGCGCGCCacg                                                   >  2:158569/1‑146 (MQ=255)                                                                                             gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTAGCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCAATGACGCAAGACACATTCTCGTCATTTAGCCTTTTAACTGGCGCGCGCCAGGCAt                                                >  1:255061/1‑147 (MQ=255)                                                                                                 cccTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACTGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGtat                                         <  4:71785/149‑1 (MQ=255)                                                                                                 cccTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGta                                          >  6:251436/1‑148 (MQ=255)                                                                                                            gtagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGc                                                  <  8:195196/129‑1 (MQ=255)                                                                                                            gtagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGc                                                  >  7:195196/1‑129 (MQ=255)                                                                                                             tagtTTATTCCTGTACCGCCTCTTAGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGGTGCCATTGACGCAAGACACATTATCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATAAGCCGTTCaaa                             >  2:157388/1‑149 (MQ=255)                                                                                                                tttATTCCTGTACCGCCTCTTCGCGGGGTTA‑ACCGCGAAGAGGCCCACACAGGCAAACCACGGCCCATGCCATAGACGCAAGACACAATCTCGTCAGTTTGACTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTc                          >  1:177607/1‑149 (MQ=255)                                                                                                                  tATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGt                                           >  2:296265/1‑130 (MQ=255)                                                                                                                  tATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGt                                           <  1:296265/130‑1 (MQ=255)                                                                                                                    ttCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTg                      <  5:251436/149‑1 (MQ=255)                                                                                                                       cTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATACCCTGGAg                   <  8:79114/149‑1 (MQ=255)                                                                                                                       cTGTACCGCCTCTTAGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTATCACTGGCGCGCGACAGGCATAAGGTATACGACGTTCAAATTACCtag                     >  6:183341/1‑145 (MQ=255)                                                                                                                          ttccGCCTCTTCGCGGGTGTA‑CCAGAGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTAAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGTAGCt                 <  8:146791/146‑1 (MQ=255)                                                                                                                            ccGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTc              <  1:158569/149‑1 (MQ=255)                                                                                                                              gCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCAGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCAt            <  7:282604/149‑1 (MQ=255)                                                                                                                                ctctTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGt          <  2:17186/149‑1 (MQ=255)                                                                                                                                 tctTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGtat                                         >  7:256519/1‑117 (MQ=255)                                                                                                                                 tctTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGtat                                         <  8:256519/117‑1 (MQ=255)                                                                                                                                        ggggTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCaaa                             <  6:61801/122‑1 (MQ=255)                                                                                                                                        ggggTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCaaa                             >  5:61801/1‑122 (MQ=255)                                                                                                                                        ggggTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaaccaa  <  7:223228/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                CGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAA  >  NC_002947/3312549‑3312832

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