Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,645,313	G→C	6.7%	L222L (CTG→CTC) ‡	PP_4954 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,645,313	0	G	C	6.7%	238.4 / 6.3	74	L222L (CTG→CTC) ‡	PP_4954	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base C (3/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

GTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGA  >  NC_002947/5645167‑5645455                                                                                                                                                   |                                                                                                                                               gTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGaa                                                                                                                                              >  5:15038/1‑149 (MQ=255)       cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGc                                                                                                                                        >  8:1200/1‑149 (MQ=255)          agCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCTGCGTTCAGCGCGGCCac                                                                                                                                     <  5:104965/149‑2 (MQ=255)             cTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGga                                                                                                                                  <  2:76218/149‑1 (MQ=255)                gCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGTTGGCCGCGCTGAACGCAGCCAAg                                                                                                                                    >  3:75466/1‑144 (MQ=255)                    gACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGatcatc                                                                                                                           <  6:224426/149‑1 (MQ=255)                    gACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGatcatc                                                                                                                           <  8:118208/149‑1 (MQ=255)                      cTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGATCGCAGCCAAGGAGATCATcc                                                                                                                          <  8:288887/148‑1 (MQ=255)                            ttcGAGCTGGGCAGCTACCTGAAGCTCATGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCg                                                                                                                   >  1:213814/1‑149 (MQ=255)                            ttcGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCg                                                                                                                   >  8:48165/1‑149 (MQ=255)                             tcGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCgg                                                                                                                  >  4:175904/1‑149 (MQ=255)                               gAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGt                                                                                                                <  5:216607/149‑1 (MQ=255)                                 gCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTgg                                                                                                              <  7:292633/149‑1 (MQ=255)                                 gCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTgg                                                                                                              >  2:183588/1‑149 (MQ=255)                                  cTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTggg                                                                                                             <  2:367932/149‑1 (MQ=255)                                    gggCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGt                                                                                                                <  1:100956/144‑1 (MQ=255)                                    gggCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGt                                                                                                                >  2:100956/1‑144 (MQ=255)                                     ggCAGCTAACTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCa                                                                                                           >  7:352405/1‑148 (MQ=255)                                                   gCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACt                                                                                            >  7:297746/1‑149 (MQ=255)                                                    cTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                           <  7:322262/149‑1 (MQ=255)                                                     tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACa                                                                                                        <  3:339276/135‑1 (MQ=255)                                                     tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACa                                                                                                        >  4:339276/1‑135 (MQ=255)                                                      cAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCg                                                                                          >  7:185356/1‑148 (MQ=255)                                                       aaGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCaaa                                                                                              >  4:346093/1‑143 (MQ=255)                                                       aaGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCaaa                                                                                              <  3:346093/143‑1 (MQ=255)                                                       aaGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                           >  8:322262/1‑146 (MQ=255)                                                       aaGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGaa                                                                                        <  7:48165/149‑1 (MQ=255)                                                       aaGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGaa                                                                                        <  8:238375/149‑1 (MQ=255)                                                         ggCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAgg                                                                                      >  6:226151/1‑149 (MQ=255)                                                          gCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAgg                                                                                      >  7:23950/1‑148 (MQ=255)                                                           ccAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGt                                                                                     >  1:343987/1‑148 (MQ=255)                                                              gCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGc                                                                                 >  4:152476/1‑149 (MQ=255)                                                              gCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGc                                                                                 <  6:114911/149‑1 (MQ=255)                                                                   tttGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTa                                                                            >  3:331533/1‑149 (MQ=255)                                                                    ttGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTAc                                                                           <  8:301492/149‑1 (MQ=255)                                                                       aaCCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCt                                                                         >  2:135930/1‑148 (MQ=255)                                                                       aaCCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCt                                                                         >  6:305738/1‑148 (MQ=255)                                                                       aaCCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCt                                                                         >  6:300750/1‑148 (MQ=255)                                                                               cGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                 <  5:54296/148‑1 (MQ=255)                                                                               cGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                 <  1:183588/148‑1 (MQ=255)                                                                               cGCACGCTCGACGTCGAGGCACAGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTa                                                                >  5:147163/1‑149 (MQ=255)                                                                                 cACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACt                                                              >  3:140652/1‑149 (MQ=255)                                                                                    gCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTcc                                                           >  6:267571/1‑149 (MQ=255)                                                                                    gCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTcc                                                           <  7:124941/149‑1 (MQ=255)                                                                                    gCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTcc                                                           <  3:287449/149‑1 (MQ=255)                                                                                     cacGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCt                                                          <  8:310602/147‑1 (MQ=255)                                                                                       cGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTgg                                                        <  6:15038/149‑1 (MQ=255)                                                                                          cGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGccc                                                     <  8:23950/149‑1 (MQ=255)                                                                                           gTCGAGGCACAGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCa                                                    <  6:147163/149‑1 (MQ=255)                                                                                              gAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAgg                                                  <  7:1200/148‑1 (MQ=255)                                                                                                ggCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCt                                                          >  7:330125/1‑138 (MQ=255)                                                                                                ggCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCt                                                          <  8:330125/138‑1 (MQ=255)                                                                                                 gCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGatat                                              <  8:185356/149‑1 (MQ=255)                                                                                                      gcgcCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACg                                         <  3:175904/149‑1 (MQ=255)                                                                                                       cgcCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGa                                        <  8:352405/149‑1 (MQ=255)                                                                                                         ccTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGc                                      <  4:140652/149‑1 (MQ=255)                                                                                                               aCTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCa                                <  1:135930/149‑1 (MQ=255)                                                                                                                    cccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGcc                           >  7:102463/1‑149 (MQ=255)                                                                                                                     ccAGCAGAACGGCAAGGTGGTGGCCTCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCt                          <  5:267571/149‑1 (MQ=255)                                                                                                                     ccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCt                          >  4:96902/1‑149 (MQ=255)                                                                                                                      cagcagAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTc                         >  4:232388/1‑149 (MQ=255)                                                                                                                        gcagAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGc                       <  8:297746/149‑1 (MQ=255)                                                                                                                        gcagAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCACAGGATATCCACGAGCGTGTCAGTGCCTCg                        >  8:53658/1‑148 (MQ=255)                                                                                                                         cagAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCGGGATATCCACGAGCGTGTCAGTGCCTCGCa                      <  2:343987/149‑1 (MQ=255)                                                                                                                          agAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCAc                     >  7:307246/1‑149 (MQ=255)                                                                                                                                  aaGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTa             <  5:300750/149‑1 (MQ=255)                                                                                                                                  aaGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTa             <  5:305738/149‑1 (MQ=255)                                                                                                                                       ggtggCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCg                        >  4:227598/1‑133 (MQ=255)                                                                                                                                       ggtggCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCg                        <  3:227598/133‑1 (MQ=255)                                                                                                                                       cttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAg                                                                                       <  5:10551/68‑1 (MQ=25)                                                                                                                                       cttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAg                                                                                       >  6:10551/3‑70 (MQ=25)                                                                                                                                        ttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg        >  8:316882/2‑148 (MQ=255)                                                                                                                                        ttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg        >  3:20300/2‑148 (MQ=255)                                                                                                                                             cgcgCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCAACGAGCGTGTCAGTGCCTCGCACTCCCCCTACAACGCACTGa  >  2:322391/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                               GTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGA  >  NC_002947/5645167‑5645455

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