Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,433,112	Δ1 bp	5.9%	coding (333/351 nt)	PP_5648 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,433,112	0	A	.	5.9%	216.6 / 7.0	68	coding (333/351 nt)	PP_5648	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base . (2/2);  total (34/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 = 9.55e-01

ACCGCGCGCGCTGCGACACAACTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGG  >  NC_002947/4432964‑4433256                                                                                                                                                     |                                                                                                                                                  aCCGCGCGCGCTGCGACACAACTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTgcagca                                                                                                                                                   <  3:152203/149‑1 (MQ=255)        gcgcTGCGACACAACTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTc                                                                                                                                            <  2:235939/149‑1 (MQ=255)          gcTGCGACACAACTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCgg                                                                                                                                          >  6:351676/1‑149 (MQ=255)                acacAACTGCCTTGAAAAAACCTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGc                                                                                                                                   <  3:87115/149‑1 (MQ=255)                    aaCTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCagag                                                                                                                               >  1:357818/1‑149 (MQ=255)                      cTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGc                                                                                                                              >  4:208105/1‑148 (MQ=255)                      cTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGc                                                                                                                              >  3:50423/1‑148 (MQ=255)                           ttGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGcc                                                                                                                        >  6:22962/1‑149 (MQ=255)                              aaaaaaaCTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCcgtc                                                                                                                     >  7:58576/1‑149 (MQ=255)                               aaaaaaCTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCcgtcg                                                                                                                    <  5:45125/149‑1 (MQ=255)                                aaaaaCTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGC‑GCCGAGCGGCTGCTGCaac                                                                                                                                >  5:64001/1‑134 (MQ=255)                                aaaaaCTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGC‑GCCGAGCGGCTGCTGCaac                                                                                                                                <  6:64001/136‑3 (MQ=255)                                  aaaCTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTc                                                                                                                 <  2:304587/149‑1 (MQ=255)                                    aCTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGC‑GCCGAGCGGCTGCTGCa                                                                                                                                  <  3:33072/130‑1 (MQ=255)                                    aCTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGC‑GCCGAGCGGCTGCTGCa                                                                                                                                  >  4:33072/1‑130 (MQ=255)                                        tAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTc                                                                                                           >  3:11279/1‑149 (MQ=255)                                        tAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTc                                                                                                           >  7:7387/1‑149 (MQ=255)                                          aaGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTcgc                                                                                                         <  3:151211/149‑1 (MQ=255)                                           aGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTcgcg                                                                                                        <  3:131132/149‑1 (MQ=255)                                            gCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCACTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTcgcg                                                                                                        >  3:183781/1‑148 (MQ=255)                                                  gcgcCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAAt                                                                                                 <  3:208105/149‑1 (MQ=255)                                                   cgcCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACGCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATc                                                                                                <  8:7387/149‑1 (MQ=255)                                                   cgcCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATc                                                                                                <  8:58576/149‑1 (MQ=255)                                                            tgtgGCAACAACTGCCTTGCTCAGTACACTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGt                                                                                       >  5:24966/1‑149 (MQ=255)                                                              tgGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCTGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTc                                                                                      <  4:50423/148‑1 (MQ=255)                                                              tgGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTcg                                                                                                          <  8:152954/128‑1 (MQ=255)                                                              tgGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTcg                                                                                                          >  7:152954/1‑128 (MQ=255)                                                              tgGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTc                                                                                      <  5:22962/148‑1 (MQ=255)                                                                   acaacTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTcccgccc                                                                                <  2:357818/149‑1 (MQ=255)                                                                     aacTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGa                                                                              >  6:211223/1‑149 (MQ=255)                                                                        tGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATc                                                                                                >  2:135775/1‑128 (MQ=255)                                                                        tGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTg                                                                           >  7:208724/1‑149 (MQ=255)                                                                          gcTTTCTCAGTACGCTGAATGCACGTGTCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTc                                                                                      <  4:33592/135‑1 (MQ=255)                                                                          ccggccTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGc                                                                          <  2:277129/143‑1 (MQ=255)                                                                          ccTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCgtag                                                                                                     <  8:296497/121‑1 (MQ=255)                                                                          ccTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCgtag                                                                                                     >  7:296497/1‑121 (MQ=255)                                                                          ccTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTc                                                                                      >  3:33592/1‑136 (MQ=255)                                                                          ccTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCg                                                                         >  8:81929/1‑149 (MQ=255)                                                                          ccTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGAAATCCACACGCGTCCCGCCCGACTGCg                                                                         >  8:86918/1‑149 (MQ=255)                                                                           cTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCgg                                                                        <  4:11279/149‑1 (MQ=255)                                                                           cTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCgg                                                                        <  5:27146/149‑1 (MQ=255)                                                                             tGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCggg                                                                       >  2:127016/1‑148 (MQ=255)                                                                               cTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTaa                                                                    <  6:24966/149‑1 (MQ=255)                                                                                       ccTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCAc                                                            >  7:240075/1‑149 (MQ=255)                                                                                        cTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCAc                                                            <  1:127016/148‑1 (MQ=255)                                                                                           aaTGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAAc                                                        <  5:351676/149‑1 (MQ=255)                                                                                           aaTGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCACTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGaa                                                         <  4:183781/148‑1 (MQ=255)                                                                                                  gggCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGCCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAAc                                                        <  2:161897/142‑1 (MQ=255)                                                                                                  gggCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGCCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAAc                                                        >  1:161897/1‑142 (MQ=255)                                                                                                           cAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCggg                                                                       <  5:341153/118‑1 (MQ=255)                                                                                                           cAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCggg                                                                       >  6:341153/1‑118 (MQ=255)                                                                                                              gAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTcccgc                                                                                  >  7:69493/1‑104 (MQ=255)                                                                                                                  ggcgATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAacac                                 >  6:15989/1‑149 (MQ=255)                                                                                                                    cgATAGCGCCAGCCGTTACAGGCTGGTTGCAGTAGGCGCTCTTC‑TCTGCAGAGCCACGCCGTCGTTATGCTTCGAGTATTAATCCACTCGGGTCCGGCCCGTCTGAGGGTAAACGTCCAGTAAAGCTTCGGCCTCATGGGTATacaagt                               <  8:208724/149‑4 (MQ=255)                                                                                                                     gATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGt                              <  5:15989/149‑1 (MQ=255)                                                                                                                      aTAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGt                              >  2:248280/1‑148 (MQ=255)                                                                                                                         gcgcCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAAt                          <  5:211223/149‑1 (MQ=255)                                                                                                                            ccAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACgg                                         >  6:206394/1‑131 (MQ=255)                                                                                                                            ccAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACgg                                         <  5:206394/131‑1 (MQ=255)                                                                                                                            ccAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCg                       >  6:199892/1‑149 (MQ=255)                                                                                                                             cAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGc                      >  8:73660/1‑149 (MQ=255)                                                                                                                               gCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGt                              <  7:300489/139‑1 (MQ=255)                                                                                                                               gCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGt                              >  8:300489/1‑139 (MQ=255)                                                                                                                                       aggcaggTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTAgcg            >  7:82017/1‑149 (MQ=255)                                                                                                                                        ggcaggTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTgcgc           <  5:199892/149‑1 (MQ=255)                                                                                                                                            ggTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCgg                                                                        <  2:123990/84‑1 (MQ=255)                                                                                                                                            ggTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCgg                                                                        >  1:123990/1‑84 (MQ=255)                                                                                                                                             gTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGc                      >  3:280649/1‑133 (MQ=255)                                                                                                                                             gTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGc                      <  4:280649/133‑1 (MQ=255)                                                                                                                                                gcagcaGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAg   >  3:327460/1‑149 (MQ=255)                                                                                                                                                 cagcagCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTAGCGCGGCTCCAgg  <  8:82017/149‑1 (MQ=255)                                                                                                                                                  agcagcCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAgg  >  5:160787/1‑148 (MQ=255)                                                                                                                                                   gcagcCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAgg  >  1:176494/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                  ACCGCGCGCGCTGCGACACAACTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGG  >  NC_002947/4432964‑4433256

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