Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,951,594	G→A	13.3%	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,594	0	G	A	13.3%	171.6 / 14.3	69	L455L (CTG→TTG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base A (4/5);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

GGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGCG  >  NC_002947/5951451‑5951738                                                                                                                                                 |                                                                                                                                                 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                                                                                     >  2:397313/1‑136 (MQ=255)                                                                      cAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCa                                                                                     <  1:397313/136‑1 (MQ=255)                                                                        gAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCg                                                                      >  3:354373/1‑149 (MQ=255)                                                                                  catcTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGa                                                                                             <  2:228791/116‑1 (MQ=255)                                                                                  catcTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGa                                                                                             >  1:228791/1‑116 (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)                                                                                                                                          gcctGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTg     <  2:56590/146‑1 (MQ=255)                                                                                                                                            ttGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGgc   >  3:290793/1‑149 (MQ=255)                                                                                                                                             tGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGgcg  >  7:24293/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                 GGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCA‑ACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGCG  >  NC_002947/5951451‑5951738

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