Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,793,892	A→T	7.4%	L211M (TTG→ATG)	PP_3355 ←	beta‑ketothiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,793,892	0	A	T	7.4%	191.4 / 4.1	66	L211M (TTG→ATG)	PP_3355	beta‑ketothiolase
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base T (2/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGC  >  NC_002947/3793747‑3794034                                                                                                                                                  |                                                                                                                                               gCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGgcaatg                                                                                                                                              >  1:155522/1‑148 (MQ=255) gCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGgcaatg                                                                                                                                              >  1:81664/1‑148 (MQ=255)     ccgccgCGCCATCGACCACCGCGCAACTGTTGGCAGCCGTCTGCAGGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGgcaatgcatt                                                                                                                                          <  4:58569/148‑3 (MQ=255)     ccgccgCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGgcattgcatt                                                                                                                                          <  6:237842/148‑3 (MQ=255)         cgcgCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCATTGCATTGCcgc                                                                                                                                     <  3:3500/149‑1 (MQ=255)         cgcgCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCcgc                                                                                                                                     <  4:294090/149‑1 (MQ=255)          gcgcCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCATTGCATTGCcgcg                                                                                                                                    <  7:201633/149‑1 (MQ=255)           cgcCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGt                                                                                                                                   <  2:163135/149‑1 (MQ=255)              cATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCATTGCATTGCCGCGtggc                                                                                                                                >  1:292828/1‑149 (MQ=255)                     caccGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGc                                                                                                                         >  8:68819/1‑149 (MQ=255)                     caccGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGc                                                                                                                         >  4:90734/1‑149 (MQ=255)                     caccGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGc                                                                                                                         <  3:90734/149‑1 (MQ=255)                        cgcgCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCtt                                                                                                                      <  8:35914/149‑1 (MQ=255)                          cgcAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCt                                                                                                                    <  6:42807/149‑1 (MQ=255)                              aCTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCt                                                                                                                >  5:156777/1‑149 (MQ=255)                                  ttGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTcc                                                                                                              <  2:330868/147‑1 (MQ=255)                                    gCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTg                                                                                                          >  6:52637/1‑149 (MQ=255)                                     ccAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGc                                                                                                         >  1:137553/1‑149 (MQ=255)                                          cGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAAAGCGTGGGTGGCTGTCCTGGCTGCCCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTgcg                                                                                                    >  1:98437/1‑149 (MQ=255)                                             gtGCACGCCGGTGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCgg                                                                                                 >  5:152362/2‑149 (MQ=255)                                               gCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGt                                                                                                >  6:322924/1‑148 (MQ=255)                                                  cGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACg                                                                                            <  4:114127/149‑1 (MQ=255)                                                       ccGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGAc                                                                                       <  2:81664/149‑1 (MQ=255)                                                         ggCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGa                                                                                     <  5:242920/149‑1 (MQ=255)                                                          gCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGa                                                                                     <  6:152362/148‑1 (MQ=255)                                                          gCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGAt                                                                                    <  2:246275/149‑1 (MQ=255)                                                                aCCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtc                                                                              >  6:334199/1‑149 (MQ=255)                                                                 ccGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCtcg                                                                             <  6:156777/149‑1 (MQ=255)                                                                          gTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAg                                                                    >  2:172768/1‑149 (MQ=255)                                                                           tCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTa                                                                                                                   >  7:35548/1‑101 (MQ=255)                                                                           tCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTa                                                                                                                   <  8:35548/101‑1 (MQ=255)                                                                           tCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGc                                                                   >  6:250737/1‑149 (MQ=255)                                                                           tCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGc                                                                   <  2:155522/149‑1 (MQ=255)                                                                             gggCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCGCTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGccc                                                                 >  3:202196/1‑149 (MQ=255)                                                                             gggCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGcc                                                                  <  5:52637/148‑1 (MQ=255)                                                                               gCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCACCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGcc                                                                  <  2:33924/146‑1 (MQ=255)                                                                               gCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCACCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGcc                                                                  >  1:33924/1‑146 (MQ=255)                                                                                ccAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTg                                                               >  1:310632/1‑148 (MQ=255)                                                                                   gTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGTCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGc                                                                   <  5:186981/141‑1 (MQ=255)                                                                                   gTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGc                                                                   >  6:186981/1‑141 (MQ=255)                                                                                       ggCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTCCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAAc                                                                       <  1:75026/133‑1 (MQ=255)                                                                                       cgCCGCGTCGGTAGGCCCTGCGTGGCTGTCCTGCCTGAACTCACTCACGCCACCCGACAATGCAATGCCGCGTGGCTGGCAGCCATCTAACTCCAGTGCCTGCCCGGTTACGCTCACGATTTCCTCGTCGAAc                                                                       >  2:75026/2‑133 (MQ=255)                                                                                        gCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGc                                                      <  5:250737/149‑1 (MQ=255)                                                                                           gcgTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCAAGCCACGCGGCATTGCATTGCCGCGt                                                                                                                                   <  4:71250/69‑1 (MQ=37)                                                                                           gcgTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCAAGCCACGCGGCATTGCATTGCCGCGt                                                                                                                                   >  3:71250/1‑69 (MQ=37)                                                                                                 gTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGc                                             >  6:285016/1‑149 (MQ=255)                                                                                                    ggCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCg                                          <  7:68819/149‑1 (MQ=255)                                                                                                    ggCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCg                                          >  5:74458/1‑149 (MQ=255)                                                                                                     gCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGc                                         >  8:98611/1‑149 (MQ=255)                                                                                                        gtgggtggCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGa                                      <  5:322924/149‑1 (MQ=255)                                                                                                         tgggtggCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGAt                                     >  5:308053/1‑149 (MQ=255)                                                                                                           ggtggCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATgg                                   <  7:98611/149‑1 (MQ=255)                                                                                                               gCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTgc                               >  1:319100/1‑149 (MQ=255)                                                                                                                 tGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGAGTGGATGGAAGCCTTCTACATCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCATGCCCATGGGATTGCAGGGCCCTCTGATGGCTgac                              >  6:205193/1‑146 (MQ=255)                                                                                                                      tGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGc                        >  6:45273/1‑149 (MQ=255)                                                                                                                       ggCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGcc                       >  8:31088/1‑149 (MQ=255)                                                                                                                            aCCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCAc                    <  6:308053/147‑1 (MQ=255)                                                                                                                              ctcactcaCGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGt                >  1:27228/1‑149 (MQ=255)                                                                                                                                 actcacGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTcc              <  5:285016/148‑1 (MQ=255)                                                                                                                                 actcacGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCa             >  6:306328/1‑149 (MQ=255)                                                                                                                                           cgcgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcg   >  2:317722/1‑149 (MQ=255)                                                                                                                                            gcgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGtt        >  1:246142/1‑143 (MQ=255)                                                                                                                                            gcgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGtt        <  2:246142/143‑1 (MQ=255)                                                                                                                                            gcgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCAcgc  <  2:137553/149‑1 (MQ=255)                                                                                                                                             cgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTgg                                                        >  4:52534/1‑94 (MQ=255)                                                                                                                                             cgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTgg                                                        <  3:52534/94‑1 (MQ=255)                                                                                                                                             cgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGccc                                           >  8:182690/1‑107 (MQ=255)                                                                                                                                             cgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGccc                                           <  7:182690/107‑1 (MQ=255)                                                                                                                                             cgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCAc                    <  6:286701/130‑1 (MQ=255)                                                                                                                                             cgGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCAc                    >  5:286701/1‑130 (MQ=255)                                                                                                                                                  |                                                                                                                                               GCTGCCGCCGCGCCATCGACCACCGCGCAACTGTTGCCAGCCGTCTGCACGCCACCGGCGTGCACCGGGCGCAGTCGGGCCAGTGCGGCCGCGTCGGTAGGCCGTGGGTGGCTGTCCTGGCTGACCTCACTCACGCCACGCGGCAATGCAATGCCGCGTGGCTGGCAGCCTTCTACCTCCAGTGCCTGCGCGGTGACGCTGACGATTTCCTCGTCGAACCAGCCCTGCACCTGGGCTTGCAGGGCCCGCTGATGGCTGCGCAGGGCCCACGCGTCCACCGTTTCACGC  >  NC_002947/3793747‑3794034

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