Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,857,969	T→A	7.5%	L173H (CTT→CAT) ‡	PP_5134 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,857,969	0	T	A	7.5%	233.3 / 7.9	80	L173H (CTT→CAT) ‡	PP_5134	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base A (3/3);  total (40/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 = 2.28e-01

ATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAGCGT  >  NC_002947/5857824‑5858111                                                                                                                                                    |                                                                                                                                               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                                                                                                       >  3:301225/1‑146 (MQ=255)                                          aCTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                       >  5:263531/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                                                                                                                                              <  6:303693/105‑1 (MQ=255)                                            tGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTg                                                                                                                                              >  5:303693/1‑105 (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                                                                                     >  1:29928/1‑149 (MQ=255)                                                         cAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGcc                                                                                     >  7:346681/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                                                                             >  6:183849/1‑149 (MQ=255)                                                                 gCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTAc                                                                             >  2:59695/1‑149 (MQ=255)                                                                  cAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACg                                                                            <  7:234222/149‑1 (MQ=255)                                                                   aaGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                       >  7:18642/1‑121 (MQ=255)                                                                   aaGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                       <  8:18642/121‑1 (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:42487/148‑1 (MQ=255)                                                                            ctgcCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTgg                                                                   <  5:359144/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                            <  6:219214/149‑1 (MQ=255)                                                                                                                    cgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGc                            <  5:93750/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)                                                                                                                                                    |                                                                                                                                               ATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAG‑CA‑CCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAGCGT  >  NC_002947/5857824‑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: ‑