Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,951,645	T→C	11.5%	Q456Q (CAA→CAG)	ppk ←	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,951,645	0	T	C	11.5%	161.8 / 10.2	63	Q456Q (CAA→CAG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (23/31);  new base C (3/4);  total (28/35)
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.62e-01

TGTACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGG  >  NC_002947/5951497‑5951792                                                                                                                                                      |                                                                                                                                                     tgtACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCtt                                                                                                                                                     <  1:682638/150‑1 (MQ=255)   tACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCCTGc                                                                                                                                                   >  1:976941/1‑150 (MQ=255)          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                                                                                                                                             <  1:276831/120‑1 (MQ=255)                                       cccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                             >  2:276831/1‑120 (MQ=255)                                        ccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCa                                                                                                                                                  >  2:438921/1‑114 (MQ=255)                                        ccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCa                                                                                                                                                  <  1:438921/114‑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                                                                                                                               <  1:1154475/124‑1 (MQ=255)                                                 cgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAg                                                                                                                               >  2:1154475/1‑124 (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:863911/150‑1 (MQ=255)                                                                      gCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGa                                                                               <  1:838825/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:192573/1‑150 (MQ=255)                                                                                                            ttGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTcg                                         >  2:793677/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                                                                     >  2:573515/1‑98 (MQ=255)                                                                                                                                     aCTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGcc                                                                     <  1:573515/98‑1 (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)                                                                                                                                                      |                                                                                                                                                     TGTACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCG‑TAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATC‑ACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGG  >  NC_002947/5951497‑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: ‑