Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	895,198	G→C	5.6%	intergenic (+38/+26)	PP_0777 → / ← PP_0778	glutathione peroxidase/group 1 family glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	895,198	0	G	C	5.6%	230.9 / 3.1	71	intergenic (+38/+26)	PP_0777/PP_0778	glutathione peroxidase/group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base C (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

CGAGAAGTTTCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCAT  >  NC_002947/895052‑895342                                                                                                                                                    |                                                                                                                                                  cgtgaAGTTTCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGa                                                                                                                                                  <  2:391978/145‑1 (MQ=255)  gagaAGTTTCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGt                                                                                                                                                <  2:120464/149‑1 (MQ=255)        tttCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCcgcg                                                                                                                                          >  1:318323/1‑149 (MQ=255)              gtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCgcgggcg                                                                                                                                    >  1:377082/1‑149 (MQ=255)              gtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCgcgg                                                                                                                                       <  3:214060/146‑1 (MQ=255)                  gCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGg‑cgc                                                                                                                               >  5:397949/1‑149 (MQ=255)                   cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCt                                                                                                                              >  5:407537/1‑149 (MQ=255)                   cAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCt                                                                                                                              <  8:267509/149‑1 (MQ=255)                    aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTg                                                                                                                             >  2:415865/1‑149 (MQ=255)                    aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTg                                                                                                                             <  2:324468/149‑1 (MQ=255)                     aGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTgg                                                                                                                            <  6:397949/149‑1 (MQ=255)                      ggATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGa                                                                                                                           <  1:379238/149‑1 (MQ=255)                        aTGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGAtc                                                                                                                         <  1:333730/149‑1 (MQ=255)                        aTGGCCGGGTGCTGGCGCGTTTCGCGCCGCGCAAGGCGCCGGATGACCCAGCGGTTGTGCAGGCCATCGACAAACCCCTGGCCTGAATAGTGCCTTAGCGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGAt                                                                                                                          <  1:108475/149‑1 (MQ=255)                             cGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTca                                                                                                                      >  4:292969/1‑147 (MQ=255)                             cGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGCGCGG‑CACTCGATCTcaca                                                                                                                    <  1:317442/149‑1 (MQ=38)                              gggTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGTGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAg                                                                                                                   <  2:274502/149‑1 (MQ=255)                                gTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGgc                                                                                                                 <  7:24476/149‑1 (MQ=255)                                gTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGgc                                                                                                                 <  3:18833/149‑1 (MQ=255)                                  gCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGgcgc                                                                                                               <  3:23604/149‑1 (MQ=255)                                    tGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTg                                                                                                             <  8:362784/149‑1 (MQ=255)                                      gcgcGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTg                                                                                                             <  6:16281/147‑1 (MQ=255)                                      gcgcGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGc                                                                                                            >  4:361979/1‑148 (MQ=255)                                      gcgcGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGCGCGG‑CACTCGATCTca                                                                                                                      >  1:123817/1‑138 (MQ=38)                                      gcgcGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCCAGCGCCGCCCGCGCGG‑CACTCGATCTca                                                                                                                      <  2:123817/138‑1 (MQ=38)                                        gcgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCaaa                                                                                                         <  2:377082/149‑1 (MQ=255)                                         cgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCaaa                                                                                                         >  8:312229/1‑148 (MQ=255)                                         cgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAAc                                                                                                        <  6:407537/149‑1 (MQ=255)                                          gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTgg                                                                                                                            >  5:100675/1‑128 (MQ=255)                                          gTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTgg                                                                                                                            <  6:100675/128‑1 (MQ=255)                                                 gccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGGCCGCTGGATCTCACAGGCGCTGCAAACATTg                                                                                                    <  1:40157/146‑1 (MQ=255)                                                  ccgcGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTgcg                                                                                                  >  7:201622/1‑146 (MQ=255)                                                    gcgcACGGCGCCGTTTGACCCAGCGGTGGTGCAGGCAAGCCTGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGTATCTCAAAGGCGCTGCAAACATTGCGCAGaa                                                                                             <  8:201622/149‑1 (MQ=255)                                                       cACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGc                                                                                          <  1:112634/149‑1 (MQ=255)                                                            cgcCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGc                                                                                     >  2:17191/1‑149 (MQ=255)                                                               ctgATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGt                                                                                  <  7:312229/147‑1 (MQ=255)                                                               cGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCg                                                                                   >  5:46810/1‑148 (MQ=255)                                                                 gATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTCGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAg                                                                                <  5:379922/149‑1 (MQ=255)                                                                       ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                          <  6:2196/149‑1 (MQ=255)                                                                         aGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGt                                                                        >  1:303171/1‑149 (MQ=255)                                                                              tggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGcca                                                                   <  2:303171/149‑2 (MQ=255)                                                                              tggtgCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCg                                                                   <  1:17191/149‑1 (MQ=255)                                                                                   cAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTg                                                                                                    >  3:335588/1‑111 (MQ=255)                                                                                   cAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTg                                                                                                    <  4:335588/111‑1 (MQ=255)                                                                                   cAGGCCATAGAGAAAGCCCTGGCCTGAAGACTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATAGCGCAGAACGCCCAGCCGTAGCACCAGGATGCCggactt                                                              >  5:47093/1‑145 (MQ=255)                                                                                            gagaAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                       >  5:276148/1‑147 (MQ=255)                                                                                            gagaAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCa                                                     <  6:46810/149‑1 (MQ=255)                                                                                             agaAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAg                                                    <  8:93829/149‑1 (MQ=255)                                                                                               aaaGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGa                                                   >  7:7250/1‑148 (MQ=255)                                                                                                aaGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGAt                                                  <  3:361979/148‑1 (MQ=255)                                                                                                     cTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCt                                            >  2:422050/1‑149 (MQ=255)                                                                                                          cTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAAc                                       <  3:292969/149‑1 (MQ=255)                                                                                                          cTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAAc                                       >  2:289790/1‑149 (MQ=255)                                                                                                             aaGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGt                                    >  4:109833/1‑149 (MQ=255)                                                                                                                 gTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGa                                >  1:21901/1‑149 (MQ=255)                                                                                                                            tGATGTTTTGGCAACTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCgg                      >  1:161870/1‑149 (MQ=255)                                                                                                                                gTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAg                  >  4:169905/1‑149 (MQ=255)                                                                                                                                gTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAg                  >  2:299361/1‑149 (MQ=255)                                                                                                                                 ttttGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGc                 >  8:6290/1‑149 (MQ=255)                                                                                                                                      gCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg            >  4:324308/1‑149 (MQ=255)                                                                                                                                      gCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg            >  3:411702/1‑149 (MQ=255)                                                                                                                                      gCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg            >  3:123582/1‑149 (MQ=255)                                                                                                                                      gCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTgg            >  7:312430/1‑149 (MQ=255)                                                                                                                                       cACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggc           <  4:406616/149‑1 (MQ=255)                                                                                                                                       cACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggc           >  3:406616/1‑149 (MQ=255)                                                                                                                                        aCCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcg          >  3:27787/1‑149 (MQ=255)                                                                                                                                        aCCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcg          >  3:74369/1‑149 (MQ=255)                                                                                                                                          cTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcggc        >  7:197937/1‑149 (MQ=255)                                                                                                                                               agaTCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCa   <  5:198835/149‑1 (MQ=255)                                                                                                                                                gaTCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCAt  <  2:259700/149‑1 (MQ=255)                                                                                                                                                gaTCCAGCGCCGCCCGCGCGG‑CACTCGATCTCACAGGTGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCAt  >  7:151554/1‑149 (MQ=38)                                                                                                                                                    |                                                                                                                                                  CGAGAAGTTTCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGA‑CGTGATGTTTTGGCACCTGTGAGATCGAGTGCCGCGCGGGCGG‑CGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCAT  >  NC_002947/895052‑895342

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