Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,951,650	G→A	10.4%	L455L (CTG→TTG)	ppk ←	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,951,650	0	G	A	10.4%	145.0 / 5.9	57	L455L (CTG→TTG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (20/29);  new base A (2/4);  total (24/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 = 9.96e-01

GGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGG  >  NC_002947/5951506‑5951792                                                                                                                                                  |                                                                                                                                                gggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAggcggc                                                                                                                                            >  2:278542/1‑150 (MQ=255)   gcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTg                                                                                                                                          <  1:140713/150‑1 (MQ=255)   gcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCCCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACCCCG‑TAGATCACTACCACACCAGCCGCCTGCAAGCGGCTg                                                                                                                                          <  2:688931/150‑1 (MQ=255)     gTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGTGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTggt                                                                                                                                        <  2:976941/150‑2 (MQ=255)     gTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTggt                                                                                                                                        <  2:491768/150‑2 (MQ=255)          cGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCt                                                                                                                                   <  1:718951/150‑1 (MQ=255)             cGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCC‑GCCACGGTGTTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCa                                                                                                                                >  2:749743/1‑150 (MQ=255)                 gTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGtt                                                                                                                            <  1:102043/150‑1 (MQ=255)                    gTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGa                                                                                                                         <  1:943415/150‑1 (MQ=255)                    gTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGa                                                                                                                         <  2:741845/150‑1 (MQ=255)                        ccGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGActct                                                                                                                     >  2:132235/1‑150 (MQ=255)                         cGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCtt                                                                                                                    >  1:1410015/1‑150 (MQ=255)                          ggTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTc                                                                                                                   >  1:169311/1‑150 (MQ=255)                          ggTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTc                                                                                                                   <  2:169311/150‑1 (MQ=255)                              cccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                             >  2:276831/1‑120 (MQ=255)                              cccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                             <  1:276831/120‑1 (MQ=255)                               ccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCaa                                                                                                              >  1:504493/1‑150 (MQ=255)                                  gATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAAc                                                                                                           <  1:696381/150‑1 (MQ=255)                                    tGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCATGCCTCAGGTTGGACTCTTCGTCAAAAcgc                                                                                                         <  1:749743/150‑1 (MQ=255)                                      cgcgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACgcgc                                                                                                       <  1:132235/150‑1 (MQ=255)                                        cgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAg                                                                                                                               >  2:1154475/1‑124 (MQ=255)                                        cgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAg                                                                                                                               <  1:1154475/124‑1 (MQ=255)                                         gTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGcgcacg                                                                                                    <  2:746656/150‑1 (MQ=255)                                              gCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACt                                                                                               >  1:894380/1‑150 (MQ=255)                                                      tcGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑Acc                                                                                      >  2:213238/1‑150 (MQ=255)                                                         cccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCa                                                                                   <  2:894380/150‑1 (MQ=255)                                                         cccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCa                                                                                   <  2:100192/150‑1 (MQ=255)                                                         cccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCa                                                                                   <  1:487199/150‑1 (MQ=255)                                                             gCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGa                                                                               <  1:838825/150‑1 (MQ=255)                                                             gCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGa                                                                               <  1:863911/150‑1 (MQ=255)                                                               tCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTCGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCATGAGGCTGGCTATCTGCAGGTTGGACTCTTCGTCAAAACGCGCCCGCAACTCGATC‑ACCGCAGTGccc                                                                             >  2:267906/1‑150 (MQ=255)                                                                      gcAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCAAACCACACCG‑TAAATCACTACCGCACCAGCAGGGTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑AAAGCAGTGACCTCCTTGc                                                                      >  1:1168754/1‑150 (MQ=255)                                                                          aaTCAGCATTATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGtt                                                                  <  1:61560/150‑1 (MQ=255)                                                                           aTCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTa                                                                 <  1:278542/150‑1 (MQ=255)                                                                            tCAGCATCATCTTGGCGTGGGTCTTGAAGGCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTAc                                                                >  2:433544/1‑150 (MQ=255)                                                                                catcatCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACgcgc                                                            <  1:309237/150‑1 (MQ=255)                                                                                     tcTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAt                                                       >  1:814454/1‑150 (MQ=255)                                                                                            gTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCNANTCGATC‑NCCGCAGTGACCTNCTTGCCGTTACGCGCCGCATCCACNAg                                                <  1:881755/150‑1 (MQ=255)                                                                                             tGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAgg                                               <  1:213238/150‑1 (MQ=255)                                                                                                   ttGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTcg                                         >  2:793677/1‑150 (MQ=255)                                                                                                   ttGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTcg                                         >  2:192573/1‑150 (MQ=255)                                                                                                            aCCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCg                                                                                                                  <  2:427431/69‑1 (MQ=255)                                                                                                            aCCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCg                                                                                                                  >  1:427431/1‑69 (MQ=255)                                                                                                              cacaCCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCaa                                                                                                 <  1:1010956/84‑1 (MQ=255)                                                                                                              cacaCCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCaa                                                                                                 >  2:1010956/1‑84 (MQ=255)                                                                                                              cacaCCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGa                              <  1:295357/150‑1 (MQ=255)                                                                                                              cacaCCG‑TAGATCACTACCGCACCAGCCACATTCTTTGGGCATTCTTGCTGCATGTTGGTATCTTCGTCAACACGCG‑GGGGTGCATGCTCAACTGCACTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGAGTCGACGTGTTCCGa                              >  1:1033789/1‑150 (MQ=255)                                                                                                                   cG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCACCGGCTAGTGCAGC‑AGCCCAGAGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAg                         >  2:1046804/1‑150 (MQ=255)                                                                                                                         aTCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGccccc                    >  2:959094/1‑150 (MQ=255)                                                                                                                            aCTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGcc                                                                     <  1:573515/98‑1 (MQ=255)                                                                                                                            aCTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGcc                                                                     >  2:573515/1‑98 (MQ=255)                                                                                                                            aCTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAg                 <  2:1410015/150‑1 (MQ=255)                                                                                                                             cTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGc                <  1:433544/150‑1 (MQ=255)                                                                                                                              tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTGGGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCg               >  2:911030/1‑150 (MQ=255)                                                                                                                                  tggccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTAc           >  1:304456/4‑150 (MQ=255)                                                                                                                                   ggccAGCCGCCTGCAAGCGGCTGGTCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTcc                                >  2:304330/3‑128 (MQ=255)                                                                                                                                   ggccAGCCGCCTGCAAGCGGCTGGTCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTcc                                <  1:304330/126‑1 (MQ=255)                                                                                                                                      cAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAggg       >  1:1424612/1‑150 (MQ=255)                                                                                                                                         ccGCTTGCAGGCGGCTGGCCATCTGCAGGGGGGGCTTTTCGTCAAAACGTGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCt    <  2:468460/150‑1 (MQ=255)                                                                                                                                           gCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACGGGGTCTgg  <  1:192573/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                GGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGG  >  NC_002947/5951506‑5951792

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