Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,297,132	T→A	5.5%	E77D (GAA→GAT)	TpbB ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,297,132	0	T	A	5.5%	281.1 / 3.5	91	E77D (GAA→GAT)	TpbB	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (37/49);  new base A (2/3);  total (39/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

ATCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACG  >  NC_002947/5296988‑5297279                                                                                                                                                 |                                                                                                                                                    aTCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTc                                                                                                                                                   <  6:324713/147‑1 (MQ=255) 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:28683/149‑1 (MQ=255)                                                                        ccAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGc                                                                          <  2:33657/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                                                            <  7:239414/147‑1 (MQ=255)                                                                                        gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGc                                                            <  3:242329/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      >  7:164304/1‑148 (MQ=255)                                                                                                                                             gCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGcc      >  4:235263/1‑148 (MQ=255)                                                                                                                                              cgtttCGTTGGCGGCGGCGCTGTCATCGAACACGACAGCCGCTTCTACGGTGTTGTTGATCGAGCGGGGGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCaa    <  6:331892/146‑1 (MQ=255)                                                                                                                                                aTTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACg  <  3:235263/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    ATCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCCAACG  >  NC_002947/5296988‑5297279

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