Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,645,307	C→T	6.8%	A220A (GCC→GCT)	PP_4954 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,645,307	0	C	T	6.8%	213.5 / 3.9	73	A220A (GCC→GCT)	PP_4954	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base T (3/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-02

ACCAGCCGGTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGA  >  NC_002947/5645159‑5645455                                                                                                                                                     |                                                                                                                                                     aCCAGCCGGTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGcc                                                                                                                                                      >  7:301492/1‑149 (MQ=255)  ccAGCCGGTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCcg                                                                                                                                                     >  7:118208/1‑149 (MQ=255)  ccAGCCGGTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCcg                                                                                                                                                     >  7:238375/1‑149 (MQ=255)    aGCCGGTACAGCAGAGCGTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCcgcg                                                                                                                                                   <  7:291745/149‑1 (MQ=255)         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                                                                                                                           <  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                                                                                                                <  1:100956/144‑1 (MQ=255)                                            gggCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGt                                                                                                                >  2:100956/1‑144 (MQ=255)                                             ggCAGCTAACTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCa                                                                                                           >  7:352405/1‑148 (MQ=255)                                                     ccTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGcc                                                                                                                                                      >  4:60123/1‑97 (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                                                                                              <  3:346093/143‑1 (MQ=255)                                                               aaGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCaaa                                                                                              >  4:346093/1‑143 (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                                                                                 <  6:114911/149‑1 (MQ=255)                                                                      gCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGc                                                                                 >  4:152476/1‑149 (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:300750/1‑148 (MQ=255)                                                                               aaCCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCt                                                                         >  6:305738/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                                                           <  3:287449/149‑1 (MQ=255)                                                                                            gCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTcc                                                           <  7:124941/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                                                          <  8:330125/138‑1 (MQ=255)                                                                                                        ggCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCt                                                          >  7:330125/1‑138 (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)                                                                                                                                                     |                                                                                                                                                     ACCAGCCGGTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGA  >  NC_002947/5645159‑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: ‑