Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,375,222	A→C	5.3%	G73G (GGT→GGG)	dnaJ ←	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,375,222	0	A	C	5.3%	236.0 / 2.4	75	G73G (GGT→GGG)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base A (36/35);  new base C (2/2);  total (38/37)
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.47e-01

GTAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAACGGCGGT  >  NC_002947/5375075‑5375367                                                                                                                                                    |                                                                                                                                                  gtAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCCCCCAGGCTTGGGTCCacac                                                                                                                                                  >  8:286470/1‑149 (MQ=255)    gcgcAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCt                                                                                                                                                <  3:334398/148‑1 (MQ=255)    gcgcAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTg                                                                                                                                               >  3:129522/1‑149 (MQ=255)     cgcAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGc                                                                                                                                              <  7:229819/149‑1 (MQ=255)        ccggcGCTGACGCGCTTGGTACCGCCTCGACCACCACCGCGACCAGCGCCTCCGAACAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATg                                                                                                                                           <  7:168626/145‑1 (MQ=255)          gTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCAACAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGc                                                                                                                                         >  2:223612/1‑149 (MQ=255)           tCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAACAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCAGGTGTCGACCCAAGCATgggt                                                                                                                                        <  1:310239/149‑3 (MQ=255)            cGCTGCCGCGCTGGGCACCGCCTCGACCACCACAGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCa                                                                                                                                       >  3:110384/1‑149 (MQ=255)               tGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATAc                                                                                                                                    >  7:229282/1‑149 (MQ=255)                gCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACt                                                                                                                                   >  2:296544/1‑149 (MQ=255)                 ccGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACt                                                                                                                                   <  7:286470/148‑1 (MQ=255)                        gggCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAg                                                                                                                           >  4:179401/1‑149 (MQ=255)                         ggCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAgg                                                                                                                          <  2:229740/149‑1 (MQ=255)                          gCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGc                                                                                                                         >  3:181435/1‑149 (MQ=255)                           cACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGcc                                                                                                                        >  3:332045/1‑149 (MQ=255)                              cGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCa                                                                                                                     <  4:110384/149‑1 (MQ=255)                               gCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCAc                                                                                                                    <  5:336892/149‑1 (MQ=255)                                  tCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCt                                                                                                                 >  5:85522/1‑149 (MQ=255)                                    gACCACCACCGCGACCACCGCCACCGAACAAGGCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCTTGGCCATACTGGTCGAAGGCCGCACGCt                                                                                                                 <  5:91984/147‑1 (MQ=255)                                       caccaccGCGACCACCGCCACCGAAGAAGTCGATGAACACATCTCCGAAGATGTCGGAGAAATTGGCACCAAAAAAGCCGGCACCACCGCCAACCATGCTTGGCTCGACACCAGCATGGCCATACTGGTCGAATGCCGCACGCTAGCTc                                                                                                            >  6:381005/1‑149 (MQ=255)                                        accaccGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCAGGTGTCGACCCAAGCATgggtg                                                                                                                                       >  2:186959/1‑118 (MQ=255)                                         ccaccGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTcgc                                                                                                          >  2:338848/1‑149 (MQ=255)                                            ccGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTc                                                                                                       <  7:137596/149‑1 (MQ=255)                                            ccGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTc                                                                                                       >  7:321443/1‑149 (MQ=255)                                               cgACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGa                                                                                                    <  1:296544/149‑1 (MQ=255)                                               cgACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGa                                                                                                    <  3:240696/149‑1 (MQ=255)                                               cgACCACCGCCACCGAAGAAGTCGCTGAACACACCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGa                                                                                                    >  8:220305/1‑149 (MQ=255)                                                 accaccGACACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCTCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACa                                                                                                  <  5:381005/149‑1 (MQ=255)                                                          accgAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCg                                                                                          >  1:155357/1‑148 (MQ=255)                                                           cctaagaagTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTa                                                                                        <  8:90565/146‑1 (MQ=255)                                                           ccgAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTa                                                                                        <  8:321443/149‑1 (MQ=255)                                                            cgAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAg                                                                                       <  4:129522/149‑1 (MQ=255)                                                             gaagaaGTCGCTGAACACATCGCCGAGGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAgg                                                                                      >  1:202952/1‑149 (MQ=255)                                                                gaagTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGTACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCt                                                                                                      >  3:225728/1‑130 (MQ=255)                                                                gaagTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGTACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCt                                                                                                      <  4:225728/130‑1 (MQ=255)                                                                gaagTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCt                                                                                   <  2:59973/149‑1 (MQ=255)                                                                 aagTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCCCCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTc                                                                                  >  4:157154/1‑149 (MQ=255)                                                                    tCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGtt                                                                               >  2:174193/1‑149 (MQ=255)                                                                        tGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGcc                                                                           >  1:277223/1‑149 (MQ=255)                                                                              caTCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTg                                                                     <  1:261639/149‑1 (MQ=255)                                                                                 cGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAAc                                                                  <  7:220305/149‑1 (MQ=255)                                                                                  gCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAAc                                                                  <  2:155357/148‑1 (MQ=255)                                                                                  gCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAAc                                                                  <  2:202952/148‑1 (MQ=255)                                                                                   ccGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATAc                                                                                                                                    <  1:375573/81‑1 (MQ=255)                                                                                   ccGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATAc                                                                                                                                    >  2:375573/1‑81 (MQ=255)                                                                                   ccGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACtg                                                                >  8:220828/1‑148 (MQ=255)                                                                                      aaGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGAc                                                                                                   >  6:382732/1‑111 (MQ=255)                                                                                      aaGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGAc                                                                                                   <  5:382732/111‑1 (MQ=255)                                                                                         aTGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTc                                                          <  6:85522/149‑1 (MQ=255)                                                                                          tGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTc                                                          <  1:223612/148‑1 (MQ=255)                                                                                           gTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTcc                                                         >  8:141004/1‑148 (MQ=255)                                                                                              ggAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGt                                                              >  4:141548/1‑140 (MQ=255)                                                                                              ggAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGt                                                              <  3:141548/140‑1 (MQ=255)                                                                                              ggAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGAtt                                                     <  6:111126/149‑1 (MQ=255)                                                                                              ggAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGAtt                                                     >  1:87078/1‑149 (MQ=255)                                                                                              ggAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTC‑AAGGCCGCACGCATGCGCGCGGAGGACAGCACTTCGTAAGCGTAGTTGGCCGCCGTGAACTTGTCGTACGATTc                                                    >  8:347502/1‑149 (MQ=255)                                                                                               gagaAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTc                                                    >  4:52273/1‑149 (MQ=255)                                                                                                agaAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCg                                                                                 <  3:340632/119‑1 (MQ=255)                                                                                                agaAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCg                                                                                 >  4:340632/1‑119 (MQ=255)                                                                                                   aaTTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTa                                                <  4:181435/149‑1 (MQ=255)                                                                                                     ttGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATc                                              >  7:219947/1‑149 (MQ=255)                                                                                                          accaccAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGcctc                                                                         >  5:178229/1‑117 (MQ=255)                                                                                                          accaccAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGcctc                                                                         <  6:178229/117‑1 (MQ=255)                                                                                                               caAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTg                                    >  6:158813/1‑149 (MQ=255)                                                                                                               caAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACAGGTCGAAGGCCGCACGCTTGCTGGCGTGTGACAGCACTTCGTAGGCCTCGTTGGCCTCCGTGAACTTGGCTTCCGATTCTTTTTAGCCCGGGTTg                                    >  1:273516/1‑149 (MQ=255)                                                                                                                 aaGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCg                                  >  4:215771/1‑149 (MQ=255)                                                                                                                              gtatggccATGCAGGTGTCGACCCAAGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTATGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGggtg                         >  4:48515/7‑145 (MQ=255)                                                                                                                              gtatggccATGCAGGTGTCGACCCAAGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTATGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGggtg                         <  3:48515/139‑1 (MQ=255)                                                                                                                                  aCCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCa                 <  7:141004/149‑1 (MQ=255)                                                                                                                                    ccATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATc               <  1:174193/149‑1 (MQ=255)                                                                                                                                    ccATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATc               >  3:96385/1‑149 (MQ=255)                                                                                                                                             ggTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAAcgg      <  3:179401/149‑1 (MQ=255)                                                                                                                                              gTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAAcggc     <  3:157154/149‑1 (MQ=255)                                                                                                                                               tCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAAcggcg    >  7:15713/1‑149 (MQ=255)                                                                                                                                               tCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCACCAAAcggcg    <  8:229282/149‑1 (MQ=255)                                                                                                                                                 gACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAAcggcgg   >  2:17325/1‑148 (MQ=255)                                                                                                                                                  acacCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAACGGCGGt  >  4:168134/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GTAGCGCAGGTCGCTGCCGCGCTGGGCACCGCCTCGACCACCACCGCGACCACCGCCACCGAAGAAGTCGCTGAACACATCGCCGAAGATGTCGGAGAAATTGGCACCACCAAAGCCGGCACCACCGCCACCCATGCTTGGGTCGACACCTGCATGGCCATACTGGTCGAAGGCCGCACGCTTGCTCGCGTCTGACAGCACTTCGTAGGCCTCGTTGGCCTCCTTGAACTTGTCTTCCGATTCTTTATCGCCCGGGTTGCGGTCCGGGTGGTACTTCATCGCCAAACGGCGGT  >  NC_002947/5375075‑5375367

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