Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,297,136	T→C	5.8%	N76S (AAC→AGC)	TpbB ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,297,136	0	T	C	5.8%	247.9 / 4.2	86	N76S (AAC→AGC)	TpbB	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (36/44);  new base C (2/3);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

ATCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAG  >  NC_002947/5296988‑5297284                                                                                                                                                     |                                                                                                                                                     aTCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGt                                                                                                                                                      <  5:101999/149‑1 (MQ=255) aTCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGt                                                                                                                                                      >  1:253922/1‑149 (MQ=255)  tCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCTCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGtt                                                                                                                                                     <  8:226292/149‑1 (MQ=255)  tCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGtt                                                                                                                                                     <  8:222065/149‑1 (MQ=255)  tCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGtt                                                                                                                                                     <  2:188579/149‑1 (MQ=255)    cAGCAAGGCGCTGGCAACCTGCACCCCCAGCTGCGCGAGCATGCCGGTATTGTCTCGTTGCCAGGGCGGCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTg                                                                                                                                                    <  2:298367/148‑1 (MQ=255)     aGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAACGAATCGCTggc                                                                                                                                                  >  1:123453/1‑149 (MQ=255)      gCATGGCGCTGGCAACCTGCTCCTCCAGCTGGGCGAGCATGCCGGTATTGTCGGGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACGTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCTTTggcg                                                                                                                                                 <  2:244604/149‑1 (MQ=255)      gCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTggcg                                                                                                                                                 >  2:17400/1‑149 (MQ=255)       cAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTggcgg                                                                                                                                                >  8:11053/1‑149 (MQ=255)       cAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAACGAATCGCTggcgt                                                                                                                                                <  2:123453/149‑2 (MQ=255)       cAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAACGAATCGCTggcgt                                                                                                                                                <  2:254237/149‑2 (MQ=255)        aaGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAACGAATCGCTggcgtt                                                                                                                                               <  3:122277/149‑3 (MQ=255)          ggCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTggcggcgg                                                                                                                                             >  4:277221/1‑149 (MQ=255)            cgcTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCTGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGgcggcggcg                                                                                                                                           <  4:95451/149‑1 (MQ=255)             gcTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGgcgc                                                                                                                                          >  2:135090/1‑149 (MQ=255)                 gCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGt                                                                                                                                       <  7:11053/148‑1 (MQ=255)                 gCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGt                                                                                                                                       >  8:192206/1‑148 (MQ=255)                    aCCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCTTCGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATc                                                                                                                                   <  1:273084/149‑1 (MQ=255)                          aCCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGATGACAGCGCCGCCGCCAACGaatcg                                                                                                                             <  1:323720/149‑4 (MQ=255)                           ccTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAAcacc                                                                                                                            <  1:37178/149‑1 (MQ=255)                               cAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCTAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAg                                                                                                                        <  7:27300/149‑1 (MQ=255)                                      gcgAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTc                                                                                                                 >  7:190501/1‑149 (MQ=255)                                        gAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTa                                                                                                               >  6:156655/1‑149 (MQ=255)                                              gCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGgtg                                                                                                          >  8:89285/1‑148 (MQ=255)                                                cGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAg                                                                                                       >  3:42681/1‑149 (MQ=255)                                                     ttGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGAt                                                                                                  <  8:9622/149‑1 (MQ=255)                                                        tCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGa                                                                                               <  1:59195/149‑1 (MQ=255)                                                         cgcgTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAg                                                                                              >  2:309259/1‑149 (MQ=255)                                                          gcgTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCACCGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCTCAGCCGCTTCTACGGTGTAGTTGATCGAGc                                                                                             >  8:295511/1‑149 (MQ=255)                                                          gcgTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGc                                                                                             <  4:42681/149‑1 (MQ=255)                                                           cgTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCg                                                                                            >  6:60582/1‑149 (MQ=255)                                                            gTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCgg                                                                                           <  3:128508/149‑1 (MQ=255)                                                                 cAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGTAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGAt                                                                                      >  3:309821/1‑149 (MQ=255)                                                                 cAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAACACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCTAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGAt                                                                                      <  8:56422/149‑1 (MQ=255)                                                                      cgcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAAc                                                                                  <  6:178887/148‑1 (MQ=255)                                                                      cgcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACt                                                                                 <  7:89285/149‑1 (MQ=255)                                                                       gcCAGTTGCTCGCCGTCGATGTGGAAGAGGTTGGCCTCCGCCTCTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGGGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGCCCGCTCGATCAACta                                                                                <  4:249067/149‑2 (MQ=255)                                                                        ccAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGc                                                                               <  2:33657/149‑1 (MQ=255)                                                                        ccAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGc                                                                               <  2:28683/149‑1 (MQ=255)                                                                          aGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAg                                                                             <  7:192206/149‑1 (MQ=255)                                                                           gTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATc                                                                                                 <  3:226977/128‑1 (MQ=255)                                                                           gTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATc                                                                                                 >  4:226977/1‑128 (MQ=255)                                                                           gTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAgg                                                                            <  1:135090/149‑1 (MQ=255)                                                                              gCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTa                                                                         >  3:33809/1‑149 (MQ=255)                                                                              gCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTa                                                                         >  4:49563/1‑149 (MQ=255)                                                                               cTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTAt                                                                        >  6:245809/1‑149 (MQ=255)                                                                                tCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGTCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATg                                                                       >  5:153340/1‑149 (MQ=255)                                                                                  gCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATgg                                                                      <  1:309259/148‑1 (MQ=255)                                                                                    ccTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                                   >  5:319702/1‑149 (MQ=255)                                                                                     cacGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAgg                                                                            <  1:275369/137‑1 (MQ=255)                                                                                     cTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGTCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                  <  6:153340/149‑1 (MQ=255)                                                                                     cTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAgg                                                                            >  2:275369/1‑139 (MQ=255)                                                                                     cTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                  >  4:242329/1‑149 (MQ=255)                                                                                     cTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGCCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                                  <  3:203364/149‑1 (MQ=255)                                                                                      tCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAAAGCCAGCGATTCGTTGGCGGCGGCCCTGACATCGAACACCACAGCCGCTTCTACGGTGTAGATGATCGAGCGTGCGATCAACTGCAGGTTATGGTTgc                                                                  >  6:65446/1‑147 (MQ=255)                                                                                        gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGc                                                                 <  3:242329/147‑1 (MQ=255)                                                                                        gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGc                                                                 <  7:239414/147‑1 (MQ=255)                                                                                        gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCAt                                                               <  5:156655/149‑1 (MQ=255)                                                                                         ttgttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTa                                                                                                               >  3:2717/1‑100 (MQ=255)                                                                                         ttgttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTa                                                                                                               <  4:2717/100‑1 (MQ=255)                                                                                               aaGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGacac                                                         <  2:253922/148‑1 (MQ=255)                                                                                                  aCCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAg                                                     <  3:277221/149‑1 (MQ=255)                                                                                                        gCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGa                                                <  3:132740/148‑1 (MQ=255)                                                                                                         ccTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAg                                               >  2:60692/1‑148 (MQ=255)                                                                                                               ccACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGccc                                        <  7:238344/149‑1 (MQ=255)                                                                                                               ccACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGcc                                         >  8:238344/1‑148 (MQ=255)                                                                                                                 aCTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcag                                      >  3:302022/1‑149 (MQ=255)                                                                                                                      ttcGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAggg                                 >  3:225145/1‑149 (MQ=255)                                                                                                                      ttcGGTACGGGCAATCAACGCCTGCGATTCGTTGGCGTCGCCGCTGGCATCGACCACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGGTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAggg                                 <  5:65446/149‑1 (MQ=255)                                                                                                                        cGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                                  <  5:245809/146‑1 (MQ=255)                                                                                                                        cGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGt                                <  1:60692/148‑1 (MQ=255)                                                                                                                          gTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                             >  3:177794/1‑149 (MQ=255)                                                                                                                           tACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                            <  3:49563/149‑1 (MQ=255)                                                                                                                           tACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                            <  8:190501/149‑1 (MQ=255)                                                                                                                             cTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGTGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAggg                           <  1:175845/148‑1 (MQ=255)                                                                                                                              tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGc                                                  >  7:184475/1‑124 (MQ=255)                                                                                                                              tGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGc                                                  <  8:184475/124‑1 (MQ=255)                                                                                                                               ggCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                            >  3:24654/1‑145 (MQ=255)                                                                                                                                  aaTCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAc                     <  4:33809/149‑1 (MQ=255)                                                                                                                                    tCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAccc                   >  4:74940/1‑149 (MQ=255)                                                                                                                                    tCAACGCCAGCGATTCCTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACTTTGTAGTTGATCGAGAGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAccc                   <  1:268113/149‑1 (MQ=255)                                                                                                                                    ccgccGCCAACGAATCGCTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATAccc                   >  4:99975/5‑149 (MQ=255)                                                                                                                                        cGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGc                 >  7:110535/1‑147 (MQ=255)                                                                                                                                           cAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCa                                                      <  3:314482/107‑1 (MQ=255)                                                                                                                                           cAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCa                                                      >  4:314482/1‑107 (MQ=255)                                                                                                                                           cAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGc            >  2:189447/1‑149 (MQ=255)                                                                                                                                             gCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGcc           >  4:235263/1‑148 (MQ=255)                                                                                                                                             gCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGcc           >  7:164304/1‑148 (MQ=255)                                                                                                                                              cgtttCGTTGGCGGCGGCGCTGTCATCGAACACGACAGCCGCTTCTACGGTGTTGTTGATCGAGCGGGGGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCaa         <  6:331892/146‑1 (MQ=255)                                                                                                                                                aTTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACg       <  3:235263/149‑1 (MQ=255)                                                                                                                                                     ttGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGgcaag  <  4:309821/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACGGCAAG  >  NC_002947/5296988‑5297284

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