Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,133,979	C→G	8.1%	A351P (GCA→CCA) ‡	PP_3638 ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,133,979	0	C	G	8.1%	275.0 / 10.8	86	A351P (GCA→CCA) ‡	PP_3638	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (41/38);  new base G (5/2);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

CACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTA  >  NC_002947/4133831‑4134122                                                                                                                                                       |                                                                                                                                                  cacaCCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACgcct‑gc                                                                                                                                                   >  2:178913/1‑149 (MQ=255)  acacCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACgcct‑gcc                                                                                                                                                  >  7:396760/1‑149 (MQ=255)            ggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCa                                                                                                                                      <  6:431512/149‑1 (MQ=255)                 aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCg                                                                                                                                 <  7:343618/149‑1 (MQ=255)                   gATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCA‑GGCAGGC‑‑GTTGc                                                                                                                                       >  1:58567/1‑142 (MQ=255)                   gATGCTCACGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCA‑GGCAGGC‑‑GTTGc                                                                                                                                       <  2:58567/141‑1 (MQ=255)                      gCTCACGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGgcgcg                                                                                                                            >  7:73547/1‑149 (MQ=255)                      gCTCACGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGgcgc                                                                                                                             >  7:69523/1‑148 (MQ=255)                       cTCACGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGgcgcgc                                                                                                                           <  5:87404/149‑1 (MQ=255)                          aCGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑g                                                                                                                                           <  8:62701/130‑1 (MQ=255)                          aCGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑g                                                                                                                                           >  7:62701/1‑130 (MQ=255)                          aCGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCaa                                                                                                                         <  8:119802/148‑1 (MQ=255)                             aCCTCACCTGGTGCCAGCTTAATG‑TCGATCTTCGCCAGCACCTGCTTCCAACCACCCTCACGCTGGTATCCTAAGCGGATGTCACGGGCCTCCAGCATCACCGCGCTCATGTAACTCCT‑GCCTGTC‑‑GTTGCAGTTCGGCGCGCAATTg                                                                                                                      <  2:390147/148‑1 (MQ=255)                              ccTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCAGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCa                                                                                                                    >  5:50531/1‑149 (MQ=255)                               cTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCAc                                                                                                                   >  7:120615/1‑149 (MQ=255)                                tCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCAc                                                                                                                   >  1:106147/1‑148 (MQ=255)                                  aCCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAg                                                                                                                >  5:140914/1‑149 (MQ=255)                                          ccAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                         <  3:153550/148‑1 (MQ=255)                                            aGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGt                                                                                                       >  3:323246/1‑148 (MQ=255)                                               tGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGACGt                                                                                                       >  5:126387/1‑145 (MQ=255)                                               tGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGACGt                                                                                                       <  6:126387/145‑1 (MQ=255)                                                 cAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACg                                                                                                   <  5:169940/147‑1 (MQ=255)                                                  aGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATc                                                                                                <  1:178913/149‑1 (MQ=255)                                                    g‑TCGAACTGCGCCAGCACCTGCTGCAAACCACCCTCCCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACCCCGCGCTCATGCAACGCCT‑GCCTCGC‑‑GTTGCAGTTAGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCgg                                                                                              >  6:445305/1‑149 (MQ=255)                                                       cGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCa                                                                                            <  3:77408/149‑1 (MQ=255)                                                        gAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACCCCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCa                                                                                            <  4:290337/148‑1 (MQ=255)                                                        gAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCg                                                                                               <  5:15125/145‑1 (MQ=255)                                                        gAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCg                                                                                               >  6:15125/1‑145 (MQ=255)                                                              cgcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGc                                                                                     >  3:108137/1‑149 (MQ=255)                                                                ccAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAgggtt                                                                                   >  8:228729/1‑146 (MQ=255)                                                                ccAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGcc                                                                                    >  4:427174/1‑148 (MQ=255)                                                                ccAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGCC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGcc                                                                                    >  8:124801/1‑148 (MQ=255)                                                                 cAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGa                                                                                  <  4:160264/149‑1 (MQ=255)                                                                    cACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCg                                                                                   <  7:164393/145‑1 (MQ=255)                                                                    cACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCg                                                                                   >  8:164393/1‑145 (MQ=255)                                                                      ccTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTccc                                                                             <  2:107970/149‑1 (MQ=255)                                                                       gtgctgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTGACGGGCCTCCAGCAACACGGCGCTCATGCAACGCCT‑GCCTGTC‑‑GTTGCAGTTCGGCGCGCAATTGCAACAGGCTTGGCGTGACGATCGGCACGAAGGCCGACTCCCg                                                                            <  1:71108/148‑1 (MQ=255)                                                                       ctgctgCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTccc                                                                             <  6:135531/148‑1 (MQ=255)                                                                       ctgctgCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTccc                                                                             >  5:135531/1‑148 (MQ=255)                                                                         gctgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATACAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGcc                                                                          >  3:381283/1‑149 (MQ=255)                                                                             ccAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcg                                                                      >  7:47856/1‑149 (MQ=255)                                                                               aaCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                                    >  6:307668/1‑149 (MQ=255)                                                                                accaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                                   <  3:172465/149‑1 (MQ=255)                                                                                  caccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCggg                                                                 <  8:324304/149‑1 (MQ=255)                                                                                  caccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCggg                                                                 <  3:427174/149‑1 (MQ=255)                                                                                   accCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATACAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCGATGGCggga                                                                <  4:381283/149‑2 (MQ=255)                                                                                     ccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCGATGGCgggagt                                                              <  8:120615/149‑4 (MQ=255)                                                                                       tCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGCCAGGTGCAGTTAGGCGCGCAATTGCACCAGGCAGGGCGTGACGATAGGCCCGAAGGCCGCCTACCGCCATCGGCCGGCa                                                               >  1:135661/1‑148 (MQ=255)                                                                                           gcgGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCtt                                                        >  4:11785/1‑149 (MQ=255)                                                                                             gggTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCAGGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCg                                                      >  1:96042/1‑149 (MQ=255)                                                                                               ggaGCCCAAGCTGAGGCGACGGGCCGCCAGCAACTCCGCGCTCATGCAACGCCT‑GCCTGGC‑‑TTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCAGCCATCGGCGGGAAAACCCTTCTcc                                                    <  8:197571/147‑1 (MQ=255)                                                                                               gTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                                    <  7:70386/149‑1 (MQ=255)                                                                                               gTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                                    <  8:396760/149‑1 (MQ=255)                                                                                                tAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCg                                                   >  4:332377/1‑149 (MQ=255)                                                                                                   cccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAt                                                >  1:134176/1‑149 (MQ=255)                                                                                                    ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAtt                                               >  6:210410/1‑149 (MQ=255)                                                                                                    ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTACACTTCGGCGCGCAATTGCACCAGCCTGGGCATGACGATCGGCACGAAGGCCGACTCCCCCCATCCGCCGGCCAACCCTTCGCCGTAtt                                               >  7:155208/1‑149 (MQ=255)                                                                                                            aTGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAgg                                       <  8:47856/149‑1 (MQ=255)                                                                                                             tGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCTCCGTAt                                                <  3:290813/139‑1 (MQ=255)                                                                                                             tGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCACCGTAt                                                >  4:290813/1‑139 (MQ=255)                                                                                                                cACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAg                                   >  5:42575/1‑149 (MQ=255)                                                                                                                 aCGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCg                                             >  8:6482/1‑138 (MQ=255)                                                                                                                 aCGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCg                                             <  7:6482/138‑1 (MQ=255)                                                                                                                 aCGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGc                                  <  3:11785/149‑1 (MQ=255)                                                                                                                 aCGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCACGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGCCTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGc                                  >  6:354956/1‑149 (MQ=255)                                                                                                                   gggCCTCCAGCAACACCGCGCTCATGCAACGCCTGGCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGcc                                 >  1:284561/1‑149 (MQ=255)                                                                                                                     gCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTg                              >  3:70831/1‑149 (MQ=255)                                                                                                                      ccTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgc                             >  7:191346/1‑149 (MQ=255)                                                                                                                       cTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTg                                                                                                                                        <  8:147082/41‑1 (MQ=255)                                                                                                                       cTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgcc                            <  2:284561/149‑1 (MQ=255)                                                                                                                       cTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCAACGGCGGGCCAACCCTGCGCCGTCTTCGCAGCGGTAGGGCTTgcc                            >  8:348648/1‑149 (MQ=255)                                                                                                                        tCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGAAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccg                           >  1:330953/1‑149 (MQ=255)                                                                                                                           aGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCa                        >  5:194450/1‑149 (MQ=255)                                                                                                                            gCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                       <  6:42575/149‑1 (MQ=255)                                                                                                                              aaCACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTg                     >  4:390603/1‑149 (MQ=255)                                                                                                                              aaCACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCATATTCGCAGAGGTAAGCCTTGCCGCCACTg                     >  4:395345/1‑149 (MQ=255)                                                                                                                                cacCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGc                   <  3:332377/149‑1 (MQ=255)                                                                                                                                    gcgcTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGc               <  2:96042/149‑1 (MQ=255)                                                                                                                                    gccgaactGCAACGCCA‑GGCAGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTg                >  6:54053/8‑148 (MQ=255)                                                                                                                                          ctGCAACGCCA‑GGCAGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc          >  4:146435/2‑148 (MQ=255)                                                                                                                                           tGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCa        <  1:79081/149‑1 (MQ=255)                                                                                                                                           tGCAACGCCA‑GGCAGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGc                   <  2:237934/138‑1 (MQ=255)                                                                                                                                           tGCAACGCCA‑GGCAGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGc                   >  1:237934/1‑138 (MQ=255)                                                                                                                                            gCAACGCCA‑GGCAGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg       >  7:65707/1‑149 (MQ=255)                                                                                                                                             cAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGt      >  6:140177/1‑149 (MQ=255)                                                                                                                                              aaCGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGtt     >  2:44855/1‑149 (MQ=255)                                                                                                                                               aCGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCATATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTg    <  3:395345/149‑1 (MQ=255)                                                                                                                                                cGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGTGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGt   <  4:323246/149‑1 (MQ=255)                                                                                                                                                 gcct‑gcctGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTa  <  5:307668/149‑1 (MQ=255)                                                                                                                                                 gcct‑gcctGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTa  >  6:400409/1‑149 (MQ=255)                                                                                                                                                 gcct‑gcctGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTa  >  8:8443/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                  CACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGG‑TCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCT‑GCCTGGC‑‑GTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTA  >  NC_002947/4133831‑4134122

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