Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,958,364	A→G	6.2%	M590V (ATG→GTG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,958,364	0	A	G	6.2%	185.3 / 2.3	64	M590V (ATG→GTG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (34/26);  new base G (2/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

CTGGGCTGTCTGCGGCAAAGTATCCCCCTGCCCTTGCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTGCGCTGCCTGCGCGA  >  NC_002947/5958217‑5958512                                                                                                                                                    |                                                                                                                                                      ctgGGCTGTCTGCGGCAAAGTATCCCCCTGCCCTTGCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCtgat                                                                                                                                                      <  5:42504/149‑1 (MQ=255)             cGGCAAAGTATCCCCCTGCCCTTGCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACagag                                                                                                                                            <  7:59514/147‑1 (MQ=255)              ggCAAAGTATCCCCCTGCCCTTGCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAgcg                                                                                                                                          >  7:202253/1‑148 (MQ=255)                    gTATCCCCCTGCCCTTGCCGACCGTCAGCGAAGTGCGCCTTTCGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGAc                                                                                                                                   <  5:36411/149‑1 (MQ=255)                          cccTGCCCTTTCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTg                                                                                                                             <  6:365877/149‑1 (MQ=255)                          cccTGCCCTTGCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGtt                                                                                                                              >  1:204150/1‑148 (MQ=255)                          cccTGCCCTTGCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGtt                                                                                                                              >  1:279276/1‑148 (MQ=255)                            cTGCCCTTGCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAg                                                                                                                           >  7:137530/1‑149 (MQ=255)                               cccTTGCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTa                                                                                                                        <  3:327110/149‑1 (MQ=255)                                    gCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTg                                                                                                                   >  8:175658/1‑149 (MQ=255)                                     ccGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTgg                                                                                                                  >  6:310316/1‑149 (MQ=255)                                     ccGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGGGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTg                                                                                                                   <  3:24185/148‑1 (MQ=255)                                     ccGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGCTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCAGAACGTGCTGATGACAACAGAGCGCACCGAAGCGTTGAGCTACGCTgg                                                                                                                  >  3:56877/1‑149 (MQ=255)                                      cGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGGTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGCGCGCACCGACGCGTTGAGCTACGCAGGc                                                                                                                 >  4:162249/1‑149 (MQ=255)                                          cGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGtt                                                                                                                              <  3:131264/132‑1 (MQ=255)                                          cGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGtt                                                                                                                              >  4:131264/1‑132 (MQ=255)                                              aGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTg                                                                                                                   <  4:45568/139‑1 (MQ=255)                                              aGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTg                                                                                                                   >  3:45568/1‑139 (MQ=255)                                                     tccgcCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGGGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTg                                                                                                  <  8:135567/147‑1 (MQ=255)                                                       cgcCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTgg                                                                                                 <  7:175658/148‑1 (MQ=255)                                                          cTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCt                                                                                             <  6:358381/149‑1 (MQ=255)                                                              tGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGAcc                                                                                         >  1:4198/1‑149 (MQ=255)                                                               gCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGGTGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGAccc                                                                                        >  3:63708/1‑149 (MQ=255)                                                               gCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGGTGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGAccc                                                                                        >  3:59025/1‑149 (MQ=255)                                                                cAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCCCCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCt                                                                                       >  2:264584/1‑149 (MQ=255)                                                                cAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCtgat                                                                                                                                                      >  1:312941/1‑86 (MQ=255)                                                                cAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCtgat                                                                                                                                                      <  2:312941/86‑1 (MQ=255)                                                                      aGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCCCCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGAc                                                                                   >  3:73704/1‑147 (MQ=255)                                                                      aGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGGTGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCa                                                                                 <  4:63708/149‑1 (MQ=255)                                                                      aGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGGTGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCa                                                                                 <  4:59025/149‑1 (MQ=255)                                                                        cGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGt                                                                                                >  7:134093/1‑132 (MQ=255)                                                                        cGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGt                                                                                                <  8:134093/132‑1 (MQ=255)                                                                               gATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGGTGACCGTGGAc                                                                                   <  6:353497/138‑1 (MQ=255)                                                                               gATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGAc                                                                                   >  5:353497/1‑138 (MQ=255)                                                                               gATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGt                                                                        <  4:73704/149‑1 (MQ=255)                                                                                 tGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTc                                                                       >  3:227732/1‑148 (MQ=255)                                                                                   aGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCaa                                                                     <  8:202253/148‑1 (MQ=255)                                                                                      tgttgtTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGc                                                                 >  3:16092/1‑149 (MQ=255)                                                                                       gttgttGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCt                                                                >  4:292282/1‑149 (MQ=255)                                                                                           ttgCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGaa                                                            <  5:309343/149‑1 (MQ=255)                                                                                             gCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAc                                                           >  5:267725/1‑148 (MQ=255)                                                                                              cTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGAGTTGCGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCACAAGCTGCAAC‑Ga                                                        >  6:10393/1‑149 (MQ=255)                                                                                                 gTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAg                                                       <  3:162249/147‑1 (MQ=255)                                                                                                     aCGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTg                                                 <  5:10393/149‑1 (MQ=255)                                                                                                     aCGTCCGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTg                                                 <  4:227732/149‑1 (MQ=255)                                                                                                       gTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTgg                                                <  6:258778/148‑1 (MQ=255)                                                                                                       gTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTgg                                                >  7:80526/1‑148 (MQ=255)                                                                                                       gTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTgg                                                <  8:80526/148‑1 (MQ=255)                                                                                                       gTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTgg                                                >  5:258778/1‑148 (MQ=255)                                                                                                        tCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTc                                              <  4:16092/149‑1 (MQ=255)                                                                                                         cGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTc                                              >  4:122672/1‑148 (MQ=255)                                                                                                                gACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTa                                      >  3:119822/1‑149 (MQ=255)                                                                                                                gACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTa                                      >  4:38158/1‑149 (MQ=255)                                                                                                                gACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTa                                      >  4:129948/1‑149 (MQ=255)                                                                                                                     gCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACg                                 >  2:287083/1‑149 (MQ=255)                                                                                                                     gCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACg                                 >  4:111610/1‑149 (MQ=255)                                                                                                                            caccacCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAaca                          >  6:118589/1‑149 (MQ=255)                                                                                                                            caccacCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGT‑CAGCGCTATGACGGCGAacac                         >  5:323897/1‑149 (MQ=255)                                                                                                                              ccaccGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACAcg                        <  6:267725/149‑1 (MQ=255)                                                                                                                                  cgcgACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCt                    <  3:122672/149‑1 (MQ=255)                                                                                                                                       ccTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTgc               >  5:351215/1‑149 (MQ=255)                                                                                                                                        cTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTgcg              >  5:298113/1‑149 (MQ=255)                                                                                                                                          gAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTGCGct            >  1:62933/1‑149 (MQ=255)                                                                                                                                           aaCGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTgcgc             <  3:292282/147‑1 (MQ=255)                                                                                                                                            aCGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTGCGctgc          <  5:310316/149‑1 (MQ=255)                                                                                                                                             cGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAACAGAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTGCGctgc          <  1:264584/149‑1 (MQ=255)                                                                                                                                                gCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTGCGctgcctgc      >  8:216378/1‑149 (MQ=255)                                                                                                                                                   gatgaCCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACAcgcg                      >  3:311149/1‑130 (MQ=255)                                                                                                                                                   gatgaCCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACAcgcg                      <  4:311149/130‑1 (MQ=255)                                                                                                                                                    atgaCCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGCCCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTGCGCTGCCTGCGCGa  >  8:100332/1‑149 (MQ=255)                                                                                                                                                    atgaCCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAgg                                                      <  7:304907/97‑1 (MQ=255)                                                                                                                                                    atgaCCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAgg                                                      >  8:304907/1‑97 (MQ=255)                                                                                                                                                    atgaCCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTgc               >  3:1346/1‑136 (MQ=255)                                                                                                                                                    atgaCCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTgc               <  4:1346/136‑1 (MQ=255)                                                                                                                                                    atgaCCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTGCGCTGCCTgcgc    <  3:38158/147‑1 (MQ=255)                                                                                                                                                    atgaCCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTGCGCTGCCTgcgc    <  2:204150/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CTGGGCTGTCTGCGGCAAAGTATCCCCCTGCCCTTGCCGACCGTCAGCGAAGTGCGCCTTTTGCAGCCCAGCGTTGCCGATGAGGTGTTGTTGCTGGTCAACGTCGGCGTCGACCCGCTGCGCCACCACCGCGACCTGAACGTGCTGATGACCACAGAGCGCACCGACGCGTTGAGCTACGCTGGCGTGCGCGAAAACCTGGTGCTGACCCTGGACCAGGTCACCGTCAACAGCTGGAAC‑GAGGTGCTGGTCCAGCGCTATGACGGCGAACACGCGCTGGTGCGCTGCCTGCGCGA  >  NC_002947/5958217‑5958512

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