Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,857,965	C→G	7.6%	R172G (CGG→GGG)	PP_5134 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,857,965	0	C	G	7.6%	238.7 / 9.1	79	R172G (CGG→GGG)	PP_5134	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/37);  new base G (3/3);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

TATGTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAGCGT  >  NC_002947/5857818‑5858111                                                                                                                                                      |                                                                                                                                                   tATGTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGccc                                                                                                                                                    <  8:80479/148‑1 (MQ=255)       ttCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGAGGAACTGGAAGTGCTGGCCTACTGCGGCAAGCCGTTGGTGCCGGTGGTCAACTTCGTCGCCAGCAGCCAG‑GA‑CGGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTgg                                                                                                                                             <  7:224218/148‑1 (MQ=255)           aTGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCAAGCCGGGCAAGGGCtt                                                                                                                                         <  3:193084/149‑1 (MQ=255)            tGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTc                                                                                                                                        <  1:208125/149‑1 (MQ=255)              cgcgcgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCAAGCCGGGCAAGGGCTTCAc                                                                                                                                      <  5:38015/149‑1 (MQ=255)                  cgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCAt                                                                                                                                  >  8:234222/1‑149 (MQ=255)                    aGCCGGGGCTGGCCAAGTACCGCGTCGAACTGGAAGGGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCATCAGCCAG‑CATCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCAtt                                                                                                                                 <  1:290507/149‑1 (MQ=255)                             tGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTc                                                                                                                       >  6:42487/1‑149 (MQ=255)                                   agcaCCGCGACGAACTGGAAGGGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGc                                                                                                                 <  3:278589/146‑1 (MQ=255)                                       ccGCTACGAACTGTAAGTGCTGGCCAAGTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCAAGAAGCATGCCT‑TGGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTg                                                                                                              <  1:314816/149‑1 (MQ=255)                                       ccGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACg                                                                                                                           >  1:283455/1‑135 (MQ=255)                                       ccGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACg                                                                                                                           <  2:283455/135‑1 (MQ=255)                                       ccGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTgg                                                                                                             >  3:371182/1‑149 (MQ=255)                                          cgacgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGccc                                                                                                          >  1:268811/1‑149 (MQ=255)                                           gacgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGcccc                                                                                                         <  7:40519/149‑1 (MQ=255)                                               aaCTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                       >  1:159856/1‑147 (MQ=255)                                                aCTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                       >  5:263531/1‑146 (MQ=255)                                                aCTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                       >  3:301225/1‑146 (MQ=255)                                                 cTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCa                                                                                                     >  7:167800/1‑147 (MQ=255)                                                 cTGGAAGGGCGGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCaga                                                                                                   <  5:35386/149‑1 (MQ=255)                                                  tGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTg                                                                                                                                              >  5:303693/1‑105 (MQ=255)                                                  tGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTg                                                                                                                                              <  6:303693/105‑1 (MQ=255)                                                    gAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGTTACGGTTCGACAGCGTGGCCCCGCCagag                                                                                                  <  2:315593/147‑1 (MQ=255)                                                    gAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCagag                                                                                                  <  4:301225/147‑1 (MQ=255)                                                       gTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGGCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATgg                                                                                             <  7:217863/149‑1 (MQ=255)                                                        tGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                            <  6:263531/149‑1 (MQ=255)                                                          cTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGtt                                                                                                                        >  1:167756/1‑119 (MQ=255)                                                          cTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGtt                                                                                                                        <  2:167756/119‑1 (MQ=255)                                                          cTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                            <  6:271312/147‑1 (MQ=255)                                                             gCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAgcg                                                                                       >  3:207137/1‑149 (MQ=255)                                                             gCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAgcg                                                                                       >  5:300169/1‑149 (MQ=255)                                                               cAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGcc                                                                                     >  6:359144/1‑149 (MQ=255)                                                               cAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGcc                                                                                     >  7:346681/1‑149 (MQ=255)                                                               cAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGcc                                                                                     >  1:29928/1‑149 (MQ=255)                                                                     cGGCATGCCTTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑GCGCGAGCTTGAATGGCGTGAAGCCCTTGCCCGTCTTGGGCTTCACGCTTTGGTACGGTTCGACAGCTTGTCCCCGCCAGATGATGGTGAGCGCGGttt                                                                                 <  8:177826/147‑1 (MQ=255)                                                                     cGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTg                                                                                >  6:93750/1‑148 (MQ=255)                                                                     cGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGt                                                                               <  2:126483/149‑1 (MQ=255)                                                                      ggCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTa                                                                              >  5:357870/1‑149 (MQ=255)                                                                      ggCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTa                                                                              >  3:292953/1‑149 (MQ=255)                                                                       gCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTAc                                                                             >  2:59695/1‑149 (MQ=255)                                                                       gCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTAc                                                                             >  6:183849/1‑149 (MQ=255)                                                                        cAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACg                                                                            <  7:234222/149‑1 (MQ=255)                                                                         aaGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                       <  8:18642/121‑1 (MQ=255)                                                                         aaGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                       >  7:18642/1‑121 (MQ=255)                                                                         aaGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACg                                                                            >  2:338780/1‑148 (MQ=255)                                                                              gTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTg                                                                                           >  3:110307/1‑128 (MQ=255)                                                                              gTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTg                                                                                           <  4:110307/128‑1 (MQ=255)                                                                                  ctgcCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTgg                                                                   <  5:359144/148‑1 (MQ=255)                                                                                  ctgcCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTgg                                                                   <  5:42487/148‑1 (MQ=255)                                                                                  ctgcCGGTGCTCAACTTCGGCGCCACCCGCCAG‑CC‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGgc                                                                  <  8:346681/149‑1 (MQ=255)                                                                                     ccGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgct                                                               >  5:170478/1‑149 (MQ=255)                                                                                     ccGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCAGTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgct                                                               <  8:167800/149‑1 (MQ=255)                                                                                       ggTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgctgc                                                             <  1:338780/149‑1 (MQ=255)                                                                                           cTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgctgctgct                                                         >  5:219214/1‑149 (MQ=255)                                                                                            tCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                            >  7:288021/1‑113 (MQ=255)                                                                                            tCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                            <  8:288021/113‑1 (MQ=255)                                                                                              aaCTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGctgctgctg                                                        <  7:283447/147‑1 (MQ=255)                                                                                              aaCTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGctgctgctg                                                        <  4:207137/147‑1 (MQ=255)                                                                                              aaCTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGACCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGCCAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGGTTGTACGAAAGCCTGCCGCTGCTGCTgg                                                       >  3:69506/1‑148 (MQ=255)                                                                                                    gtcgCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgc                                                >  5:251476/1‑149 (MQ=255)                                                                                                           gcagccag‑ca‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGcccgcccg                                         <  2:159856/149‑1 (MQ=255)                                                                                                             agccag‑ca‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGgc                                       <  8:178708/149‑1 (MQ=255)                                                                                                                     a‑CCGCCAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACCGATCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGCAAAACGCCCGCCCGgcgc                                     >  5:285154/1‑144 (MQ=255)                                                                                                                       ccGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAAc                               <  4:7849/149‑1 (MQ=255)                                                                                                                        cgcgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACg                              >  6:84339/1‑149 (MQ=255)                                                                                                                         gcgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACgg                             >  7:351772/1‑149 (MQ=255)                                                                                                                          cgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGc                            <  5:93750/149‑1 (MQ=255)                                                                                                                          cgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGc                            <  6:219214/149‑1 (MQ=255)                                                                                                                          cgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCAGTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGc                            <  8:351772/149‑1 (MQ=255)                                                                                                                            aGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCt                           >  6:222408/1‑148 (MQ=255)                                                                                                                             gCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGa                         <  6:251476/149‑1 (MQ=255)                                                                                                                                  aaTGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATcga                     <  6:300169/148‑1 (MQ=255)                                                                                                                                       cGTGAAGCCCAAGCCGGGCAAGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCgacgac                 >  2:115421/1‑147 (MQ=255)                                                                                                                                        gTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGAcccg              >  4:71735/1‑147 (MQ=255)                                                                                                                                        gTGAAGCCCAAGCCGGGCAAGGGCTTCACGCATTGGTACGGTTCGATAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACcag              >  3:37989/1‑149 (MQ=255)                                                                                                                                        gTGAAGCCCAAGCCGGGCAAGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCaa                                <  5:223653/131‑1 (MQ=255)                                                                                                                                        gTGAAGCCCAAGCCGGGCAAGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCaa                                >  6:223653/1‑131 (MQ=255)                                                                                                                                              cccTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGca        <  6:302533/149‑1 (MQ=255)                                                                                                                                                cTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGcaca      >  6:259611/1‑149 (MQ=255)                                                                                                                                                 ttGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAg     >  7:330187/1‑149 (MQ=255)                                                                                                                                                    cccGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAGCGt  >  4:78992/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TATGTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAGCGT  >  NC_002947/5857818‑5858111

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