Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,645,319	A→G	6.7%	A224A (GCA→GCG)	PP_4954 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,645,319	0	A	G	6.7%	216.7 / 4.7	74	A224A (GCA→GCG)	PP_4954	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (33/36);  new base G (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 = 9.52e-01

CAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGTTTT  >  NC_002947/5645173‑5645467                                                                                                                                                   |                                                                                                                                                     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                                                                                                                                       <  8:118208/149‑1 (MQ=255)              gACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGatcatc                                                                                                                                       <  6:224426/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                                                                                                                            >  2:100956/1‑144 (MQ=255)                              gggCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGt                                                                                                                            <  1:100956/144‑1 (MQ=255)                               ggCAGCTAACTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCa                                                                                                                       >  7:352405/1‑148 (MQ=255)                                             gCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACt                                                                                                        >  7:297746/1‑149 (MQ=255)                                              cTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                                       <  7:322262/149‑1 (MQ=255)                                               tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACa                                                                                                                    >  4:339276/1‑135 (MQ=255)                                               tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACa                                                                                                                    <  3:339276/135‑1 (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                                                                                                    <  8:238375/149‑1 (MQ=255)                                                 aaGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGaa                                                                                                    <  7:48165/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                                                                             <  1:183588/148‑1 (MQ=255)                                                                         cGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                             <  5:54296/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:305738/149‑1 (MQ=255)                                                                                                                            aaGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTa                         <  5:300750/149‑1 (MQ=255)                                                                                                                                 ggtggCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCg                                    >  4:227598/1‑133 (MQ=255)                                                                                                                                 ggtggCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCg                                    <  3:227598/133‑1 (MQ=255)                                                                                                                                 cttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAg                                                                                                   >  6:10551/3‑70 (MQ=25)                                                                                                                                 cttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAg                                                                                                   <  5:10551/68‑1 (MQ=25)                                                                                                                                  ttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg                    >  8:316882/2‑148 (MQ=255)                                                                                                                                  ttggCTGCGTTCAGCGCGGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg                    >  3:20300/2‑148 (MQ=255)                                                                                                                                       cgcgCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCAACGAGCGTGTCAGTGCCTCGCACTCCCCCTACAACGCACTGa              >  2:322391/1‑149 (MQ=255)                                                                                                                                                 gCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGtt    >  1:285607/1‑149 (MQ=255)                                                                                                                                                   aGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGtttt  >  1:92848/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGACCGAGGCGTTTT  >  NC_002947/5645173‑5645467

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