Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,687,693	C→G	6.5%	I223M (ATC→ATG)	PP_4148 →	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,687,693	0	C	G	6.5%	189.0 / 2.9	61	I223M (ATC→ATG)	PP_4148	ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (2/2);  total (31/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 = 4.36e-01

GTTCGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACGT  >  NC_002947/4687549‑4687830                                                                                                                                                 |                                                                                                                                           gTTCGCCATCCTGCTGATCGTGCTGTTCACCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCCCCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACc‑g                                                                                                                                       >  1:23528/1‑149 (MQ=255)  ttCGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACc‑gc                                                                                                                                      >  6:228487/1‑149 (MQ=255)       cATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTgg                                                                                                                                 >  5:144668/1‑149 (MQ=255)        aTCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGt                                                                                                                                <  7:355541/149‑1 (MQ=255)                     tGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTc                                                                                                                   <  5:168317/149‑1 (MQ=255)                     tGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTc                                                                                                                   >  3:256800/1‑149 (MQ=255)                      gCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCg                                                                                                                  >  8:7321/1‑149 (MQ=255)                       cTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGc                                                                                                                 <  1:299883/149‑1 (MQ=255)                         gTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGc                                                                                                                         <  2:276147/139‑1 (MQ=39)                         gTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGc                                                                                                                         >  1:276147/1‑139 (MQ=39)                            cGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGccacc                                                                                                             <  5:18758/148‑1 (MQ=255)                              ccGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAGCAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGCa                                                                                                                        >  1:292963/1‑135 (MQ=255)                              ccGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAGCAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATGACCAGCGC‑GGTGATCGGCa                                                                                                                        <  2:292963/135‑1 (MQ=255)                                           cTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGAcc                                                                                              <  1:6223/148‑1 (MQ=255)                                                 cTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACa                                                                                       <  5:54927/149‑1 (MQ=255)                                                  tGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTc                                                                                                                               <  8:308861/108‑1 (MQ=255)                                                  tGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTc                                                                                                                               >  7:308861/1‑108 (MQ=255)                                                       cccGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTc                                                                                  <  6:371925/148‑1 (MQ=255)                                                       cccGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTc                                                                                  <  5:169877/148‑1 (MQ=255)                                                         cGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTc                                                                               >  2:156141/1‑149 (MQ=255)                                                          gCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGc                                                                                                                 <  2:48539/114‑1 (MQ=255)                                                          gCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGc                                                                                                                 >  1:48539/1‑114 (MQ=255)                                                          gCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTc                                                                               >  1:100744/1‑148 (MQ=255)                                                                    atCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGACGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGAACCTGCTGACAAAGAACAGCTTCCTCGACGAAATcaa                                                                    >  1:242205/2‑149 (MQ=255)                                                                     tCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATcaa                                                                    >  3:177477/1‑148 (MQ=255)                                                                        ccTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCaacaaa                                                                <  2:100744/149‑1 (MQ=255)                                                                         cTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCaacaaa                                                                >  7:108345/1‑148 (MQ=255)                                                                           ccAACAACTTCGACGAACTCAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAg                                                             <  4:307969/149‑1 (MQ=255)                                                                            caacaaCTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTAATCGGCAACTTATCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGt                                                            <  5:228996/149‑1 (MQ=255)                                                                             aacaacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTa                                                           >  6:37469/1‑149 (MQ=255)                                                                             aacaacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTa                                                           <  8:108345/149‑1 (MQ=255)                                                                               caacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACg                                                         >  3:345372/1‑149 (MQ=255)                                                                                aacTTCGACGCACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTAGACACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAAGCAACAACCAGTACGt                                                        >  6:279898/1‑149 (MQ=39)                                                                                aacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCAGCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGt                                                        <  5:228487/149‑1 (MQ=255)                                                                                aacTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGt                                                        >  8:301187/1‑149 (MQ=255)                                                                                 acTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTg                                                       <  3:155424/149‑1 (MQ=255)                                                                                   ttCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGt                                                     >  1:146550/1‑149 (MQ=255)                                                                                   ttCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGt                                                     >  6:193/1‑149 (MQ=255)                                                                                      gacgaACTTAGCACCACCTGCATGGTGCTGGACTACTTCTGGCACGTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAATAACAGCTTCGTCGACGAAATCAACAAACAGTACGTGGTCAc                                                  <  5:279898/149‑1 (MQ=255)                                                                                          aaCTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTTACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTTGTCACCgcc                                              <  4:177477/149‑1 (MQ=255)                                                                                          aaCTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgc                                               >  4:371648/1‑148 (MQ=255)                                                                                               aGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggc                                         >  3:367183/1‑149 (MQ=255)                                                                                                  accaccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCgccaaggcc                                        >  8:205464/1‑147 (MQ=255)                                                                                                    caccGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGccaaggccaagg                                    >  5:203673/1‑149 (MQ=255)                                                                                                    caccGGCAAGGTGCTGGACTACTTCTGGCACATGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGATATCAACAAACAGTACGTGGTCACCCccaaggccaagg                                    >  6:283744/1‑149 (MQ=255)                                                                                                        ggCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCt                                <  1:66033/149‑1 (MQ=255)                                                                                                         gCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTTGTCACCGCCAAGGCCAAGGGCCTg                               <  6:203673/149‑1 (MQ=255)                                                                                                            aGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGc                            <  4:367183/149‑1 (MQ=255)                                                                                                             ggTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGcc                           <  7:163670/149‑1 (MQ=255)                                                                                                              gTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGcc                           >  8:86352/1‑148 (MQ=255)                                                                                                                  tGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCa                      >  3:291719/1‑149 (MQ=255)                                                                                                                    gACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACg                    >  7:286037/1‑149 (MQ=255)                                                                                                                     actactTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACgg                   <  2:23528/149‑1 (MQ=255)                                                                                                                     actactTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACgg                   <  5:193/149‑1 (MQ=255)                                                                                                                          ttCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACAAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCt              >  7:317586/1‑149 (MQ=255)                                                                                                                              ggCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTAc          >  5:181231/1‑149 (MQ=255)                                                                                                                               gCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTAc          <  3:371648/148‑1 (MQ=255)                                                                                                                                      gTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCAc    <  6:181231/147‑1 (MQ=255)                                                                                                                                      gTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCAACATGACCCTGCTGACCTAGAACAGCTTCGTCGACGAAATCAACATTCAGTACGGTGTCACCGCCAAGGCCAAGGTCCTTAGCCTATCACGGGTTCTGTACTGCCACgt  <  2:242205/149‑1 (MQ=255)                                                                                                                                             cGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACgt  <  4:126207/142‑1 (MQ=255)                                                                                                                                             cGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACgt  >  3:126207/1‑142 (MQ=255)                                                                                                                                                 |                                                                                                                                           GTTCGCCATCCTGCTGATCGTGCTGTTCGCCGGTGGCAGCTACTTCGACTGGTTCCCGCTGCGTGGCCTCACCTCCAACAACTTCGACGAACTGAGCACCACCGGCAAGGTGCTGGACTACTTCTGGCACCTGGTGCTGCCGATCACC‑GCGCTGGTCATCGGCAACTTCGCCACCATGACCCTGCTGACCAAGAACAGCTTCCTCGACGAAATCAACAAACAGTACGTGGTCACCGCCAAGGCCAAGGGCCTGAGCCGATCACGGGTGCTGTACGGCCACGT  >  NC_002947/4687549‑4687830

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