Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,649,996	C→T	10.5%	P306P (CCG→CCA)	PP_4113 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,649,996	0	C	T	10.5%	146.9 / 7.5	57	P306P (CCG→CCA)	PP_4113	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base T (2/4);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

ACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTCAAAGTCCA  >  NC_002947/4649853‑4650144                                                                                                                                                |                                                                                                                                                     aCCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCa                                                                                                                                                  >  6:173780/1‑148 (MQ=255) aCCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTAACCGAGGCAGCGATGTAGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGCCTGTGCTTGCCAAGCCCCGGCa                                                                                                                                                  >  8:91/1‑148 (MQ=255)    cGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCGGGGCTGGGCAAGCACAGCCCTGGcttcga                                                                                                                                              <  6:207499/149‑6 (MQ=255)     ggTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCGGGGCTGGGCAAGCACAGCCCTGGcttcgag                                                                                                                                             <  3:129631/149‑7 (MQ=255)           tGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCa                                                                                                                                       <  2:47867/149‑1 (MQ=255)           tGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCa                                                                                                                                       <  1:133706/149‑1 (MQ=255)              cccAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCAccc                                                                                                                                    >  2:173651/1‑149 (MQ=255)              cccAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCAccc                                                                                                                                    >  2:302003/1‑149 (MQ=255)              cccAGCTCCAGTTCAGCACCCGCACGCCATCATCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGcc                                                                                                                                           >  1:89483/1‑142 (MQ=255)              cccAGCTCCAGTTCAGCACCCGCACGCCATCATCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGcc                                                                                                                                           <  2:89483/142‑1 (MQ=255)                cAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTg                                                                                                                                  <  3:159277/149‑1 (MQ=255)                  gCTCCAGTGCAGCCCGCGCACTCCAGCTTCGACAAGGTTCACCGAGGGAGGGTGGTTGGCGTTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCGCGGCAAAAAGCCCTCACCCTgcg                                                                                                                                <  3:238209/149‑1 (MQ=255)                      cAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGc                                                                                                                            >  3:9101/1‑149 (MQ=255)                      cAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGc                                                                                                                            <  3:301752/149‑1 (MQ=255)                          tCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCgag                                                                                                                         >  3:13191/1‑148 (MQ=255)                                 ccGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                                 <  1:173651/149‑1 (MQ=255)                                 ccGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTg                                                                                                                 >  1:68635/1‑149 (MQ=255)                                       gCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCg                                                                                                                           >  3:105950/1‑133 (MQ=255)                                       gCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCg                                                                                                                           <  4:105950/133‑1 (MQ=255)                                         cATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATgg                                                                                                         <  1:302003/149‑1 (MQ=255)                                            cTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTg                                                                                                      >  3:69329/1‑149 (MQ=255)                                                 aCCAGGTTCTCCGAGGGAGCGATGTTGGCGGTGATACCGGCTTCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCTGCAAAAAGCGCTCACCCTGCGGAGCGAGAATTCCTGCGACATGGCTGCCGTg                                                                                                 <  7:171690/149‑1 (MQ=255)                                                  ccAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACAATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTgg                                                                                                >  2:264242/1‑149 (MQ=255)                                                              aGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCt                                                                                                                                         <  1:16039/96‑1 (MQ=255)                                                              aGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCt                                                                                                                                         >  2:16039/1‑96 (MQ=255)                                                                 cAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTa                                                                                 <  2:68635/149‑1 (MQ=255)                                                                                 gATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCCGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTaga                                                                 <  7:308622/149‑1 (MQ=255)                                                                                      cGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGc                                                                                               >  2:291314/1‑114 (MQ=255)                                                                                      cGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCAt                                                            >  5:160140/1‑149 (MQ=255)                                                                                           tCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTg                                                       <  1:264242/149‑1 (MQ=255)                                                                                                         tACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCa                                           <  5:173780/147‑1 (MQ=255)                                                                                                           ccATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAg                                          <  6:160140/146‑1 (MQ=255)                                                                                                           ccATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGcc                                        <  4:171013/148‑1 (MQ=255)                                                                                                                aCCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTc                                   >  4:131633/1‑148 (MQ=255)                                                                                                                aCCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCg                                  <  4:9101/149‑1 (MQ=255)                                                                                                                 ccTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGa                                 >  3:195323/1‑149 (MQ=255)                                                                                                                    cGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGCGCGGAGAGGCAGGTCTGGCCTTGTTCGCAGCCCCCTCGATaa                              >  4:201510/1‑149 (MQ=255)                                                                                                                     gAAGCAAGGGATGTGCTTGCCCATCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTTCCGGGCGGTTACCGTCGCGAGCGTAGAGGAATGTCTTGCCTTGTTCGCAGCCCCCTCGATAAt                             <  7:91/149‑1 (MQ=255)                                                                                                                       aGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATgg                           <  4:13191/149‑1 (MQ=255)                                                                                                                        gCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                              <  7:188988/113‑1 (MQ=255)                                                                                                                        gCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                              >  8:188988/1‑113 (MQ=255)                                                                                                                         ccAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt                         >  1:278336/1‑149 (MQ=255)                                                                                                                          cggggCTGGGCAAGCACAGCCGCGGCAATAAGCCCTCACGCTGCGCAGCGAGAATGCCTGCGAGATGGCTGGCGTGGCTTGCCGGGCGGTTACGGTCGCGTGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt                         <  8:154078/146‑1 (MQ=255)                                                                                                                          cggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                              <  3:115850/109‑1 (MQ=255)                                                                                                                          cggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                              >  4:115850/3‑111 (MQ=255)                                                                                                                          cggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt                         >  4:89813/3‑148 (MQ=255)                                                                                                                          cggggCTGGGCAAGCACAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATCGCTGCCGTGGCATGCCGGGCGGGTACCGGCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCt                         >  7:154078/3‑148 (MQ=255)                                                                                                                           ggggCTGGGCAAGCACAGCCCTGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGa                       <  3:89813/148‑1 (MQ=255)                                                                                                                             ggCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAAc                     >  4:121119/1‑149 (MQ=255)                                                                                                                                tgtgCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAAcccc                  >  7:320851/1‑149 (MQ=255)                                                                                                                                  tgCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAAcccc                  <  2:278336/147‑1 (MQ=255)                                                                                                                                    cTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                              <  6:85048/101‑1 (MQ=255)                                                                                                                                    cTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGc                                                              >  5:85048/1‑101 (MQ=255)                                                                                                                                    cTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCggg               <  4:195323/148‑1 (MQ=255)                                                                                                                                       gCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCg            >  4:305085/1‑148 (MQ=255)                                                                                                                                       gCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCg            <  3:305085/148‑1 (MQ=255)                                                                                                                                        cccAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGt           >  1:328319/1‑148 (MQ=255)                                                                                                                                           agacccGGCAAAAATCCGTCTACCTGCGCAGCGAGATTGCCTGCGACATTGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTCaaa       <  8:292459/146‑1 (MQ=255)                                                                                                                                                cGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTCAAAGTCCa  >  1:156888/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     ACCCGGTGAATGCCCCAGCTCCAGTTCAGCACCCGCACGCCATCTTCGACCAGGTTCACCGAGGCAGCGATGTTGGCGGTGATACCGGCATCCGAGTTGCGCTCTACCATCACCTCGAAGCCAGGGCTGTGCTTGCCCAGCCCCGGCAAAAAGCCCTCACCCTGCGCAGCGAGAATGCCTGCGACATGGCTGCCGTGGCTTGCCGGGCGGTTACCGTCGCGAGCGTAGAGGCATGTCTGGCCTTGTTCGCAGCCCCCTCGATAATGGCTGAACCCCGGGGCGTCAAAGTCCA  >  NC_002947/4649853‑4650144

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