Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,762,930:1	+C	5.6%	coding (173/723 nt)	PP_3326 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,762,930	1	.	C	5.6%	304.6 / 8.8	72	coding (173/723 nt)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (36/32);  new base C (2/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAT  >  NC_002947/3762783‑3763078                                                                                                                                                      |                                                                                                                                                     cccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑c                                                                                                                                                     >  1:36231/1‑149 (MQ=255)  ccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑c                                                                                                                                                     <  7:285701/148‑1 (MQ=255)      ttGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGa                                                                                                                                                 <  1:61381/148‑1 (MQ=255)        ggCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATAGACGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGAcc                                                                                                                                               >  1:97782/1‑149 (MQ=255)           ccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCggg                                                                                                                                            >  6:187352/1‑148 (MQ=255)           ccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCggg                                                                                                                                            >  5:91950/1‑148 (MQ=255)           ccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCggg                                                                                                                                            >  3:219814/1‑148 (MQ=255)           ccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCggg                                                                                                                                            >  8:133433/1‑148 (MQ=255)              gATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACa                                                                                                                                        >  7:152147/1‑149 (MQ=255)              gATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGAGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACa                                                                                                                                        >  2:68392/1‑149 (MQ=255)               aTGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACa                                                                                                                                        <  3:310570/148‑1 (MQ=255)                 ggTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGAACGGGAACATg                                                                                                                                      >  4:261361/1‑148 (MQ=255)                    cAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATcc                                                                                                                                  >  5:99912/1‑149 (MQ=255)                      gCGTCAACGCAGGCGCCCACGGCCTGGGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCg                                                                                                                                <  3:110085/149‑1 (MQ=255)                         tCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTc                                                                                                                             <  1:48068/149‑1 (MQ=255)                          cAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCa                                                                                                                            >  7:256585/1‑149 (MQ=255)                                 ggCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGa                                                                                                                     <  1:86102/149‑1 (MQ=255)                                      ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGcc                                                                                                                >  6:58730/1‑149 (MQ=255)                                        aCGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTg                                                                                                              >  7:55645/1‑149 (MQ=255)                                        aCGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTg                                                                                                              <  2:97782/149‑1 (MQ=255)                                           gCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCg                                                                                                           <  8:152147/149‑1 (MQ=255)                                               gtgGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGAt                                                                                                        >  3:120012/1‑148 (MQ=255)                                                  gACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGc                                                                                                            >  3:51723/1‑141 (MQ=255)                                                  gACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGc                                                                                                            <  4:51723/141‑1 (MQ=255)                                                       aggagCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCgg                                                                                                          <  2:239242/138‑1 (MQ=255)                                                       aggagCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGAAAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCgg                                                                                                          >  1:239242/1‑138 (MQ=255)                                                                  ggCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACt                                                                                    >  4:265017/1‑149 (MQ=255)                                                                     aaTA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAc                                                                                 >  3:19690/1‑149 (MQ=255)                                                                       ta‑tCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAccc                                                                               >  4:6633/1‑149 (MQ=255)                                                                        a‑tCGAGGCCGTGTACGGCCACGGCGCCGAGGAAAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGATCGA‑CGTGACCGGGAACATGATCCCGGTCACCGTCGAGCGCCTGGCGGATCATGCCGAGCTGCAACCACATAcccc                                                                              >  6:68318/1‑149 (MQ=255)                                                                          tCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAcccc                                                                              <  4:42274/148‑1 (MQ=255)                                                                          tCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCg                                                                             >  4:292153/1‑149 (MQ=255)                                                                             tggCCGTACACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCt                                                                          >  1:66296/2‑149 (MQ=255)                                                                              ggCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTg                                                                         >  7:118191/1‑149 (MQ=255)                                                                                 cgtacaCGGCCTCGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCgg                                                                      <  1:210354/144‑1 (MQ=255)                                                                                   tgtACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCa                                                                    <  8:256585/149‑1 (MQ=255)                                                                                       cGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCAC‑GTCGAacgc                                                                                                                 >  2:40558/1‑100 (MQ=21)                                                                                       cGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCAC‑GTCGAacgc                                                                                                                 <  1:40558/100‑1 (MQ=21)                                                                                          ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCa                                                                                      >  3:66417/1‑124 (MQ=255)                                                                                          ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCa                                                                                      <  4:66417/124‑1 (MQ=255)                                                                                           cACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTc                                                                                                                       >  1:176574/1‑90 (MQ=255)                                                                                           cACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTc                                                                                                                       <  2:176574/90‑1 (MQ=255)                                                                                           cACGGCCTCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                            <  8:245478/149‑1 (MQ=255)                                                                                                 gcCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTgctgct                                                      <  2:36231/149‑1 (MQ=255)                                                                                                  ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGctgctg                                                     <  8:118191/149‑1 (MQ=255)                                                                                                  ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCC‑AGGCGCTCGACCGTGACCGGGATCATGTTcc                                                                                                                                  >  6:69232/1‑72 (MQ=17)                                                                                                  ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCC‑AGGCGCTCGACCGTGACCGGGATCATGTTcc                                                                                                                                  <  5:69232/72‑1 (MQ=17)                                                                                                    gAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGa                                                   >  4:213135/1‑149 (MQ=255)                                                                                                    gAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGa                                                   <  2:263346/149‑1 (MQ=255)                                                                                                                     ggCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTca                                   >  3:179089/1‑148 (MQ=255)                                                                                                                      gCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                 <  1:162839/149‑1 (MQ=255)                                                                                                                        ggAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGt                                     >  5:235248/1‑143 (MQ=255)                                                                                                                        ggAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGt                                     <  6:235248/143‑1 (MQ=255)                                                                                                                         gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTAACGGTCGAGCTCCTGGCGGATCATGCCGTGCTGCAACCACTTACACCGCCTGCGGCACTTTCCGTGCTGATGGACAACGGGCCTTGGTAACAgcg                              >  6:164198/1‑149 (MQ=255)                                                                                                                            aaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGccc                           >  6:242881/1‑149 (MQ=255)                                                                                                                               aTCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGa                        >  5:12811/1‑149 (MQ=255)                                                                                                                                  aaGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACaa                     <  3:6633/149‑1 (MQ=255)                                                                                                                                  aaGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACaa                     <  4:219814/149‑1 (MQ=255)                                                                                                                                   aGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAg                    >  4:38417/1‑149 (MQ=255)                                                                                                                                       cgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAt                <  3:265017/149‑1 (MQ=255)                                                                                                                                         ccTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg              <  7:133433/149‑1 (MQ=255)                                                                                                                                         ccTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg              <  4:19690/149‑1 (MQ=255)                                                                                                                                          cTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCg                          <  8:196915/136‑1 (MQ=255)                                                                                                                                          cTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCg                          >  7:196915/1‑136 (MQ=255)                                                                                                                                          cTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATggg             >  3:30475/1‑149 (MQ=255)                                                                                                                                          cTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg              >  8:35713/1‑148 (MQ=255)                                                                                                                                          cTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg              <  7:35713/148‑1 (MQ=255)                                                                                                                                           tgtgCGTTCGT‑CGTGACCTGTAACTTGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCGTGCGGCACTTGCCGTGCTGCTGGACATCGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt            <  4:308793/146‑1 (MQ=255)                                                                                                                                             ggCGTTCGA‑CGTGACCGGGAAAATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCAAGCCGTGCAGCAACCACTTACCCCGCCTGCGGAACTTGCCCTGCTGCTGGACAACGCGCCTTGGTCACAGCGCCCGGACAAGCAATCGGt            >  3:308793/1‑147 (MQ=255)                                                                                                                                                 ttCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCa      <  3:38417/149‑1 (MQ=255)                                                                                                                                                   cGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAg     <  6:91950/148‑1 (MQ=255)                                                                                                                                                    gA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCTCTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAg                                <  1:218138/120‑1 (MQ=255)                                                                                                                                                    gA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAg                                >  2:218138/1‑120 (MQ=255)                                                                                                                                                    gA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGa   <  5:58730/149‑1 (MQ=255)                                                                                                                                                     a‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATAATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt  >  5:247649/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATA‑TCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAT  >  NC_002947/3762783‑3763078

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