Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	188,616	C→T	5.7%	T242I (ACC→ATC) ‡	PP_0165 →	GGDEF domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	188,616	0	C	T	5.7%	263.9 / 3.6	88	T242I (ACC→ATC) ‡	PP_0165	GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (42/41);  new base T (2/3);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

CTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAG  >  NC_002947/188468‑188764                                                                                                                                                     |                                                                                                                                                             cTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAc                                                                                                                                                              >  6:173415/1‑149 (MQ=255) cTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAc                                                                                                                                                              >  1:291747/1‑149 (MQ=255)  tCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGAcc                                                                                                                                                             >  2:1504/1‑149 (MQ=255)       cTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCt                                                                                                                                                        <  3:13887/149‑1 (MQ=255)       cTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCt                                                                                                                                                        >  1:149929/1‑149 (MQ=255)        tGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAAGTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTg                                                                                                                                                       >  3:245143/1‑149 (MQ=255)          ccAGAGTTGCCACGCACGCCGGAGCTACGCCTCTTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGt                                                                                                                                                      <  2:253286/148‑1 (MQ=255)          ccAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGt                                                                                                                                                      <  6:15211/148‑1 (MQ=255)          ccAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTc                                                                                                                                                     <  3:71167/149‑1 (MQ=255)           cAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGTAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcc                                                                                                                                                    <  2:329189/149‑2 (MQ=255)           cAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCCGGTCAGGCTGTcc                                                                                                                                                    <  3:202032/149‑2 (MQ=255)              agTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAc                                                                                                                                               <  2:291747/149‑1 (MQ=255)                  gCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAc‑c‑gc‑‑c                                                                                                                                       >  5:379513/1‑149 (MQ=255)                       gcacgcCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTAt                                                                                                                                  >  2:147178/1‑149 (MQ=255)                       gcacgcCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTCAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTAt                                                                                                                                  >  5:370237/1‑149 (MQ=255)                             cGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGaa                                                                                                                          <  6:370237/149‑1 (MQ=255)                                    aCGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGc                                                                                                                   >  5:296871/1‑149 (MQ=255)                                     cgccgcGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAg                                                                                                                    <  7:6274/147‑1 (MQ=255)                                       ccgcGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTg                                                                                                                 >  8:74606/1‑148 (MQ=255)                                        cgcgTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGaa                                                                                                               >  7:143767/1‑149 (MQ=255)                                        cgcgTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGaa                                                                                                               >  1:329651/1‑149 (MQ=255)                                             aGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCGAGTT‑‑TGAAATGCAGCGGAACAATc                                                                                                          >  8:162351/1‑149 (MQ=255)                                              gTGCAAGCTATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCTGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATc                                                                                                          <  7:68337/148‑1 (MQ=255)                                                gCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGaa                                                                                                               >  1:321717/1‑141 (MQ=255)                                                gCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGaa                                                                                                               <  2:321717/141‑1 (MQ=255)                                                 cAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCggg                                                                                                       >  3:111699/1‑148 (MQ=255)                                                  aaccGATGAACAAGATGGTCGAGAAACTCAAAGCCCTGTTAACCGAGCAGGCAGACCGCAGTGAGAGGCTGCGTGCGGAGTCCTACAAGGACAGCATGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGt                                                                                                      >  2:339483/4‑148 (MQ=255)                                                  aaGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGt                                                                                                      >  1:330800/1‑148 (MQ=255)                                                  aaGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGt                                                                                                      >  4:70720/1‑148 (MQ=255)                                                   aGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTaa                                                                                                    <  6:109416/149‑1 (MQ=255)                                                   aGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTaa                                                                                                    <  6:296871/149‑1 (MQ=255)                                                        tGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAAc                                                                                               <  5:173415/149‑1 (MQ=255)                                                               gATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GCAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAg                                                                                        >  4:86012/1‑149 (MQ=255)                                                               gATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTA‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTgg                                                                                           >  7:163818/1‑146 (MQ=255)                                                                      gagaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgcgc                                                                                  <  8:163818/148‑1 (MQ=255)                                                                       agaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGaa                                                                                                                          >  2:373611/1‑107 (MQ=255)                                                                       agaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGaa                                                                                                                          <  1:373611/107‑1 (MQ=255)                                                                       agaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAcgcgcgc                                                                                <  7:74606/149‑1 (MQ=255)                                                                        gaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAc                                                                                      <  7:148356/142‑1 (MQ=255)                                                                        gaAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGAc                                                                                      >  8:148356/1‑142 (MQ=255)                                                                         aaaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACgcgcgcgc                                                                              >  5:122159/1‑149 (MQ=255)                                                                         aaaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACgcgcgcgc                                                                              >  1:241263/1‑149 (MQ=255)                                                                         aaaCTCAAAGCCCTGTTCACCGAGCAGGCCGAGAGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACgcgcgcgc                                                                              <  7:97414/149‑1 (MQ=255)                                                                           aCTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCg                                                                            <  8:124377/149‑1 (MQ=255)                                                                               aaaGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGAGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTa                                                                        >  6:333375/1‑149 (MQ=255)                                                                                   cccTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACctg                                                                    <  6:379513/149‑1 (MQ=255)                                                                                   cccTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGACCTACCAGGACAGCCTGACAGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAACACGCGCGCGCCGGTTACctg                                                                    >  6:93923/1‑149 (MQ=255)                                                                                    ccTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACctgc                                                                   <  8:143767/149‑1 (MQ=255)                                                                                        ttCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACctgctg                                                                 <  4:111699/147‑1 (MQ=255)                                                                                         tCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTtg                                                              >  1:10289/1‑149 (MQ=255)                                                                                           aCCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTtgct                                                            >  4:162686/1‑149 (MQ=255)                                                                                            ccGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACgcgcgcgc                                                                              <  2:196352/130‑1 (MQ=255)                                                                                            ccGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTtgctg                                                           >  4:261894/1‑149 (MQ=255)                                                                                             cGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTCTCTGGGAC‑C‑GC‑‑CGCTATTT‑‑TGAAAAGCAGCTGAACAATCGCGTAAGCAACCTGGAAGACGCGCGCGCCGCTTACCTGCTGTtgct                                                            >  1:118390/1‑149 (MQ=255)                                                                                                gCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCggg                                                       >  3:216051/1‑149 (MQ=255)                                                                                                      cgtgcgcAGTGACCGGCTGCCTGCGCACTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑TAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATTCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAgg                                                 <  2:149929/146‑1 (MQ=255)                                                                                                       gAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAggg                                                <  7:162915/149‑1 (MQ=255)                                                                                                          cgcCGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTCTTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGGTACCTGCTGTTGCTGCGGGGTCAGGGCct                                             <  3:86012/149‑1 (MQ=255)                                                                                                          cgcAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCct                                             <  1:1504/149‑1 (MQ=255)                                                                                                           gcAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGGTGTC‑‑GAACTCAGCATCCCTTTTTACTGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAg                                                  <  1:311596/149‑1 (MQ=255)                                                                                                                 aGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAACACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCtggctggc                                      <  5:285080/149‑1 (MQ=255)                                                                                                                   cGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTc                                                    >  8:336837/1‑133 (MQ=255)                                                                                                                   cGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTc                                                    <  7:336837/133‑1 (MQ=255)                                                                                                                   cGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCt                                    <  1:147178/149‑1 (MQ=255)                                                                                                                           caggctgTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgccc                            >  8:201768/7‑149 (MQ=255)                                                                                                                            gcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACct                                                                     >  5:354091/1‑107 (MQ=255)                                                                                                                            gcgGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACct                                                                     <  6:354091/107‑1 (MQ=255)                                                                                                                                 gTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTg                      >  1:116080/1‑149 (MQ=255)                                                                                                                                 gTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGGTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTg                      <  7:162351/149‑1 (MQ=255)                                                                                                                                  tCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgccc                            <  8:192516/142‑1 (MQ=255)                                                                                                                                  tCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgccc                            >  7:192516/1‑142 (MQ=255)                                                                                                                                     tACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcgg                  >  7:93955/1‑149 (MQ=255)                                                                                                                                      aCCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTggcggc                 >  1:234535/1‑149 (MQ=255)                                                                                                                                         tcgACAGCCCGGTCAGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTAGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAg              >  3:116473/3‑149 (MQ=255)                                                                                                                                          ggACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGc             >  1:376853/1‑149 (MQ=255)                                                                                                                                          cgACAGCCCGGTCAGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGc             >  2:359315/2‑149 (MQ=255)                                                                                                                                           gACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCg            >  7:57248/1‑149 (MQ=255)                                                                                                                                              aGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTAc         <  7:201768/149‑1 (MQ=255)                                                                                                                                              aGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTAc         <  2:10289/149‑1 (MQ=255)                                                                                                                                                ccTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCTACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCg       >  2:148406/1‑149 (MQ=255)                                                                                                                                                ccTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTCCCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgc                              >  8:90498/1‑126 (MQ=255)                                                                                                                                                ccTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAAcgc                              <  7:90498/126‑1 (MQ=255)                                                                                                                                                ccTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTCGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTAcc        >  2:318254/1‑148 (MQ=255)                                                                                                                                                ccTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATAT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCg       <  4:245143/149‑1 (MQ=255)                                                                                                                                                 cTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGa      >  4:198007/1‑149 (MQ=255)                                                                                                                                                   gACCGGGCTGCC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGt                                                                >  4:29212/1‑89 (MQ=255)                                                                                                                                                   gACCGGGCTGCC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGt                                                                <  3:29212/89‑1 (MQ=255)                                                                                                                                                    aCCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGAc     <  5:93923/147‑1 (MQ=255)                                                                                                                                                    aCCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGAc     <  4:216051/147‑1 (MQ=255)                                                                                                                                                    aCCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGAc     <  2:329651/147‑1 (MQ=255)                                                                                                                                                    aCCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCa   <  8:53858/149‑1 (MQ=255)                                                                                                                                                     ccGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAg  <  5:333375/149‑1 (MQ=255)                                                                                                                                                     ccGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAg  <  3:162686/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                             CTCAGCCTGCCAGAGTTGCCACGCACGCCGGAGCTACGCCGCGTAGTGCAAGCGATGAACCAGATGGTCGAGAAACTCAAAGCCCTGTTCACCGAGCAGGCCGAGCGCAGTGAGCGGCTGCGTGCGGAGTCCTACCAGGACAGCCTGACCGGGCTGTC‑‑GAAC‑C‑GC‑‑CGCTATTT‑‑TGAAATGCAGCTGAACAATCGGGTAAGCAACCTGGAAGACGCGCGCGCCGGTTACCTGCTGTTGCTGCGGGTTCAGGGCCTGGCTGGCCTGAACGCCCGCCTTGGCGGCCAGCGTACCGACCAG  >  NC_002947/188468‑188764

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑