Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,433,122:1	+T	5.3%	coding (323/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,122	1	.	T	5.3%	304.1 / 11.0	75	coding (323/351 nt)	PP_5648	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (36/35);  new base T (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 = 7.98e-01

ACACAACTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGG  >  NC_002947/4432979‑4433266                                                                                                                                                 |                                                                                                                                                 acacAACTGCCTTGAAAAAACCTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGc                                                                                                                                             <  3:87115/149‑1 (MQ=255)     aaCTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCagag                                                                                                                                         >  1:357818/1‑149 (MQ=255)       cTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGc                                                                                                                                        >  3:50423/1‑148 (MQ=255)       cTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGc                                                                                                                                        >  4:208105/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                                                                                                                                            >  4:33072/1‑130 (MQ=255)                     aCTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGC‑GCCGAGCGGCTGCTGCa                                                                                                                                            <  3:33072/130‑1 (MQ=255)                         tAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTc                                                                                                                     >  7:7387/1‑149 (MQ=255)                         tAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTc                                                                                                                     >  3:11279/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                                                                                                                    >  7:152954/1‑128 (MQ=255)                                               tgGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTcg                                                                                                                    <  8:152954/128‑1 (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                                                                                                               >  7:296497/1‑121 (MQ=255)                                                           ccTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCgtag                                                                                                               <  8:296497/121‑1 (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                                                                                  <  5:27146/149‑1 (MQ=255)                                                            cTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCgg                                                                                  <  4:11279/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                                        >  8:300489/1‑139 (MQ=255)                                                                                                                gCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGt                                        <  7:300489/139‑1 (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)                                                                                                                                           cTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAg             >  7:330741/1‑139 (MQ=255)                                                                                                                                           cTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAg             <  8:330741/139‑1 (MQ=255)                                                                                                                                             cGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTgcgc                     >  3:317841/1‑129 (MQ=255)                                                                                                                                             cGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTgcgc                     <  4:317841/129‑1 (MQ=255)                                                                                                                                              ggC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTgg  <  1:248280/147‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 ACACAACTGCCTTGAAAAAAACTGTAAAGCCCGCGCGCCTCTGCTGTGGCAACAACTGCCTTGCTCAGTACCCTGAATGCACGGGCCCTTACAGGAGCGGCGATAGCGCCAGCCGTTACAGGCAGGTTGCAGCAGCCGCTCGGC‑GCTGCAGAGCCACGCCGTCGTTCTGCTTCGCGTAGTAATCCACACGCGTCCCGCCCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGG  >  NC_002947/4432979‑4433266

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