Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,948,173	C→G	12.9%	R494R (CGG→CGC)	PP_1749 ←	acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,948,173	0	C	G	12.9%	170.7 / 15.9	62	R494R (CGG→CGC)	PP_1749	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/28);  new base G (5/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGATCAAC  >  NC_002947/1948029‑1948320                                                                                                                                                   |                                                                                                                                                    cccGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGcgccg                                                                                                                                                   <  4:97706/147‑1 (MQ=255)   cGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACg                                                                                                                                               >  8:188720/1‑149 (MQ=255)      cATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGGCGCACGGcc                                                                                                                                            <  4:540/149‑1 (MQ=255)            gCACCATGAAGTCCAGCCCCACCACCGGAAATTCCAGTGCC‑AGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGAGCGCCGCACGGCCAGgcgc                                                                                                                                      >  7:158497/1‑149 (MQ=255)             cACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCCGGCCGGCCAGC‑GCCTGGCCGTGCGGCGCACGGCCaacctgc                                                                                                                                     >  8:193704/1‑143 (MQ=37)             cACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGCTGGCCGTGCGGCGCACGGCCaacc                                                                                                                                        >  7:4756/1‑143 (MQ=255)             cACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGCTGGCCGTGCGGCGCACGGCCaacc                                                                                                                                        <  8:4756/146‑4 (MQ=255)               ccATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTgg                                                                                                                                   <  1:121846/149‑1 (MQ=255)                cATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCCACGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTgg                                                                                                                                   >  1:221157/1‑149 (MQ=255)                  tGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCTCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTggccg                                                                                                                                <  1:119342/149‑1 (MQ=255)                         cAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCcggccggc                                                                                                                          >  7:85045/1‑148 (MQ=255)                           gCCCCACCACCGGAATTTCCAGTGCC‑GGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCTGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAac                                                                                                                       <  8:137247/149‑1 (MQ=255)                                accaccGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcgtc                                                                                                                  <  7:280818/149‑1 (MQ=255)                                  caccGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGgccggccg                                                                                                                            <  6:194309/137‑1 (MQ=255)                                  caccGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGgccggccg                                                                                                                            >  5:194309/1‑137 (MQ=255)                                  caccGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcgtcgt                                                                                                                <  2:221157/149‑1 (MQ=255)                                   accGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGCGCCTggccgtgc                                                                                                                             >  6:44798/1‑132 (MQ=255)                                   accGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGCGCCTggccgtgc                                                                                                                             <  5:44798/135‑4 (MQ=255)                                   accGGAATTTCCACTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGCTGACCTCCTCCAGGGTGCC‑GCCGGTGTGCAGCTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGCCAACAcgtcct                                                                                                                >  8:253361/1‑146 (MQ=255)                                          tttCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGcc                                                                                                         >  6:133771/1‑148 (MQ=255)                                             ccAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAAc                                                                                                      <  7:64506/148‑1 (MQ=255)                                               aGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGcc                                                                                                   >  7:291724/1‑149 (MQ=255)                                                 tGCC‑CGTGCCGCGCGCACGGCAGCATCGGTCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcg                                                                                                                    <  8:80772/130‑1 (MQ=255)                                                 tGCC‑CGTGCCGCGCGCACGGCAGCATCGGTCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACAcg                                                                                                                    >  7:80772/1‑130 (MQ=255)                                                      cGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTgc                                                                                             >  1:56890/1‑149 (MQ=255)                                                       ctGCCGTGCGCGCGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTgcg                                                                                            >  4:35137/2‑149 (MQ=255)                                                          ccGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTgcgcg                                                                                          >  4:187224/1‑148 (MQ=255)                                                             cgcgcACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAg                                                                                      >  8:137587/1‑149 (MQ=255)                                                               cgcACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGt                                                                                    <  2:56890/149‑1 (MQ=255)                                                                      cagcaTCGGCGAGGAGGGGGTGCAGGCGTTCGGTGACGTCCTCCAGTGTGCC‑GCCGGTGTTCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAAGAGGTCGTCATAGCCATGGCCGGCTGCGCGCAAGGTGCGCTCg                                                                             <  8:158497/149‑1 (MQ=255)                                                                      cagcaTCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCg                                                                             <  6:34341/149‑1 (MQ=255)                                                                             ggCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtcg                                                                       <  3:190226/148‑1 (MQ=255)                                                                              gCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGAAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtcgtc                                                                     >  7:260168/1‑149 (MQ=255)                                                                               ccAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtcgtcg                                                                    >  1:218975/1‑149 (MQ=255)                                                                                cAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTcgtcgt                                                                   <  2:97080/149‑1 (MQ=255)                                                                                           gCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATg                                                        >  6:22543/1‑149 (MQ=255)                                                                                              ggCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCgg                                                     <  8:291724/149‑1 (MQ=255)                                                                                              ggCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCgg                                                     <  7:253361/149‑1 (MQ=255)                                                                                               gCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGt                                                                           <  1:305095/126‑1 (MQ=255)                                                                                               gCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGt                                                                           >  2:305095/1‑126 (MQ=255)                                                                                               gCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGc                                                    >  1:164124/1‑149 (MQ=255)                                                                                                cGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGct                                                   <  2:218975/149‑1 (MQ=255)                                                                                                     ggTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGc                                              <  2:164124/149‑1 (MQ=255)                                                                                                              ctccAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTgg                                     <  5:154452/149‑1 (MQ=255)                                                                                                              ctccAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTg                                      <  7:16286/148‑1 (MQ=255)                                                                                                              ctccAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTg                                      >  8:16286/1‑148 (MQ=255)                                                                                                               tccAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtcgtc                                                                     <  6:62748/116‑1 (MQ=255)                                                                                                               tccAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTtcgtc                                                                     >  5:62748/1‑116 (MQ=255)                                                                                                                     gTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTg                              <  6:30175/149‑1 (MQ=255)                                                                                                                                tgtgCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACgg                    <  5:22543/149‑1 (MQ=255)                                                                                                                                    cAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGctct                >  8:182372/1‑149 (MQ=255)                                                                                                                                    cAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGctct                >  3:247801/1‑149 (MQ=255)                                                                                                                                     aGGTTGGCCGTGCGCCGCACTGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTg               <  8:85045/149‑1 (MQ=255)                                                                                                                                     aGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTg               >  7:273140/1‑149 (MQ=255)                                                                                                                                      ggTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACTGGTGGCAGCCTGGCGTCGACGGCTCTgc              <  6:52963/149‑1 (MQ=255)                                                                                                                                      gcctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCgg                                        >  2:108168/4‑123 (MQ=255)                                                                                                                                      gcctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCgg                                        <  1:108168/120‑1 (MQ=255)                                                                                                                                       gTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTgcg             >  4:110714/1‑149 (MQ=255)                                                                                                                                       gTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACCTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTgcg             >  1:228626/1‑149 (MQ=255)                                                                                                                                        ctGGCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTgcgc            >  8:52206/2‑149 (MQ=255)                                                                                                                                          ggCCGTGCGGCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCtt          >  4:174681/1‑149 (MQ=255)                                                                                                                                              gTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCTTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCTAt      <  3:187224/149‑1 (MQ=255)                                                                                                                                              gTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGAt      <  5:58255/149‑1 (MQ=255)                                                                                                                                                  gccgcACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGATCAAc  >  1:146693/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCCGGCATCACGCACCATGAAGTCCAGCCCCACCACCGGAATTTCCAGTGCC‑CGTGCCGCGCGCACGGCAGCATCGGCCAGCACCGGGTGCAGGCGTTCGGTGACGTCCTCCAGGGTGCC‑GCCGGTGTGCAGGTTGGCCGTGCGCCGCACGGCCAGGCGCTGGCCGGCCGGCAACACGTCGTCATAGCCAACGCCCGCTGCGCGCAAGGTGCGCTCGGTTTCGTCGTCCAGCGGAATGCGGCTCTCGCCACCGGTGGCAGCCTGGCGTCGACGGCTCTGCGCTTCGATCAAC  >  NC_002947/1948029‑1948320

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