Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,248,512	C→G	5.1%	L505L (CTG→CTC)	PP_3724 ←	acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,248,512	0	C	G	5.1%	260.2 / 3.5	79	L505L (CTG→CTC)	PP_3724	acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (32/43);  new base G (2/2);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

GTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTC  >  NC_002947/4248367‑4248659                                                                                                                                                            |                                                                                                                                                      gTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAg                                                                                                                                                     >  8:238708/1‑148 (MQ=255) gTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGc                                                                                                                                                    <  5:185448/149‑1 (MQ=255)        ttCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCg                                                                                                                                             >  3:90194/1‑149 (MQ=255)        ttCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCg                                                                                                                                             >  3:90778/1‑149 (MQ=255)          ccAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCCTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCg                                                                                                                                             <  8:206423/147‑1 (MQ=255)              ggCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGCTGCGTTCGGGCGTACTCAGCACGCCCGCAGGc                                                                                                                                        <  1:260816/149‑1 (MQ=255)                  cgGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑gg                                                                                                                                 <  2:247572/149‑1 (MQ=255)                      cATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGa                                                                                                                             <  4:90194/149‑1 (MQ=255)                      cATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGa                                                                                                                             >  6:14113/1‑149 (MQ=255)                      cATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGa                                                                                                                             >  2:69015/1‑149 (MQ=255)                       aTCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGAc                                                                                                                            <  7:238708/149‑1 (MQ=255)                        tCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGAcc                                                                                                                           <  4:178456/149‑1 (MQ=255)                           cgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGa                                                                                                                             <  2:32052/144‑1 (MQ=255)                           cgcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGa                                                                                                                             >  1:32052/1‑144 (MQ=255)                            gcAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAg                                                                                                                       <  7:99857/149‑1 (MQ=255)                              aGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCa                                                                                                                                                   <  6:42056/121‑1 (MQ=255)                              aGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCa                                                                                                                                                   >  5:42056/1‑121 (MQ=255)                              aGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGgc                                                                                                                     <  1:121985/149‑1 (MQ=255)                               ggTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGgcg                                                                                                                    <  2:323322/149‑1 (MQ=255)                                 ttCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGgcgcg                                                                                                                  >  5:366426/1‑149 (MQ=255)                                 ttCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGCC‑GCGTTCGGGCGTACTCAGCACGCCCGCACGCACC‑‑GGCACACCGAAGgcgcg                                                                                                                  >  5:299125/1‑149 (MQ=255)                                    aaCGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGTGGAAGATGTCAGGGCACCCACGCC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGCCACC‑‑GGCAGACCGac                                                                                                                        >  1:359452/1‑148 (MQ=255)                                    aaCGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAg                                                                                                               >  6:305053/1‑149 (MQ=255)                                      cGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAggg                                                                                                             <  6:326041/149‑1 (MQ=255)                                        ccACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGcc                                                                                                           <  3:167700/149‑1 (MQ=255)                                          acaGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTg                                                                                                         >  8:128231/1‑149 (MQ=255)                                           caGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGc                                                                                                        <  2:58937/149‑1 (MQ=255)                                           caGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGc                                                                                                        <  6:366426/149‑1 (MQ=255)                                             tgaCCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACGGCCTTCAGTAGCTAGGGGTAACCGACGCTGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACTGC‑GCGTTCGGGCGTACTCAGCACGCCGGCAGGCACC‑‑GGGAGACCTAAGGCGCGCAGGGCCTGCtt                                                                                                      <  6:362610/147‑1 (MQ=255)                                             gCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCtt                                                                                                      <  8:309549/149‑1 (MQ=255)                                               aCCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTg                                                                                                     >  5:44021/1‑148 (MQ=255)                                                        ggTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGAt                                                                                                  <  5:65936/142‑1 (MQ=255)                                                        ggTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGAt                                                                                                  >  6:65936/1‑142 (MQ=255)                                                         gTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTc                                                                                          <  6:337909/149‑1 (MQ=255)                                                         gTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTc                                                                                          <  3:120971/149‑1 (MQ=255)                                                            ttGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTc                                                                                       >  1:77654/1‑149 (MQ=255)                                                            ttGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTc                                                                                       >  7:214464/1‑149 (MQ=255)                                                              gtgCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGAGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTcc                                                                                      <  3:234834/148‑1 (MQ=255)                                                               tgCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCTGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCaa                                                                                    <  1:345770/149‑1 (MQ=255)                                                                 cGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCaa                                                                                    >  1:396029/1‑147 (MQ=255)                                                                    cAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTc                                                                               <  2:34315/149‑1 (MQ=255)                                                                     aGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCg                                                                              <  2:233834/149‑1 (MQ=255)                                                                                  cggCCTTCAGTACGAAGGGGTAAGCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGGTTGTATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTg                                                                 <  1:237571/147‑1 (MQ=255)                                                                                     ccTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAgg                                                              <  2:61672/149‑1 (MQ=255)                                                                                     ccTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAgg                                                              >  3:185401/1‑149 (MQ=255)                                                                                               cAAGGGGTAACAGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCggg                                                     >  6:112723/1‑148 (MQ=255)                                                                                               cAAGGGGAAACCGACGCGGGCCG‑‑‑‑CCTCTG‑‑‑‑‑CCACGGC‑GCCTTCGGCCGAACTCAGCACGCCCCCAAGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCACCTTCGAATCCCATAAGTCCAACACTCGGCAATACCCCGTCAAGGCCAGCggg                                                     >  8:149551/1‑148 (MQ=255)                                                                                                 aGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGca                                                   >  7:58310/1‑148 (MQ=255)                                                                                                 aGGGGTAACCAACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGca                                                   >  3:139809/1‑148 (MQ=255)                                                                                                  ggggTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGcaca                                                 <  2:396029/149‑1 (MQ=255)                                                                                                   gggTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAg                                                <  5:14113/149‑1 (MQ=255)                                                                                                    ggTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGc                                               <  6:44021/149‑1 (MQ=255)                                                                                                     gTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCg                                              <  2:282860/149‑1 (MQ=255)                                                                                                          cGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTAcc                                         <  5:305053/149‑1 (MQ=255)                                                                                                          cGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTAcc                                         >  6:15721/1‑149 (MQ=255)                                                                                                            aCGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCg                                       >  5:439000/1‑149 (MQ=255)                                                                                                            aCGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCg                                       <  4:185401/149‑1 (MQ=255)                                                                                                                   cgg‑‑‑‑cgTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGc                                      <  3:188480/143‑1 (MQ=255)                                                                                                                   cgg‑‑‑‑cgTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGc                                      >  4:188480/1‑143 (MQ=255)                                                                                                                    gg‑‑‑‑cgTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCt                               <  2:359452/149‑1 (MQ=255)                                                                                                                     g‑‑‑‑cgTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTc                              <  5:68643/149‑1 (MQ=255)                                                                                                                     g‑‑‑‑cgTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTc                              >  4:204965/1‑149 (MQ=255)                                                                                                                     g‑‑‑‑cgTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTc                              >  5:302573/1‑149 (MQ=255)                                                                                                                     g‑‑‑‑cgTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTc                              <  2:89601/149‑1 (MQ=255)                                                                                                                                     ccACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTg                        <  5:15721/148‑1 (MQ=255)                                                                                                                                         ggC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCa                   <  2:77654/149‑1 (MQ=255)                                                                                                                                           gtgcctgcGGGCGTGCTGAGTACGCCCGAACGCGCCGTGGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTcc                    <  5:297999/142‑1 (MQ=37)                                                                                                                                           c‑gcgTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATc                 <  4:139809/149‑1 (MQ=255)                                                                                                                                             gcctgcgggcgggCTGAGTACGCCCGAACGCGCC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTAccc                                        >  6:439220/6‑125 (MQ=38)                                                                                                                                             gcctgcGGGCGAGCTGAGTACGCCCGAACGCGCC‑‑GGCAGACCGACGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTAccc                                        <  5:439220/120‑1 (MQ=37)                                                                                                                                                  cGGGCGTGCTGAGTACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGc           >  2:34616/1‑149 (MQ=255)                                                                                                                                                       gTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCg      <  7:149551/149‑1 (MQ=255)                                                                                                                                                       gTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCg      <  3:204965/149‑1 (MQ=255)                                                                                                                                                         aCTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCg                                                       >  5:164547/1‑98 (MQ=255)                                                                                                                                                         aCTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCg                                                       <  6:164547/98‑1 (MQ=255)                                                                                                                                                         aCTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGt    >  3:243668/1‑149 (MQ=255)                                                                                                                                                          cTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGt    >  4:356146/1‑148 (MQ=255)                                                                                                                                                           tCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTc  >  3:289947/1‑149 (MQ=255)                                                                                                                                                           tCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTc  >  5:95207/1‑149 (MQ=255)                                                                                                                                                            |                                                                                                                                                      GTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGG‑‑‑‑CGTCTG‑‑‑‑‑CCACGGC‑GCGTTCGGGCGTACTCAGCACGCCCGCAGGCACC‑‑GGCAGACCGAAGGCGCGCAGGGCCTGCTTGGATTGCCATTCGTCCAACAGTCGGCTTTCGCCCGTGAGGGCCTGCGGGCACAGCGGTACCCGCACAGCCTCGCCCTGGTCCATCAACCGCTGGCGGTTC  >  NC_002947/4248367‑4248659

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