Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,297,128	G→A	5.6%	L79L (CTG→TTG)	TpbB ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,297,128	0	G	A	5.6%	272.5 / 3.0	88	L79L (CTG→TTG)	TpbB	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (35/48);  new base A (2/3);  total (37/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ACCGGCTCATCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCC  >  NC_002947/5296980‑5297275                                                                                                                                                     |                                                                                                                                                    aCCGGCTCATCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAg                                                                                                                                                     >  4:203364/1‑149 (MQ=255)         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                                                                                                                              >  8:192206/1‑148 (MQ=255)                         gCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGt                                                                                                                              <  7:11053/148‑1 (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                                                        <  7:239414/147‑1 (MQ=255)                                                                                                gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGc                                                        <  3:242329/147‑1 (MQ=255)                                                                                                gttgttGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCAt                                                      <  5:156655/149‑1 (MQ=255)                                                                                                 ttgttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTa                                                                                                      <  4:2717/100‑1 (MQ=255)                                                                                                 ttgttgAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTa                                                                                                      >  3:2717/1‑100 (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)                                                                                                                                                     |                                                                                                                                                    ACCGGCTCATCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGGCCACTTCTTCGGTACTGGCAATCAACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAGGCC  >  NC_002947/5296980‑5297275

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