Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,950,776	G→A	8.2%	S498N (AGC→AAC)  R727R (CGC→CGT)	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,776	0	G	A	8.2%	210.5 / 7.0	73	S498N (AGC→AAC)  R727R (CGC→CGT)	ppx ppk	exopolyphosphatase polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base A (3/3);  total (37/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 = 8.55e-01

CAGCAGATGCCAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGC  >  NC_002947/5950631‑5950918                                                                                                                                                  |                                                                                                                                               cagcagATGCCAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCg                                                                                                                                              >  8:167883/1‑148 (MQ=255) cagcagATGCCAACGGTAGAACTGCACGGCGGGCACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCg                                                                                                                                              >  2:217111/1‑148 (MQ=255)   gcagATGCCAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTAc                                                                                                                                           >  2:352665/1‑149 (MQ=255)        tGCCAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGa                                                                                                                                      <  6:208106/149‑1 (MQ=255)          ccAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGAAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGTGTCGGCTTCGTCCTCAGCGTACGTTGAg                                                                                                                                     <  5:325713/148‑1 (MQ=255)           cAACGGTAGAACTGCACGGCGTGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAACGTACGCTGAGGa                                                                                                                                   <  3:184268/149‑1 (MQ=255)           cAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTTGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCg                                                                                                                                              >  5:265593/1‑138 (MQ=255)           cAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTTGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCg                                                                                                                                              <  6:265593/138‑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                                                                                                         <  1:217111/149‑1 (MQ=255)                                     aacaGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAgg                                                                                                         >  5:116351/1‑149 (MQ=255)                                     aacaGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAACGTACGCTGAGGAcgaa                                                                                                                               >  2:262783/1‑124 (MQ=255)                                     aacaGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAACGTACGCTGAGGAcgaa                                                                                                                               <  1:262783/127‑4 (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                                                      >  8:227469/1‑147 (MQ=255)                                                                                          ccAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACg                                                      >  4:264551/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)                                                                                                                                                  |                                                                                                                                               CAGCAGATGCCAACGGTAGAACTGCACGGCGGGGACAACAGCCTGGATGTAGCCTTCCCTGAAGGTTGGCTGCAGCAGAACCAGCTGACCCAGGCCGACTTCGCCAACGAGGCGGAGTGGCTGGCCCGGGTCGGCTTCGTCCTCAGCGTACGTTGAGGACCGGGTTGCTCAGGCGCTCCAGCAGGGTCGCCTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGC  >  NC_002947/5950631‑5950918

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