Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,857,961	T→A	5.9%	L170L (CTT→CTA)	PP_5134 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,857,961	0	T	A	5.9%	208.7 / 2.8	68	L170L (CTT→CTA)	PP_5134	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/36);  new base A (2/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CTGTATGTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAG  >  NC_002947/5857815‑5858108                                                                                                                                                   |                                                                                                                                                    cTGTATGTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGc                                                                                                                                                   >  8:19187/1‑149 (MQ=255)      tgtgATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCgg                                                                                                                                               <  2:86964/148‑1 (MQ=255)      tgtgATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCgg                                                                                                                                               <  5:313403/148‑1 (MQ=255)        tgATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCtt                                                                                                                                            >  8:172236/1‑149 (MQ=255)          aTCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTgg                                                                                                                                          <  1:246677/149‑1 (MQ=255)          aTCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTg                                                                                                                                           <  7:100443/148‑1 (MQ=255)          aTCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTg                                                                                                                                           <  1:64027/148‑1 (MQ=255)           tCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTggg                                                                                                                                         >  6:172058/1‑149 (MQ=255)             gATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGc                                                                                                                                        <  8:119812/148‑1 (MQ=255)                 cgcgcgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCAAGCCGGGCAAGGGCTTCAc                                                                                                                                   <  4:279161/149‑1 (MQ=255)                      gAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCAtt                                                                                                                              >  5:194945/1‑149 (MQ=255)                       aGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTg                                                                                                                             <  5:277296/149‑1 (MQ=255)                          cGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTa                                                                                                                          <  8:248094/149‑1 (MQ=255)                          cGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTa                                                                                                                          <  1:112404/149‑1 (MQ=255)                           ggTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTAc                                                                                                                         <  5:126175/149‑1 (MQ=255)                                          ccGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCAc                                                                                                                                   <  6:283487/124‑1 (MQ=255)                                          ccGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCAc                                                                                                                                   >  5:283487/1‑124 (MQ=255)                                          ccGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTgg                                                                                                          >  4:11028/1‑149 (MQ=255)                                           cgcgACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGc                                                                                                         >  1:58650/1‑149 (MQ=255)                                              gacgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGcccc                                                                                                      >  6:10153/1‑149 (MQ=255)                                              gacgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGcccc                                                                                                      <  6:194945/149‑1 (MQ=255)                                                cgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTTGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                    <  7:172236/149‑1 (MQ=255)                                                cgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                    <  5:172058/149‑1 (MQ=255)                                                cgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGc                                                                                                    <  6:36837/149‑1 (MQ=255)                                                 gaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGa                                                                                                                  >  6:215938/1‑134 (MQ=255)                                                 gaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGa                                                                                                                  <  5:215938/134‑1 (MQ=255)                                                 gaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGTCAGCGTGGCCCCGcc                                                                                                   <  8:89049/149‑1 (MQ=255)                                                       gAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCagag                                                                                               <  3:249099/147‑1 (MQ=255)                                                       gAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGa                                                                                             >  1:84424/1‑149 (MQ=255)                                                           tGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                         <  2:84424/149‑1 (MQ=255)                                                             cTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCtt                                                                                                                                            <  2:104986/96‑1 (MQ=255)                                                             cTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGa                                                                                       <  7:19187/149‑1 (MQ=255)                                                                gCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAg                                                                                      <  2:58650/147‑1 (MQ=255)                                                                     cTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGtt                                                                               >  1:256726/1‑149 (MQ=255)                                                                          gCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTAc                                                                          <  2:256726/149‑1 (MQ=255)                                                                              gCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGaaa                                                                      <  6:316644/149‑1 (MQ=255)                                                                              gCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGaaa                                                                      <  5:10153/149‑1 (MQ=255)                                                                                cGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGc                                                                    >  6:279426/1‑149 (MQ=255)                                                                                       gcCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGgcgc                                                             >  1:191398/1‑149 (MQ=255)                                                                                         cGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgctg                                                           <  3:11028/149‑1 (MQ=255)                                                                                               tCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGTCAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGctgctgctg                                                     >  7:205195/1‑149 (MQ=255)                                                                                                cAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGgc                                                               >  7:96659/1‑138 (MQ=255)                                                                                                cAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGgc                                                               <  8:96659/138‑1 (MQ=255)                                                                                                 aaCTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTa                                                                           >  5:275274/1‑125 (MQ=255)                                                                                                 aaCTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTa                                                                           <  6:275274/125‑1 (MQ=255)                                                                                                 aaCTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGctgctgctg                                                     <  8:270866/147‑1 (MQ=255)                                                                                                 aaCTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTgg                                                    <  7:19655/148‑1 (MQ=255)                                                                                                     tcgtcgCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgctg                                                           >  1:153950/1‑137 (MQ=255)                                                                                                     tcgtcgCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgctg                                                           <  2:153950/137‑1 (MQ=255)                                                                                                        tcgCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCAGGCGCTGCTGCTGGAAAAcgcc                                            >  1:278652/1‑149 (MQ=255)                                                                                                           ccagcagccagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgcccgc                                         >  7:195710/1‑149 (MQ=255)                                                                                                                   cagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCt                                 >  3:24704/1‑149 (MQ=255)                                                                                                                     gcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTg                                >  8:44951/1‑148 (MQ=255)                                                                                                                     gcaCCGCGAGCCTGAATGGCGTGAAACCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgc                                             >  3:229000/1‑135 (MQ=255)                                                                                                                     gcaCCGCGAGCCTGAATGGCGTGAAACCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgc                                             <  4:229000/135‑1 (MQ=255)                                                                                                                         cgcgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACg                           >  1:69308/1‑149 (MQ=255)                                                                                                                           cgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGc                         <  2:191398/149‑1 (MQ=255)                                                                                                                               ccTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAACGCCTGGCGCTGCTGCTGGAAACCGCCCGCCCGGCGCTGCAACGGCTGag                     >  8:233593/1‑148 (MQ=255)                                                                                                                                cTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATc                    <  1:223847/149‑1 (MQ=255)                                                                                                                                    aatgCGTGAAGCCCAAGCCGGGCAAGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGttt                                                                              <  8:93391/84‑1 (MQ=25)                                                                                                                                    aatgCGTGAAGCCCAAGCCGGGCAAGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGttt                                                                              >  7:93391/4‑87 (MQ=25)                                                                                                                                      tgCGTGAAGCCCAAGCCGGGCAAGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCgacgac              >  2:21716/2‑149 (MQ=255)                                                                                                                                            aaGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACcagcag        >  5:121598/1‑149 (MQ=255)                                                                                                                                             aGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgccc                                           <  8:28429/113‑1 (MQ=255)                                                                                                                                             aGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgccc                                           >  7:28429/1‑113 (MQ=255)                                                                                                                                               cccTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATc                    >  4:257441/1‑134 (MQ=255)                                                                                                                                               cccTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATc                    <  3:257441/134‑1 (MQ=255)                                                                                                                                               cccTTGACCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGAGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTCCAACGGCTGATCGACGACCAGCAGGca     >  7:170711/1‑149 (MQ=255)                                                                                                                                                  ttGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAg  >  4:121703/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CTGTATGTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAG  >  NC_002947/5857815‑5858108

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