Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,457,559	T→A	6.7%	I378F (ATC→TTC)	bhbP ←	D‑beta‑hydroxybutyrate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,457,559	0	T	A	6.7%	173.7 / 2.9	60	I378F (ATC→TTC)	bhbP	D‑beta‑hydroxybutyrate permease
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base A (2/2);  total (30/30)
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.60e-01

ACCGCCAGCAAGGTGATCACCGCGCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTA  >  NC_002947/3457411‑3457706                                                                                                                                                      |                                                                                                                                                    aCCGCCAGCAAGGTGATCACCGCGCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑t                                                                                                                                                     <  6:246464/149‑1 (MQ=255) aCCGCCAGCAAGGTGATCACCGCGCCATTGTGCAGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑t                                                                                                                                                     >  1:49179/1‑149 (MQ=255)   cGCCAGCAAGGTGATCACCGCGCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGc                                                                                                                                                   <  1:329674/149‑1 (MQ=255)   cGCCAGCAAGGTGATCACCGCGCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGc                                                                                                                                                   <  1:333918/149‑1 (MQ=255)        gCAAGGTGATCACCGCGCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATg                                                                                                                                              <  8:232510/149‑1 (MQ=255)                   aCCGCGCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAg                                                                                                                                    <  8:5265/148‑1 (MQ=255)                           aTTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCgg                                                                                                                           <  6:266104/149‑1 (MQ=255)                            ttGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGt                                                                                                                          <  4:16031/149‑1 (MQ=255)                            ttGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGt                                                                                                                          <  3:249149/149‑1 (MQ=255)                                  ggCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGcc                                                                                                                    >  3:342250/1‑149 (MQ=255)                                     aGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGc                                                                                                                 <  8:42660/149‑1 (MQ=255)                                     aGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGc                                                                                                                 <  2:49179/149‑1 (MQ=255)                                              aTGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑AGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGt                                                                                                        >  4:322196/1‑149 (MQ=255)                                              aTGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑AGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGt                                                                                                        >  4:327561/1‑149 (MQ=255)                                                   aCCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑GGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACgg                                                                                                   >  7:27591/1‑149 (MQ=255)                                                        tGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATg                                                                                              >  8:258912/1‑149 (MQ=255)                                                           ccATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCt                                                                                           <  4:260958/149‑1 (MQ=255)                                                             aTGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGa                                                                                                                                <  5:290855/110‑1 (MQ=255)                                                             aTGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGa                                                                                                                                >  6:290855/1‑110 (MQ=255)                                                                  ggcTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCaa                                                                                                           >  3:124238/1‑126 (MQ=255)                                                                  ggcTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCaa                                                                                                           <  4:124238/126‑1 (MQ=255)                                                                     tACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCAcc                                                                                  >  5:213069/1‑148 (MQ=255)                                                                      aCCCGGTGTAGCACTTCCAGCTGGATATCGGCCGCGTGGGCGGCGGAGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑GGCTCATGCCGCCCTAAGCGGAGCCGGTGATGCCGGCCAGCAAGGGCACGGTGATGGCTTCGTTCACCaa                                                                                <  7:131647/149‑1 (MQ=255)                                                                       cccGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAg                                                                               >  5:242248/1‑149 (MQ=255)                                                                       cccGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAg                                                                               >  5:236773/1‑149 (MQ=255)                                                                        ccGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAgg                                                                              >  4:175642/1‑149 (MQ=255)                                                                        ccGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAgg                                                                              >  7:119640/1‑149 (MQ=255)                                                                        ccGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGAGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGATTCTTTCACCAAgg                                                                              >  6:184797/1‑149 (MQ=255)                                                                              tAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTg                                                                         <  2:191180/148‑1 (MQ=255)                                                                               aGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTggg                                                                       >  6:113795/1‑149 (MQ=255)                                                                                gCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTggg                                                                       <  2:113671/148‑1 (MQ=255)                                                                                gCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTggg                                                                       >  1:98694/1‑148 (MQ=255)                                                                                 cACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGTGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGa                                                                     >  1:269635/1‑149 (MQ=255)                                                                                 cACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGa                                                                     >  4:302278/1‑149 (MQ=255)                                                                                  aCTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGAt                                                                    >  3:186730/1‑149 (MQ=255)                                                                                   cTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGAtg                                                                   <  4:186730/149‑1 (MQ=255)                                                                                   cTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGAtg                                                                   <  6:213069/149‑1 (MQ=255)                                                                                         gCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTgcg                                                              >  3:321466/1‑148 (MQ=255)                                                                                          cGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑AGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCa                                                            <  3:327561/149‑1 (MQ=255)                                                                                          cGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑AGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCa                                                            <  3:322196/149‑1 (MQ=255)                                                                                                atCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCAt                                                      >  7:184372/1‑149 (MQ=255)                                                                                                  cGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCg                                                    >  3:222612/1‑149 (MQ=255)                                                                                                     ccGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGcc                                                 <  5:113795/149‑1 (MQ=255)                                                                                                          tGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTgg                                                                        <  1:104305/121‑1 (MQ=255)                                                                                                          tGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTgg                                                                        >  2:104305/1‑121 (MQ=255)                                                                                                            ggcggcggcGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCa                                          <  4:342250/149‑1 (MQ=255)                                                                                                               ggcggcgATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGa                                        >  6:328458/1‑148 (MQ=255)                                                                                                                    cgATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAAccc                                  <  4:222612/149‑1 (MQ=255)                                                                                                                      aTGAACGTGTCGCCCATGGCTGCCAAGGCGATTGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCg                                 >  1:39349/1‑149 (MQ=255)                                                                                                                       tGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCgg                                >  7:262936/1‑148 (MQ=255)                                                                                                                         aaCGTTTCGCCCATGGCTGCCAAGGCTA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGACGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCTGCCAGCACCAGAAAACCCGGCa                              <  5:184797/148‑1 (MQ=255)                                                                                                                              gtCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAgg                        >  3:147502/1‑149 (MQ=255)                                                                                                                               tCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGc                       <  2:269635/149‑1 (MQ=255)                                                                                                                               tCGCCCATGACTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTggg                                                                       <  7:216940/101‑1 (MQ=255)                                                                                                                               tCGCCCATGACTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTggg                                                                       >  8:216940/1‑101 (MQ=255)                                                                                                                                  cccATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGAt                    <  2:98694/149‑1 (MQ=255)                                                                                                                                        gCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACcgc              <  6:236773/149‑1 (MQ=255)                                                                                                                                         cTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGTCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACcgc              <  6:242248/148‑1 (MQ=255)                                                                                                                                               aGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCaacg       >  6:19539/1‑147 (MQ=255)                                                                                                                                               aGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAg       <  4:134997/149‑1 (MQ=255)                                                                                                                                               aGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAg       >  3:134997/1‑149 (MQ=255)                                                                                                                                                    a‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTa  >  3:122378/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ACCGCCAGCAAGGTGATCACCGCGCCATTGTGCGGCAGGGTGTCCATGCCACCGCTGGCCATGGCGGCTACCCGGTGTAGCACTTCCAGCGGGATATCGGCCGCGTGGGCGGCGGCGATGAACGTGTCGCCCATGGCTGCCAAGGCGA‑TGCTCATGCCGCCCGAAGCGGAGCCGGTGATGCCGGCCAGCAAGGTCACGGTGATGGCTTCGTTCACCAAGGGGTTGGGGATGCTGCGCAGGGCATCGGCCAGCACCAGAAAACCCGGCAGCGAGGCGATTACCGCGCCAAAGCCGTA  >  NC_002947/3457411‑3457706

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