Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,713,789	G→T	5.5%	W674L (TGG→TTG)	PP_4171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,713,789	0	G	T	5.5%	235.2 / 3.2	73	W674L (TGG→TTG)	PP_4171	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base T (2/2);  total (35/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 = 4.53e-01

GGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTT  >  NC_002947/4713642‑4713935                                                                                                                                                    |                                                                                                                                                   ggttGCGCGAATGGGGTGCAGGCTGGTCGATCTCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGATCGGCTGGACGTGCGCCTGGCGATGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCtgtg                                                                                                                                                    <  4:251833/145‑1 (MQ=255)  gATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTgg                                                                                                                                                   >  8:62402/1‑148 (MQ=255)      cgcgAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgcagc                                                                                                                                               <  2:296525/148‑1 (MQ=255)      cgcgAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGgcag                                                                                                                                                >  5:195901/1‑147 (MQ=255)        cgAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCagcagc                                                                                                                                            <  6:195901/149‑1 (MQ=255)        cgAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCagcagc                                                                                                                                            <  8:273689/149‑1 (MQ=255)          aaTGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAgcgc                                                                                                                                          >  6:232673/1‑149 (MQ=255)          aaTGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAgcgc                                                                                                                                          >  1:203623/1‑149 (MQ=255)            tGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCt                                                                                                                                        <  7:30744/149‑1 (MQ=255)              gggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTg                                                                                                                                       <  5:65071/148‑1 (MQ=255)              gggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGCTGCCACAGCGccgcc                                                                                                                                      <  8:100668/149‑4 (MQ=255)               ggTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGa                                                                                                                                     >  2:38429/1‑149 (MQ=255)                   ccGTCTGGTCGAGCGCGCTGTTATATTACGTTGACGATGCCCTCATCCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGAt                                                                                                                                 <  8:332476/149‑1 (MQ=255)                   ccGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGAt                                                                                                                                 <  8:240836/149‑1 (MQ=255)                     ggCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATg                                                                                                                                >  8:170552/1‑148 (MQ=255)                         ggTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGAGGCGCTGCTGCCACAGCGccgccagc                                                                                                                                   >  8:219567/1‑135 (MQ=255)                         ggTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGAGGCGCTGCTGCCACAGCGccgccagc                                                                                                                                   <  7:219567/141‑7 (MQ=255)                             gAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGCACCCACCt                                                                                                                       >  1:157211/1‑149 (MQ=255)                             gAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCt                                                                                                                       <  1:271415/149‑1 (MQ=255)                                   gcTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGccc                                                                                                                 >  8:303311/1‑149 (MQ=255)                                       ttCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGc                                                                                                              <  7:62402/148‑1 (MQ=255)                                           aTTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAAc                                                                                                         >  1:37446/1‑149 (MQ=255)                                            ttACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGa                                                                                                           <  5:232673/146‑1 (MQ=255)                                            ttACGGTGACGATGCCCTCAACCAGTACCGCTAAGAGAGTCTGGACCTGCGCATGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAAcc                                                                                                        >  6:190573/1‑149 (MQ=255)                                               cGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGc                                                                                                      >  1:198778/1‑148 (MQ=255)                                                    aCGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGa                                                                                                 >  8:14998/1‑148 (MQ=255)                                                     cGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCATCCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACa                                                                                               <  1:360895/149‑1 (MQ=255)                                                      gATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAg                                                                                              <  2:203623/149‑1 (MQ=255)                                                      gATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCAGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCCGCGCTGTCCCAGCAGCGCCTGGACGATGACCCCACCTCCGCCCTCACCAACCCCTACGACAg                                                                                              >  8:131232/1‑149 (MQ=255)                                                               aaCCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCg                                                                                                       <  3:194652/131‑1 (MQ=255)                                                               aaCCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCg                                                                                                       >  4:194652/1‑131 (MQ=255)                                                                aCCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTATATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGc                                                                                    >  7:247653/1‑149 (MQ=255)                                                                    gTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCt                                                                                  >  4:87677/1‑147 (MQ=255)                                                                       ccGCTCCGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagc                                                                              <  4:305462/148‑1 (MQ=255)                                                                        cGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGc                                                                                                      <  6:167074/123‑1 (MQ=255)                                                                        cGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGc                                                                                                      >  5:167074/1‑123 (MQ=255)                                                                          cTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgt                                                                            >  4:343261/1‑147 (MQ=255)                                                                           tACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCtgagcgtgag                                                                         >  8:325985/1‑149 (MQ=255)                                                                              gAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTAc                                                                                                   >  8:318492/1‑120 (MQ=255)                                                                              gAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTAc                                                                                                   <  7:318492/120‑1 (MQ=255)                                                                               aGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGcgcg                                                                     >  5:134098/1‑149 (MQ=255)                                                                                gCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCgg                                                                    <  4:119267/149‑1 (MQ=255)                                                                                    cTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCt                                                                <  7:170552/149‑1 (MQ=255)                                                                                       gACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCGGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGa                                                             <  7:131232/149‑1 (MQ=255)                                                                                               ccTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGAt                                                     <  5:190573/149‑1 (MQ=255)                                                                                               ccTGGCCAAGCGTTTTCGGGTCTATGGCAGCAACCTGGAGCGGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGAt                                                     <  8:247653/149‑1 (MQ=255)                                                                                                cTGGCCAAGCGTTTTCGGGTGTATGGCATCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATg                                                    <  7:303311/149‑1 (MQ=255)                                                                                                   gCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATgg                                                   <  3:56005/147‑1 (MQ=255)                                                                                                   gCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGcc                                                 <  7:14998/149‑1 (MQ=255)                                                                                                     cAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTg                                                               >  1:182784/1‑133 (MQ=255)                                                                                                     cAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTg                                                               <  2:182784/133‑1 (MQ=255)                                                                                                     cAAGCGTTTTCCCGTCTATGGCAGCAACCTGGAGCTGGCTGAGCTGTGACAGCAGCGCCTAGACGATGAACCAAACTCGGACCTGCAAAACCGCAACGACACCCCACACCGCCTAAGCCTGAGCGCGGAGCTCGAGTTCTGATgaacc                                                >  4:3820/1‑144 (MQ=255)                                                                                                      aaGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctg                                              >  5:311136/1‑149 (MQ=255)                                                                                                        gCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgct                                            >  5:32720/1‑149 (MQ=255)                                                                                                                gggTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGcc                                    <  2:37446/149‑1 (MQ=255)                                                                                                                gggTGTATGACAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGcc                                    >  2:192466/1‑149 (MQ=255)                                                                                                                       tGGCAGCACCCTGGAGCTGGCGGCGCGGGGGCAGCTGCGCCTGGACGATGAACCCACCTCTGCCCTGCCGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtgct                             <  2:157211/149‑1 (MQ=255)                                                                                                                       tGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtgct                             >  3:298541/1‑149 (MQ=255)                                                                                                                       tGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtgc                              <  2:198778/148‑1 (MQ=255)                                                                                                                         gcagcaACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTtgctgc                           >  6:258721/1‑149 (MQ=255)                                                                                                                           agcaACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAcc                                                                                      >  8:81414/1‑88 (MQ=255)                                                                                                                           agcaACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCAcc                                                                                      <  7:81414/88‑1 (MQ=255)                                                                                                                            gcaACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGctct                        <  6:134098/149‑1 (MQ=255)                                                                                                                             caACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGCCGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTg                       >  4:357578/1‑149 (MQ=255)                                                                                                                              aaCCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGc                      <  5:258721/149‑1 (MQ=255)                                                                                                                              aaCCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGc                      <  3:343261/149‑1 (MQ=255)                                                                                                                                   ggAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCtgtg                 >  1:224910/1‑149 (MQ=255)                                                                                                                                    gAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCtgtg                 <  4:298541/148‑1 (MQ=255)                                                                                                                                        tccaggcgCTGCTGCCACAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGt            >  8:280230/5‑149 (MQ=255)                                                                                                                                         ggcggcGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGtt           >  1:372628/1‑149 (MQ=255)                                                                                                                                          gcggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTg          <  6:32720/149‑1 (MQ=255)                                                                                                                                          gcggcgCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTg          <  3:87677/149‑1 (MQ=255)                                                                                                                                               gcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgctg                                           <  6:117784/111‑1 (MQ=255)                                                                                                                                               gcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCctgctg                                           >  5:117784/1‑111 (MQ=255)                                                                                                                                               gcTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGcc     >  5:268552/1‑149 (MQ=255)                                                                                                                                                   tgGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCttt  <  1:38429/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTATTACGGTGACGATGCCCTCAACCAGTACCGCTACGAGCGTCTGGACCTGCGCCTGGCCAAGCGTTTTCGGGTGTATGGCAGCAACCTGGAGCTGGCGGCGCTGTGGCAGCAGCGCCTGGACGATGAACCCACCTCGGCCCTGCAGAACCGCTACGACAGCCGTCACCGCCTGAGCGTGAGCGCGGAGCTGGAGTTCTGATGGCCCTGCTGCTGCGCCTGTTGCTGCTCTGCCTGTGGCCGTTGGGGCCTTT  >  NC_002947/4713642‑4713935

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