Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,320,485	Δ1 bp	5.8%	coding (756/945 nt)	PP_2041 ←	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,320,485	0	G	.	5.8%	220.0 / 6.6	69	coding (756/945 nt)	PP_2041	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base . (2/2);  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 = 2.23e-01

CCTTTCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCT  >  NC_002947/2320339‑2320630                                                                                                                                                   |                                                                                                                                                   ccTTTCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCt                                                                                                                                                  >  8:270475/1‑149 (MQ=255)  cTTTCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTg                                                                                                                                                 >  1:250899/1‑149 (MQ=255)  cTTTCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCCCACCAGGTCCCCCATGCAGGCCGCGGCGAACGACAGCGGCTTTTTCTGCGAGTGGGCCAAGGCGATTGTGCGGGTGATTCGCGGTTCCACTATCAGCCGCGCGCTg                                                                                                                                                 >  1:256382/1‑149 (MQ=255)   tttCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTgg                                                                                                                                                >  6:83110/1‑149 (MQ=255)       cAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAg                                                                                                                                             <  5:366695/148‑1 (MQ=255)        aGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑gc                                                                                                                                          >  3:228747/1‑149 (MQ=255)         gcgcTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑gcg                                                                                                                                         >  6:316525/1‑149 (MQ=255)          cgctcgctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATAAGCCGCGCGCTGGCAG‑gcgc                                                                                                                                        >  6:331999/1‑149 (MQ=255)           gctcgctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑Gcgcg                                                                                                                                       >  4:215405/1‑149 (MQ=255)            ctcgctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGC‑CTGCCAGCGCGCgg                                                                                                                                      <  5:137682/149‑1 (MQ=255)              cgctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGcc                                                                                                                                    <  7:160208/149‑1 (MQ=255)               gctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGcc                                                                                                                                    <  5:316525/148‑1 (MQ=255)                ctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGc                                                                                                                                     >  7:92665/1‑146 (MQ=255)                   cgcACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGC‑CTGCCAGCGCGCGGctgata                                                                                                                                <  4:123724/148‑6 (MQ=255)                    gcACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATcc                                                                                                                              <  8:29521/149‑1 (MQ=255)                    gcACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATcc                                                                                                                              >  6:352570/1‑149 (MQ=255)                     cACCGCATCCACCAGCCGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCca                                                                                                                             <  6:258459/149‑1 (MQ=255)                     cACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCca                                                                                                                             <  8:344232/149‑1 (MQ=255)                      aCCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGAGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCcag                                                                                                                            >  4:330128/1‑149 (MQ=255)                        cGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCcagca                                                                                                                          <  1:115598/149‑1 (MQ=255)                              caccaGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAgcggcg                                                                                                                    <  7:270475/149‑1 (MQ=255)                              caccaGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAgcggc                                                                                                                     <  2:362405/148‑1 (MQ=255)                               cacacccgTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGcggcgg                                                                                                                   <  5:397723/142‑1 (MQ=255)                                    cagTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGt                                                                                                              >  3:279103/1‑149 (MQ=255)                                    cagTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGt                                                                                                              >  3:282519/1‑149 (MQ=255)                                          gcgcACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGa                                                                                                        <  3:148158/149‑1 (MQ=255)                                           cgcACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAT‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGAt                                                                                                       <  2:252217/149‑1 (MQ=255)                                               ccAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGAtggtg                                                                                                   >  6:13152/1‑149 (MQ=255)                                                   gTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGgccg                                                                                               >  3:115540/1‑149 (MQ=255)                                                   gTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGgccg                                                                                               >  1:280941/1‑149 (MQ=255)                                                      gcgcATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCgcca                                                                                            <  7:268561/149‑1 (MQ=255)                                                       cgcATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGGTGGTGGCCgccag                                                                                           <  4:228747/149‑1 (MQ=255)                                                        gcATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCgccagg                                                                                          >  2:209507/1‑149 (MQ=255)                                                         cATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCgccaggc                                                                                         >  2:338993/1‑149 (MQ=255)                                                                 cggcggcgAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACgg                                                                                 <  4:279103/149‑1 (MQ=255)                                                                 cggcggcgAACGACAGCGGCTTGGGCTGCGAGGGGCCCAAGGCGATGGTGCGGGCGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACgg                                                                                 <  4:115540/149‑1 (MQ=255)                                                                   gcggcgAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCTCCAGGCCAGCACGGAc                                                                               <  2:256382/149‑1 (MQ=255)                                                                   gcggcgAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGAc                                                                               <  2:250899/149‑1 (MQ=255)                                                                    cggcgAACGACAGCTGATTGTGCTGCGTGTGGGCCAAGGGTATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCTGCCAGATCCAGCTGCTGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGAct                                                                              <  5:69888/149‑2 (MQ=255)                                                                            gACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGt                                                                      <  8:92665/149‑1 (MQ=255)                                                                             aCAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGt                                                                      >  6:107807/1‑148 (MQ=255)                                                                                     ttGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCggg                                                             >  3:102569/1‑149 (MQ=255)                                                                                        tgctgaGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGc                                                          <  1:264668/149‑1 (MQ=255)                                                                                           tgcGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGa                                                       >  1:241648/1‑149 (MQ=255)                                                                                             cGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGa                                                       <  4:282519/147‑1 (MQ=255)                                                                                             cGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGAct                                                     >  4:287852/1‑149 (MQ=255)                                                                                                gTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGActgct                                                  <  4:102569/149‑1 (MQ=255)                                                                                                gTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGActgct                                                  <  5:83110/149‑1 (MQ=255)                                                                                                    gCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTc                                              <  5:13152/149‑1 (MQ=255)                                                                                                    gCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTc                                              >  4:7280/1‑149 (MQ=255)                                                                                                     ccAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCg                                             >  3:164511/1‑149 (MQ=255)                                                                                                      cAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGt                                            <  1:338993/149‑1 (MQ=255)                                                                                                       aaGGCGATGGTGCGGGGGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGGGATACCGGGGGCCGACTGCTGTTCGTa                                           <  2:266466/149‑1 (MQ=255)                                                                                                        aGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAg                                          >  5:235033/1‑149 (MQ=255)                                                                                                           cGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTg                                       <  3:7280/149‑1 (MQ=255)                                                                                                                 tGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTc                                  >  3:290662/1‑148 (MQ=255)                                                                                                                  gCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCa                                <  7:267810/149‑1 (MQ=255)                                                                                                                    gggTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAgt                              >  4:364734/1‑149 (MQ=255)                                                                                                                     ggTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAgtg                             >  5:81892/1‑149 (MQ=255)                                                                                                                         aTGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTcac                         <  5:352570/149‑1 (MQ=255)                                                                                                                             gcgGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCg                     <  3:287852/149‑1 (MQ=255)                                                                                                                                gTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCa                                >  3:390617/1‑135 (MQ=255)                                                                                                                                gTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCa                                <  4:390617/135‑1 (MQ=255)                                                                                                                                   ccACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTc                >  2:113171/1‑148 (MQ=255)                                                                                                                                    cACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCa              <  6:318745/149‑1 (MQ=255)                                                                                                                                     aCTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCAt             >  1:362374/1‑149 (MQ=255)                                                                                                                                       tATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATgg           <  3:215405/149‑1 (MQ=255)                                                                                                                                          tggCCGCGC‑CTGCCAGCGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGCCTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAg        >  3:91951/3‑149 (MQ=255)                                                                                                                                          tggCCGCGC‑CTGCCAGCGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAg        >  4:10888/3‑149 (MQ=255)                                                                                                                                              cgcgcgCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTg    >  7:260510/1‑149 (MQ=255)                                                                                                                                               gcgcgcTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGc   >  2:316595/1‑149 (MQ=255)                                                                                                                                               gcgcgcTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGc   <  3:330128/149‑1 (MQ=255)                                                                                                                                                cgcgcTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCt  <  2:280941/149‑1 (MQ=255)                                                                                                                                                cgcgcTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCt  <  3:364734/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CCTTTCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAG‑GCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCT  >  NC_002947/2320339‑2320630

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