Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,882,349	A→G	5.0%	I304V (ATC→GTC) ‡	serA →	bifunctional D‑3‑phosphoglycerate dehydrogenase/alpha‑ketoglutarate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,882,349	0	A	G	5.0%	296.3 / 3.7	99	I304V (ATC→GTC) ‡	serA	bifunctional D‑3‑phosphoglycerate dehydrogenase/alpha‑ketoglutarate reductase
Reads supporting (aligned to +/- strand):  ref base A (42/52);  new base G (3/2);  total (45/54)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GGACAAGCACCTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGAAAACATCCCGGGCGTGCTCA  >  NC_002947/5882201‑5882496                                                                                                                                                     |                                                                                                                                                     ggACAAGCACCTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACa                                                                                                                                                      <  4:241584/149‑1 (MQ=255)   aCAAGCACCTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATc                                                                                                                                                    <  1:134165/149‑1 (MQ=255)      aGCACCTGATCGGCGCCACCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGc                                                                                                                                                 <  2:55767/149‑1 (MQ=255)       gCACCTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCg                                                                                                                                                   <  5:353855/146‑1 (MQ=255)       gCACCTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCg                                                                                                                                                   >  6:353855/1‑146 (MQ=255)        cACCTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCt                                                                                                                                               >  6:443360/1‑149 (MQ=255)        cACCTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCt                                                                                                                                               <  4:389442/149‑1 (MQ=255)         aCCTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCGATGTTGGCCTg                                                                                                                                              <  2:43865/149‑1 (MQ=255)          ccTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTgg                                                                                                                                             <  4:186723/149‑1 (MQ=255)          ccTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTg                                                                                                                                              >  3:58246/1‑148 (MQ=255)           cTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTgg                                                                                                                                             >  8:100887/1‑148 (MQ=255)             gATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAg                                                                                                                                          >  5:124489/1‑149 (MQ=255)                cGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTgg                                                                                                                                       >  6:435119/1‑149 (MQ=255)                  gcgcCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGcc                                                                                                                                     <  8:303434/149‑1 (MQ=255)                  gcgcCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGcc                                                                                                                                     >  4:304100/1‑149 (MQ=255)                  gcgcCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACAACGGCCTGGAAGTGGcc                                                                                                                                     >  5:340312/1‑149 (MQ=255)                    gccgccATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGcc                                                                                                                                     >  2:186203/1‑147 (MQ=255)                    gccgccATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCg‑a                                                                                                                                  <  6:164867/149‑1 (MQ=255)                    gccgccATCGACGTGTCCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCg‑a                                                                                                                                  <  3:340270/149‑1 (MQ=255)                      cgccATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCg‑ag                                                                                                                                 >  8:364641/1‑148 (MQ=255)                        ccATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGaaa                                                                                                                              <  5:129846/149‑1 (MQ=255)                          aTCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACt                                                                                                                            >  2:182695/1‑149 (MQ=255)                           tCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTg                                                                                                                           <  7:100887/149‑1 (MQ=255)                             gACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTgg                                                                                                                          <  2:310112/148‑1 (MQ=255)                             gACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGt                                                                                                                         <  1:111829/149‑1 (MQ=255)                             gACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGt                                                                                                                         >  3:286786/1‑149 (MQ=255)                              aCGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTc                                                                                                                        <  5:435119/149‑1 (MQ=255)                              aCGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTc                                                                                                                        <  1:402410/149‑1 (MQ=255)                              aCGTGGTCCCGGTCGAGCCACGCTCGAACGACGTAGAGTGCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGTTGGTTCCACCGCCGAAGCCCAGGCCAACATCTGCCTGGAAGTGGCCG‑AGAAACTgggc                                                                                                                        <  6:340312/149‑3 (MQ=255)                                gtgtTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCaa                                                                                                                      <  3:109758/149‑1 (MQ=255)                                gtgtTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCaa                                                                                                                      <  8:126642/149‑1 (MQ=255)                                 tgtTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAg                                                                                                                     <  8:134201/149‑1 (MQ=255)                                 tgtTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAg                                                                                                                     >  2:258729/1‑149 (MQ=255)                                 tgtTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAg                                                                                                                     <  8:146081/149‑1 (MQ=255)                                   ttCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTgggggagc                                                                                                                    >  1:151350/1‑142 (MQ=255)                                   ttCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTa                                                                                                                   >  7:6198/1‑149 (MQ=255)                                   ttCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTa                                                                                                                   <  4:106774/149‑1 (MQ=255)                                   ttCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTa                                                                                                                   <  4:148103/149‑1 (MQ=255)                                     cccGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTc                                                                                                                        <  3:100146/142‑1 (MQ=255)                                     cccGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTc                                                                                                                        >  4:100146/1‑142 (MQ=255)                                              gCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACg                                                                                                        <  1:182695/149‑1 (MQ=255)                                               ccACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACgg                                                                                                       >  5:170101/1‑149 (MQ=255)                                               ccACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACgg                                                                                                       <  4:58246/149‑1 (MQ=255)                                                cACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACACCgg                                                                                                       >  6:113454/1‑148 (MQ=255)                                                 aCGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTa                                                                                                     <  7:364641/149‑1 (MQ=255)                                                       cAACGACGTAGAGTTCGAGAGGCGGCTGCGCGGCCTCGACTACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGGCCACGCCAACATCGTCCTGGAAGTGGCCG‑AGATACTGGTCATGTACATCGAGAACGTTACCTCg                                                                                                <  2:276963/148‑1 (MQ=255)                                                       cAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCg                                                                                                <  3:58135/148‑1 (MQ=255)                                                        aaCGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGgt                                                                                              <  8:75325/149‑1 (MQ=255)                                                         acgacgAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGgtg                                                                                             >  6:20432/1‑149 (MQ=255)                                                              gaAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATGCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCAAACATCGGCCTGGAAGTGGCCGCAGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCt                                                                                         <  1:84526/149‑1 (MQ=255)                                                               aaGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCg                                                                                       >  7:97899/1‑149 (MQ=255)                                                               aaGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCg                                                                                       >  1:362516/1‑149 (MQ=255)                                                                agagTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACACCGGTACCTCGGTGTCCTggg                                                                                      >  4:51628/1‑149 (MQ=255)                                                                agagTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCgg                                                                                      <  1:186203/149‑1 (MQ=255)                                                                agagTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCgg                                                                                      <  1:258729/149‑1 (MQ=255)                                                                 gagTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGt                                                                                     <  6:170101/149‑1 (MQ=255)                                                                 gagTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGt                                                                                     <  3:13766/149‑1 (MQ=255)                                                                  agTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCTTCGGTc                                                                                    <  3:19501/149‑1 (MQ=255)                                                                    ttCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCaa                                                                                  <  7:43602/149‑1 (MQ=255)                                                                      cGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAATGGTACCTCGGTGTCCTCGGTCAAc                                                                                 <  4:49904/148‑1 (MQ=255)                                                                       gagaGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGGACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACtt                                                                               >  5:176628/1‑149 (MQ=255)                                                                         gagCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTcc                                                                             <  3:304100/149‑1 (MQ=255)                                                                          agCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTccc                                                                            <  7:388156/149‑1 (MQ=255)                                                                          agCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTccc                                                                            <  6:124489/149‑1 (MQ=255)                                                                          agCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTccc                                                                            <  8:97899/149‑1 (MQ=255)                                                                            cccGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCgg                                                                          <  4:352062/149‑1 (MQ=255)                                                                                   cgcgGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGc                                                                   <  8:6198/149‑1 (MQ=255)                                                                                        ccTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGc                                                              <  6:131968/149‑1 (MQ=255)                                                                                         cTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGc                                                              >  4:55427/1‑148 (MQ=255)                                                                                            gACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAgc                                                          >  7:363756/1‑149 (MQ=255)                                                                                              cAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAgcg                                                         >  1:21012/1‑148 (MQ=255)                                                                                               aaCGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAgc                                                          >  3:34473/1‑146 (MQ=255)                                                                                               aaCGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAgc                                                          <  4:34473/146‑1 (MQ=255)                                                                                               aaCGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCa                                                       >  3:67087/1‑149 (MQ=255)                                                                                                        cTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGACCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCa                                               >  2:9536/1‑148 (MQ=255)                                                                                                          gACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCTGGCaaa                                             >  1:170858/1‑148 (MQ=255)                                                                                                               cgcACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGt                                       <  5:443360/149‑1 (MQ=255)                                                                                                                    aTTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTctgct                                  >  5:293009/1‑149 (MQ=255)                                                                                                                    aTTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTctgct                                  >  2:430128/1‑149 (MQ=255)                                                                                                                      tggtggTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTctgctg                                 >  1:427556/1‑148 (MQ=255)                                                                                                                        gtggtTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGca                               >  4:8077/1‑148 (MQ=255)                                                                                                                         tggtTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGcaca                             >  1:165453/1‑149 (MQ=255)                                                                                                                         tggtTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGcaca                             <  3:55427/149‑1 (MQ=255)                                                                                                                             tCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAg                                                  >  4:388686/1‑124 (MQ=255)                                                                                                                             tCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAg                                                  <  3:388686/124‑1 (MQ=255)                                                                                                                               cACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCAc                        >  8:189705/1‑148 (MQ=255)                                                                                                                               cACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACg                       <  4:67087/149‑1 (MQ=255)                                                                                                                                  cgccgAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGaa                     <  2:436945/148‑1 (MQ=255)                                                                                                                                  cgccgAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGaa                     >  1:436945/1‑148 (MQ=255)                                                                                                                                     ccacttccAGGCCGATGTTGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTctgctg                                 >  6:264395/7‑133 (MQ=255)                                                                                                                                     ccacttccAGGCCGATGTTGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTctgctg                                 <  5:264395/127‑1 (MQ=255)                                                                                                                                         ttccAGGCCGATGTTGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAACCACCGTCTGCTGCACATTCACGAAAACATccc             >  5:177379/3‑149 (MQ=255)                                                                                                                                         ttccAGGCCGATGTTGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGAAAACATccc             >  5:184159/3‑149 (MQ=255)                                                                                                                                          cccAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGAAAACATCCCg            >  7:301218/1‑149 (MQ=255)                                                                                                                                            cAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGAAAACATCCCggg          <  5:20432/149‑1 (MQ=255)                                                                                                                                            cAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGAAAACATCCCggg          >  7:156153/1‑149 (MQ=255)                                                                                                                                              ggCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGAAAACATCCCGGGCg        <  2:427556/149‑1 (MQ=255)                                                                                                                                               gCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGAAAACATCCCGGGCGt       >  8:39333/1‑149 (MQ=255)                                                                                                                                                 cAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCACACACCGTCTGCTGCACATTCACGAAAACATCCCGGGCGTGc     >  7:34283/1‑149 (MQ=255)                                                                                                                                                    cATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGAAAACATCCCGGGCGTGCTCa  <  7:39333/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGACAAGCACCTGATCGGCGCCGCCATCGACGTGTTCCCGGTCGAGCCACGCTCCAACGACGAAGAGTTCGAGAGCCCGCTGCGCGGCCTCGACAACGTGATCCTGACGCCGCACATTGGTGGTTCCACCGCCGAAGCCCAGGCCAACATCGGCCTGGAAGTGGCCG‑AGAAACTGGTCAAGTACAGCGACAACGGTACCTCGGTGTCCTCGGTCAACTTCCCGGAAGTGGCCCTGCCAGCGCACCCAGGCAAACACCGTCTGCTGCACATTCACGAAAACATCCCGGGCGTGCTCA  >  NC_002947/5882201‑5882496

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