Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,950,783	T→C	8.7%	R500R (CGT→CGC)  N725S (AAC→AGC)	ppx → ppk ←	exopolyphosphatase polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,950,783	0	T	C	8.7%	192.9 / 6.7	69	R500R (CGT→CGC)  N725S (AAC→AGC)	ppx ppk	exopolyphosphatase polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base C (3/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTT  >  NC_002947/5950638‑5950931                                                                                                                                                  |                                                                                                                                                     tGCCAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGa                                                                                                                                                   <  6:208106/149‑1 (MQ=255)   ccAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGAAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGTGTCGGCTTCGTCCTCAGCGTACGTTGAg                                                                                                                                                  <  5:325713/148‑1 (MQ=255)    cAACGGTAGAACTGCACGGCGTGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAACGTACGCTGAGGa                                                                                                                                                <  3:184268/149‑1 (MQ=255)     aaCGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGCCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGAc                                                                                                                                               >  5:266528/1‑149 (MQ=255)         gTAGAACTGCACGGCGGGTACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCTGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCggg                                                                                                                                           <  6:151194/149‑1 (MQ=255)         gTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCggg                                                                                                                                           <  3:59091/149‑1 (MQ=255)             aaCTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGc                                                                                                                                       >  2:52622/1‑149 (MQ=255)               cTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTc                                                                                                                                     <  2:189991/149‑1 (MQ=255)                          ggACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAACGTACGCTGAGGAcga                                                                                                                                             <  7:119091/130‑3 (MQ=255)                          ggACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAACGTACGCTGAGGAcga                                                                                                                                             >  8:119091/1‑128 (MQ=255)                              aacaGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCggg                                                                                                                                           <  4:195973/128‑1 (MQ=255)                              aacaGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCggg                                                                                                                                           >  3:195973/1‑128 (MQ=255)                              aacaGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAgg                                                                                                                      >  5:116351/1‑149 (MQ=255)                              aacaGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAgg                                                                                                                      <  1:217111/149‑1 (MQ=255)                              aacaGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAACGTACGCTGAGGAcgaa                                                                                                                                            <  1:262783/127‑4 (MQ=255)                              aacaGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAACGTACGCTGAGGAcgaa                                                                                                                                            >  2:262783/1‑124 (MQ=255)                                caGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGt                                                                                                                    >  8:70868/1‑149 (MQ=255)                                   ccTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGc                                                                                                                 <  7:70868/149‑1 (MQ=255)                                    cTGGATGTAGCCTTCCCCGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGcc                                                                                                                <  1:352665/149‑1 (MQ=255)                                           tAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGc                                                                                                         <  7:167883/149‑1 (MQ=255)                                               cTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGcc                                                                                                                <  5:386405/138‑1 (MQ=255)                                               cTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGcc                                                                                                                >  6:386405/1‑138 (MQ=255)                                                    cTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGt                                                                                                >  2:77752/1‑149 (MQ=255)                                                    cTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGt                                                                                                >  7:178060/1‑149 (MQ=255)                                                     tGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGtt                                                                                               >  2:56134/1‑149 (MQ=255)                                                       aaGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCt                                                                                             <  1:199965/149‑1 (MQ=255)                                                          gTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGt                                                                                          >  2:67513/1‑149 (MQ=255)                                                              gCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGcc                                                                                      >  5:150953/1‑149 (MQ=255)                                                              gCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGc                                                                                       <  1:52622/148‑1 (MQ=255)                                                                  cagcagAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTc                                                                                  >  5:313455/1‑149 (MQ=255)                                                                   agcagAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGc                                                                                                                 <  6:206710/117‑1 (MQ=255)                                                                   agcagAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGc                                                                                                                 >  5:206710/1‑117 (MQ=255)                                                                    gcagAACCAGCTGACCAAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACAGGGTTGCTCAGCCGCTCCAGCAGGTTCGCCTGAACGCTGCGCGGGTTCTGGTTGCAGCTCgt                                                                                >  2:159016/1‑148 (MQ=255)                                                                     cagAACCAGCTGACGCAGGCGGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCggg                                                                               <  6:266528/149‑1 (MQ=255)                                                                     cagAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCggg                                                                               >  7:44639/1‑149 (MQ=255)                                                                     cagAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCAGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCggg                                                                               <  8:325302/149‑1 (MQ=255)                                                                       gAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGt                                                                             >  7:32386/1‑149 (MQ=255)                                                                          ccAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGtgct                                                                          >  7:343883/1‑149 (MQ=255)                                                                          ccAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGtgct                                                                          <  2:114834/149‑1 (MQ=255)                                                                            aGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGtgctgc                                                                        >  1:300535/1‑149 (MQ=255)                                                                                   ccAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACg                                                                   >  4:264551/1‑147 (MQ=255)                                                                                   ccAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACg                                                                   >  8:227469/1‑147 (MQ=255)                                                                                               cGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTgg                                                     >  3:251577/1‑149 (MQ=255)                                                                                               cGCCAACGAGGCGGAGGGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGCGCTACGTGCGCAGCAccccgac                                                       <  7:339080/147‑5 (MQ=255)                                                                                                  cAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTg                                                  >  7:330090/1‑149 (MQ=255)                                                                                                    aCGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCTCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCa                                                <  5:92799/149‑1 (MQ=255)                                                                                                      gAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAg                                              <  2:300535/149‑1 (MQ=255)                                                                                                         gCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTc                                           <  6:313455/149‑1 (MQ=255)                                                                                                           ggAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCa                                         <  1:67513/149‑1 (MQ=255)                                                                                                                 gctggcCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTg                                   >  5:369176/1‑149 (MQ=255)                                                                                                                 gctggcCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGAGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTCGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTg                                   >  1:264454/1‑149 (MQ=255)                                                                                                                  ctggcCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGaa                                                                                                            <  6:118443/75‑1 (MQ=255)                                                                                                                  ctggcCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGaa                                                                                                            >  5:118443/1‑75 (MQ=255)                                                                                                                  ctggcCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTgg                                  <  8:32386/149‑1 (MQ=255)                                                                                                                   tggcCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTggg                                 >  7:279588/1‑149 (MQ=255)                                                                                                                           gTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCg                         <  8:44639/149‑1 (MQ=255)                                                                                                                            tCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCgg                        >  6:79677/1‑149 (MQ=255)                                                                                                                              ggCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTc                                           >  4:141063/1‑128 (MQ=255)                                                                                                                              ggCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTc                                           <  3:141063/128‑1 (MQ=255)                                                                                                                              cgCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTc                      <  3:133635/148‑1 (MQ=255)                                                                                                                               gCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCa                     <  6:116351/149‑1 (MQ=255)                                                                                                                                   cGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCCGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGt                  <  1:77752/148‑1 (MQ=255)                                                                                                                                    gTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTa                 <  1:56134/148‑1 (MQ=255)                                                                                                                                     tCCTCAACGTACGCTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGc               >  4:74363/1‑149 (MQ=255)                                                                                                                                       cTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTgg                                  <  2:114897/128‑1 (MQ=255)                                                                                                                                       cTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGAGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTgg                                  >  1:114897/1‑128 (MQ=255)                                                                                                                                        tCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCg                         >  1:268676/1‑136 (MQ=255)                                                                                                                                        tCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCg                         <  2:268676/136‑1 (MQ=255)                                                                                                                                             gTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGc       <  8:343883/149‑1 (MQ=255)                                                                                                                                                  ttGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCctt  <  5:79677/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TGCCAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTT  >  NC_002947/5950638‑5950931

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