Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,951,589	T→C	14.3%	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,589	0	T	C	14.3%	165.0 / 19.0	70	Q456Q (CAA→CAG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base C (5/5);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

AGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGC  >  NC_002947/5951446‑5951737                                                                                                                                                 |                                                                                                                                                     aGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCCTGc                                                                                                                                                   >  7:52211/1‑147 (MQ=255)      ggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAagcggc                                                                                                                                            <  5:210677/149‑1 (MQ=255)       gcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCt                                                                                                                                           >  2:136927/1‑149 (MQ=255)        cgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTg                                                                                                                                          <  5:259683/149‑1 (MQ=255)         gTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTgg                                                                                                                                         <  4:176730/149‑1 (MQ=255)         gTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTgg                                                                                                                                         <  8:52211/149‑1 (MQ=255)            gCCGGCGTTGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGGTCGCGGCGCAGAATCTGCTTCATCTTGGCGTGTGTCTTGAAGCCA‑AGCACAGCGTAGATCACTACCTCACGAGCCGCTTGCTGGCGGCTGGcc                                                                                                                                       <  8:173857/148‑1 (MQ=255)                     gTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGt                                                                                                                             >  2:201098/1‑149 (MQ=255)                      tAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGtt                                                                                                                            <  1:201098/149‑1 (MQ=255)                         ttGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGa                                                                                                                         <  2:126429/149‑1 (MQ=255)                         ttGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGa                                                                                                                         >  7:284193/1‑149 (MQ=255)                                  cccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGt                                                                                                                 <  7:55800/148‑1 (MQ=255)                                   ccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGt                                                                                                                 >  8:231396/1‑147 (MQ=255)                                   ccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGt                                                                                                                 <  7:231396/147‑1 (MQ=255)                                     aGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCaa                                                                                                              >  6:335079/1‑148 (MQ=255)                                       aTGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAAc                                                                                                           <  3:168092/149‑1 (MQ=255)                                           gcgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAAcgcg                                                                                                        >  5:56076/1‑148 (MQ=255)                                              tAACGCACGATCTCGCCCTGCCCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGcgcacg                                                                                                    >  8:7946/1‑149 (MQ=255)                                                aCGCACGATCTCGACCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGGCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGcgcacgc                                                                                                   <  7:387413/148‑1 (MQ=255)                                                   cACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggcggc                                                                                                                                            <  1:368760/104‑1 (MQ=255)                                                   cACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggcggc                                                                                                                                            >  2:368760/1‑104 (MQ=255)                                                   cACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACt                                                                                               <  8:284193/149‑1 (MQ=255)                                                          tcGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCAc                                                                                        >  8:90356/1‑149 (MQ=255)                                                            gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCAc                                                                                        <  6:56076/147‑1 (MQ=255)                                                            gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                      >  8:83062/1‑149 (MQ=255)                                                            gCCCTGCCCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCAc                                                                                        <  7:7946/147‑1 (MQ=255)                                                             cccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                      >  1:266661/1‑148 (MQ=255)                                                             cccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCAGGCCATCTGCAGGTTGTACTCTTCGGCAAAACGCGAACGCAACt                                                                                               >  3:55364/1‑139 (MQ=255)                                                             cagtGCTCGCGGCGCAGAATCAGATTCATGTTGGCGTGGGTCTTGAAGGCA‑ACGAGACCGTATATCACTACCGGACGAGCGGCTTGATGGCGGCTGGAGATCTTAAGGTTGGACTCTTGGTGAAAAGGCGCACGCAACt                                                                                               <  4:55364/136‑1 (MQ=255)                                                              ccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                      >  4:35309/1‑147 (MQ=255)                                                                   tCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTg                                                                                                                                          <  3:331884/90‑1 (MQ=255)                                                                   tCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTg                                                                                                                                          >  4:331884/1‑90 (MQ=255)                                                                     gcggcgCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcct                                                                             <  5:335079/149‑1 (MQ=255)                                                                       ggcgCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCACTTTACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcctc                                                                            <  1:136927/149‑1 (MQ=255)                                                                           cAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCa                                                                                    <  1:397313/136‑1 (MQ=255)                                                                           cAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCa                                                                                    >  2:397313/1‑136 (MQ=255)                                                                             gAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCg                                                                     >  3:354373/1‑149 (MQ=255)                                                                                       catcTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGa                                                                                            >  1:228791/1‑116 (MQ=255)                                                                                       catcTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGa                                                                                            <  2:228791/116‑1 (MQ=255)                                                                                       catcTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGcgccg                                                           >  4:407661/1‑149 (MQ=255)                                                                                        atcTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGcgccgc                                                          <  7:83062/149‑1 (MQ=255)                                                                                        atcTTGGCGTGGGTCTTGAAGCCA‑ACCACACCCTAGAACACAACCGCACCAGCCCCTTGCACCCGGCTGGCCATAAGCACGTTGGACTCATCGTCAAAACCCGCACGCAACTCGAACACCGCAGTGACCACCTTGCCGTTACCcgccc                                                           >  6:171489/1‑147 (MQ=255)                                                                                         tcTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGtt                                                                                                                            >  4:293787/1‑82 (MQ=255)                                                                                         tcTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGtt                                                                                                                            <  3:293787/82‑1 (MQ=255)                                                                                               cGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATccacc                                                   <  3:35309/149‑1 (MQ=255)                                                                                                     tCTTGAAGCCA‑ACCACACCCTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCCCGCAACTCGATCACCGCAGTGACCTCCTGGCCGTTACGCGCCGCATCCACCCGGGCg                                             >  2:348124/1‑149 (MQ=255)                                                                                                       ttGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTc                                           >  6:94503/1‑149 (MQ=255)                                                                                                        tGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTc                                           >  1:242615/1‑148 (MQ=255)                                                                                                                  ccACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTcc                                 <  2:242615/149‑1 (MQ=255)                                                                                                                    acacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCg                                                                     >  1:288401/1‑111 (MQ=255)                                                                                                                    acacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCg                                                                     <  2:288401/111‑1 (MQ=255)                                                                                                                     cacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATc                                                       >  4:225395/1‑124 (MQ=255)                                                                                                                     cacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATc                                                       <  3:225395/124‑1 (MQ=255)                                                                                                                     cacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTc                                  <  8:292587/145‑1 (MQ=255)                                                                                                                     cacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTc                                  >  7:292587/1‑145 (MQ=255)                                                                                                                        cGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGCCCTCCTTGCCGTTACGAGCCGCATCCACCAGGGCGTCGCCGATTTCCGAGt                             >  4:111329/1‑147 (MQ=255)                                                                                                                        cGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGt                             >  5:182872/1‑147 (MQ=255)                                                                                                                          tAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGc                         <  5:171489/149‑1 (MQ=255)                                                                                                                            gATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGcc                        <  4:354373/148‑1 (MQ=255)                                                                                                                                    ccGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCgg               >  4:156704/1‑149 (MQ=255)                                                                                                                                    ccGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCgg               >  3:197309/1‑149 (MQ=255)                                                                                                                                       cACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACACTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAg                          <  5:167626/135‑1 (MQ=255)                                                                                                                                       cACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACACTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAg                          >  6:167626/1‑135 (MQ=255)                                                                                                                                         ccAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTAccg          >  4:331978/1‑147 (MQ=255)                                                                                                                                         ccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAg                          >  3:71777/1‑133 (MQ=255)                                                                                                                                         ccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAg                          <  4:71777/133‑1 (MQ=255)                                                                                                                                         ccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAg          >  1:56590/1‑149 (MQ=255)                                                                                                                                          cAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg         >  2:258465/1‑149 (MQ=255)                                                                                                                                            cgcGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGt       <  3:376595/147‑1 (MQ=255)                                                                                                                                              cGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCt     <  7:90356/149‑1 (MQ=255)                                                                                                                                                 ttGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGgc  >  3:290793/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     AGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGC  >  NC_002947/5951446‑5951737

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