Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,846,345:1	+A	7.0%	coding (594/1377 nt)	pfeS‑II →	histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,846,345	1	.	A	7.0%	217.9 / 10.4	57	coding (594/1377 nt)	pfeS‑II	histidine kinase
Reads supporting (aligned to +/- strand):  ref base . (26/27);  new base A (2/2);  total (28/29)
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.09e-01

CAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAATGCGCACGCC  >  NC_002947/1846199‑1846494                                                                                                                                                      |                                                                                                                                                      cAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCg                                                                                                                                                    >  1:339633/1‑150 (MQ=255)   gCACCGCTGGCGGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGA‑CAGCTGATCGGc                                                                                                                                                  <  1:80130/150‑1 (MQ=255)     cccGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGccc                                                                                                                                                <  2:868525/149‑1 (MQ=255)     aCCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGA‑CAGCTGATCGGccc                                                                                                                                                <  1:1023370/150‑1 (MQ=255)      ccGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCg                                                                                                                                               >  2:846119/1‑150 (MQ=255)         cTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCt                                                                                                                                            <  1:50769/150‑1 (MQ=255)            gCTGTGGCGCCTTGTCACCAACGGCATCATCCCGGGGCTGTTCACATTGCT‑CAGT‑TGTCGACCTTCAGAGGACATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAACGCCTGCCGGGCCGATCAGCTG‑TCGGCCCGCCTGGa                                                                                                                                         <  1:684693/150‑1 (MQ=255)             cTGTGGCGCATTGTCACCAACGGCATCATCCCAGGGGTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGAc                                                                                                                                        <  1:555849/150‑1 (MQ=255)                  gcgcATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTcgcg                                                                                                                                   >  1:765316/1‑150 (MQ=255)                  gcgcATTGTCACAAAGGGCATCATCCCGGGGCTGTTCACATTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTcgcg                                                                                                                                   >  1:317396/1‑150 (MQ=255)                     cATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCAc                                                                                                                                <  2:828358/150‑1 (MQ=255)                          tCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAGCTTAGTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGATG‑AGGGCCCGCCGGGACTCCCGCACCATCg                                                                                                                           >  1:515903/1‑150 (MQ=255)                            aCCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGcc                                                                                                                         <  2:339633/150‑1 (MQ=255)                             ccAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGccc                                                                                                                        >  2:719802/1‑150 (MQ=255)                             ccAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGccc                                                                                                                        <  1:610341/150‑1 (MQ=255)                                 cGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCa                                                                                                                    <  1:499080/150‑1 (MQ=255)                                  ggCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGCCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCAt                                                                                                                   <  1:814971/150‑1 (MQ=255)                                      tcatcCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACg                                                                                                               >  1:322999/1‑150 (MQ=255)                                          cccGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCa                                                                                                                    >  2:413262/1‑141 (MQ=255)                                          cccGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCa                                                                                                                    <  1:413262/141‑1 (MQ=255)                                           ccGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTg                                                                                                          >  1:514132/1‑150 (MQ=255)                                            cGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGCGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTgg                                                                                                         <  1:74500/150‑1 (MQ=255)                                            cGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTgg                                                                                                         >  1:469111/1‑150 (MQ=255)                                                       aCCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTggg                                                                                                        <  1:81639/140‑1 (MQ=255)                                                       aCCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTggg                                                                                                        >  2:81639/1‑140 (MQ=255)                                                       aCCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGgcg                                                                                              >  2:791060/1‑150 (MQ=255)                                                               gttgtGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCtt                                                                                      <  2:322999/150‑1 (MQ=255)                                                                 tgtgtgTAGGCCTGTACCG‑‑CAGGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTg                                                                                    <  1:476438/150‑1 (MQ=255)                                                                     tgtAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCCACCGCGGCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCa                                                                                >  1:674496/1‑150 (MQ=255)                                                                      gtAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAg                                                                               >  2:995304/1‑150 (MQ=255)                                                                       tAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGa                                                                              >  2:391104/1‑150 (MQ=255)                                                                            cTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGcc                                                                         >  1:22857/1‑150 (MQ=255)                                                                             tGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTTGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TAGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCg                                                                        >  1:532864/1‑150 (MQ=255)                                                                                       aTGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTg                                                                <  1:846119/150‑1 (MQ=255)                                                                                        tGCTGATCGTGCCGCTGAACCAACTTCGCGGGCAGGCCGTTGCGTGGCGGGCGTATCAGCTG‑TCGGCCCTCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGc                                                               <  2:514132/150‑1 (MQ=255)                                                                                          cTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGACtctaa                                                             <  2:593348/150‑5 (MQ=255)                                                                                                 ggCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGTCGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACgg                                                      <  2:189699/149‑1 (MQ=255)                                                                                                  gCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGt                                                     <  1:538409/150‑1 (MQ=255)                                                                                                             aaCTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGaa                                                                                                    >  2:500710/1‑92 (MQ=255)                                                                                                             aaCTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGaa                                                                                                    <  1:500710/92‑1 (MQ=255)                                                                                                              aCTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATgcagca                                         >  2:171701/1‑150 (MQ=255)                                                                                                              aCTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATgcagca                                         >  2:499101/1‑150 (MQ=255)                                                                                                                   gcgAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGt                                    <  1:995304/150‑1 (MQ=255)                                                                                                                   gcgAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGt                                    <  2:515903/150‑1 (MQ=255)                                                                                                                         aGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGTGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTtgctgc                              <  1:791060/150‑1 (MQ=255)                                                                                                                           gCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTggg                                                                                                        >  2:294252/1‑74 (MQ=255)                                                                                                                           gCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTggg                                                                                                        <  1:294252/74‑1 (MQ=255)                                                                                                                           gCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTgcgc                            >  1:28408/1‑150 (MQ=255)                                                                                                                           gCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTgcgc                            >  1:18789/1‑150 (MQ=255)                                                                                                                             cAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGc                                                                     >  2:1085435/1‑107 (MQ=255)                                                                                                                             cAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGc                                                                     <  1:1085435/107‑1 (MQ=255)                                                                                                                                 tccaggCGGGCCGA‑CAGCTGATCGGCCCGCCAGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCATGGGACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTg                      >  1:10739/5‑150 (MQ=255)                                                                                                                                  ccTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGt                     <  1:171701/150‑1 (MQ=255)                                                                                                                                        gggCCGA‑CAGCTGATCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACg               >  2:306582/1‑150 (MQ=255)                                                                                                                                              aTCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGTTTGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAATgc         <  2:469111/150‑1 (MQ=255)                                                                                                                                               tCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAATgcg        >  1:756000/1‑150 (MQ=255)                                                                                                                                                 aGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGCATTGGGCGAACCGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGGAACGCGCa      >  2:963008/1‑150 (MQ=255)                                                                                                                                                     g‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAATGCGCAcgcc  >  1:685269/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CAGCACCGCTTGCTGTGGCGCATTGTCACCAACGGCATCATCCCGGGGCTGTTCACCTTGCTGTTGTGTGTAGGCCTGTACCG‑‑CATGCTGATCGTGCCGCTGAACCAACTTCGCGAGCAGGCCAATGCCTGGCGGGCCGATCAGCTG‑TCGGCCCGCCTGGACTCGCGCACCATCGCCCGGCATGACGAATTGGGCGAACTGGCGCGTGCCTTTGACCAGATGGCCGAGCGGCTGCAAGGCACGGTGGCTATGCAGCAGCAGTTGCTGCGCGACCTGTCCCACGAAATGCGCACGCC  >  NC_002947/1846199‑1846494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑