Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	421,161	T→A	5.5%	*929C (TGA→TGT)	PP_0348 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	421,161	0	T	A	5.5%	270.5 / 4.8	91	*929C (TGA→TGT)	PP_0348	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/43);  new base A (2/3);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

CCAGCACGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAGT  >  NC_002947/421014‑421309                                                                                                                                                    |                                                                                                                                                       ccAGCACGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGt                                                                                                                                                        <  3:29229/148‑1 (MQ=255)    gCACGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGt                                                                                                                                                    >  6:75899/1‑149 (MQ=255)      aCGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTgg                                                                                                                                                  <  6:285679/149‑1 (MQ=255)        gCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑Gc                                                                                                                                               <  5:327323/149‑1 (MQ=255)        gCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑Gc                                                                                                                                               >  1:53027/1‑149 (MQ=255)             ccAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGC‑CCACTGACAGCGGAGCcggcg                                                                                                                                          <  6:97408/149‑1 (MQ=255)              cagctcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCcggcgg                                                                                                                                         <  8:270621/149‑1 (MQ=255)                gctcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGt                                                                                                                                                    >  2:74054/1‑137 (MQ=255)                gctcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCggtg                                                                                                                                       <  1:284496/149‑1 (MQ=255)                  tcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCggtggt                                                                                                                                     >  2:151270/1‑149 (MQ=255)                    agctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTg                                                                                                                                   <  2:273660/149‑1 (MQ=255)                    agctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGACCGCCGGC‑CCACTGACAGCGGAGCcggcggcctc                                                                                                                                     <  7:245852/147‑5 (MQ=38)                           ccccGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTGGCCAAAACTGGTCGACTCAGGCCGCCGGCTCCGCTGTCTGGGG‑GCCGGCGGTGGTTGATCTGGc                                                                                                                            <  3:316158/149‑1 (MQ=255)                            cccGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCATTGG‑GCCGGCGGTGGTTGATCTgg                                                                                                                             >  4:81932/1‑147 (MQ=255)                            cccGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTgg                                                                                                                             >  5:300082/1‑147 (MQ=255)                             ccGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCTGTGGTTTATCTGGCAg                                                                                                                          <  8:99584/149‑1 (MQ=255)                                  cTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCa                                                                                                                     >  4:94258/1‑149 (MQ=255)                                   tATGCCCGTGACCGGGCACGGCTGCTACGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAg                                                                                                                    >  2:282372/1‑149 (MQ=255)                                     tGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑c                                                                                                                  <  5:75899/148‑1 (MQ=255)                                     tGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCACCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑c                                                                                                                  <  8:248944/148‑1 (MQ=255)                                        ccGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGAGTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAgt                                                                                                                         <  4:232566/139‑1 (MQ=255)                                        ccGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAgt                                                                                                                         >  3:232566/1‑139 (MQ=255)                                         cGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGt                                                                                                              >  6:195357/1‑148 (MQ=255)                                         cGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCAGCTCCGCTGTCAGTGG‑GCCGGCCGTGGTTGATCTGGCAGTGTCAG‑CCAGt                                                                                                              >  3:76312/1‑148 (MQ=255)                                             aCCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCActg                                                                                                         <  2:281679/149‑1 (MQ=255)                                              ccGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCActgc                                                                                                        >  8:118554/1‑149 (MQ=255)                                               cGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGGCAGTGG‑GCCGGCGGTGGTTTATCTGGCAGTGTCAG‑CCAGTCActgct                                                                                                       >  6:115333/1‑149 (MQ=255)                                                   cACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTgcgc                                                                                                   <  7:118554/149‑1 (MQ=255)                                                    aCGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑ACAGTCACTGATGCGcc                                                                                                  >  6:256201/1‑149 (MQ=255)                                                      ggCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCcct                                                                                                >  1:221954/1‑149 (MQ=255)                                                        ctgctgCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCcctgc                                                                                              <  2:221954/149‑1 (MQ=255)                                                          gctgcGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCcctgcct                                                                                            >  6:74644/1‑149 (MQ=255)                                                           ctgcGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCcctgcctg                                                                                           <  2:133131/149‑1 (MQ=255)                                                             gcgcAACTCCCCCACTGGCAGTGCCTTCAATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑CCCGGCGGTGGTTGATCTCGCAGTCTCAG‑CCAGTCACTGCTGCCCcctcaatggc                                                                                         >  8:252887/1‑142 (MQ=255)                                                              cgcAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTgg                                                                                          <  1:282372/147‑1 (MQ=255)                                                                 aaCTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGTTGTTGCCAATACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAAc                                                                                     <  5:323803/149‑1 (MQ=255)                                                                   cTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAAcgc                                                                                   <  5:74644/149‑1 (MQ=255)                                                                   cTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAAcgc                                                                                   <  8:107918/149‑1 (MQ=255)                                                                   cTGGCACACTGGCAGTGCCTTGCAGCAGCGGATGGACCGGCTGTTGCCAAACCTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGGCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAAcgc                                                                                   <  5:115333/149‑1 (MQ=255)                                                                          aCTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAg                                                                            >  7:183987/1‑149 (MQ=255)                                                                           cTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGc                                                                           >  5:137624/1‑149 (MQ=255)                                                                              gCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGATCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCAc                                                                        >  1:231759/1‑149 (MQ=255)                                                                                aGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTg                                                                      <  1:151270/149‑1 (MQ=255)                                                                                 gTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTgg                                                                     <  2:82697/149‑1 (MQ=255)                                                                                         cATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTc                                                             <  3:94258/149‑1 (MQ=255)                                                                                           tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCcctgcctg                                                                                           >  2:57606/1‑117 (MQ=255)                                                                                           tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCcctgcctg                                                                                           <  1:57606/117‑1 (MQ=255)                                                                                           tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCAc                                                                        >  4:139269/1‑136 (MQ=255)                                                                                           tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCAc                                                                        <  3:139269/136‑1 (MQ=255)                                                                                           tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGt                                                           >  3:36715/1‑149 (MQ=255)                                                                                           tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCATATCAGCCACTGGCCGTTTTCGt                                                           >  5:78933/1‑149 (MQ=255)                                                                                            cAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCATATCAGCCACTGGCCGTTTTCGt                                                           <  6:78933/148‑1 (MQ=255)                                                                                             aGCTGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTgg                                                         <  2:53027/149‑1 (MQ=255)                                                                                             aGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCAGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTgg                                                         <  4:76312/149‑1 (MQ=255)                                                                                                 gCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCGCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCgca                                                     <  6:300082/149‑1 (MQ=255)                                                                                                   tGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCgcagc                                                   >  3:245101/1‑149 (MQ=255)                                                                                                     gACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCgcagca                                                  <  5:195357/148‑1 (MQ=255)                                                                                                      aCCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCAcc                                                >  2:64902/1‑149 (MQ=255)                                                                                                        cGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTg                                                          <  6:218423/137‑1 (MQ=255)                                                                                                        cGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTg                                                          >  5:218423/1‑137 (MQ=255)                                                                                                              ttGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCAc                                          <  1:323482/147‑1 (MQ=255)                                                                                                                gCCAAAACTGGTCGACTGAGGCCGCCGCCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCACTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGca                                      >  7:37934/1‑149 (MQ=255)                                                                                                                 ccAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGca                                      <  8:284827/148‑1 (MQ=255)                                                                                                                    aaaCTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGcaccac                                  <  3:77696/149‑1 (MQ=255)                                                                                                                     aaCTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGACTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCcccccc                                 >  2:212964/1‑149 (MQ=255)                                                                                                                         ggTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCgg                             >  2:320860/1‑149 (MQ=255)                                                                                                                          gTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GACGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACAGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCCTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCCGc                            >  7:187752/1‑149 (MQ=255)                                                                                                                           tCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCg                           >  8:108691/1‑149 (MQ=255)                                                                                                                             gACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGtt                         >  2:142900/1‑149 (MQ=255)                                                                                                                             gACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGtt                         >  5:222801/1‑149 (MQ=255)                                                                                                                              aCTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑GCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATGAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTg                        <  3:135067/149‑1 (MQ=255)                                                                                                                              aCTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTg                        <  7:252887/149‑1 (MQ=255)                                                                                                                                tGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGcc                      >  1:322047/1‑148 (MQ=255)                                                                                                                                 gAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTg                     >  6:255312/1‑149 (MQ=255)                                                                                                                                 accaccgccgGC‑CCACTGACAGCGGAGCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTg                     >  4:130639/5‑149 (MQ=255)                                                                                                                                 accaccgccgGC‑CCACTGACAGCGGAGCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTg                     <  3:130639/145‑1 (MQ=255)                                                                                                                                  ccaccgccgGC‑CCACTGACAGCGGAGCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                    >  2:42797/4‑149 (MQ=255)                                                                                                                                  aGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                    <  2:303786/149‑1 (MQ=255)                                                                                                                                  aGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                    >  3:305039/1‑149 (MQ=255)                                                                                                                                  aGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                    >  3:310041/1‑149 (MQ=255)                                                                                                                                   ggCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGAGCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                    >  1:84916/1‑149 (MQ=255)                                                                                                                                     ccgccgGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCaa                 <  5:88287/149‑1 (MQ=255)                                                                                                                                     ccgcagGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGAACCGCCCACGGCACCACCCCGGCGTTGCTGGCac                 >  6:106351/1‑148 (MQ=255)                                                                                                                                       gccgGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCCAGCAACg               >  6:91630/1‑149 (MQ=255)                                                                                                                                         cgGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGc             <  4:245101/149‑1 (MQ=255)                                                                                                                                          ggCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCCAGCAACgccg            <  2:322047/149‑4 (MQ=255)                                                                                                                                           gCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCCAGCAACgccgg           <  5:91630/149‑5 (MQ=255)                                                                                                                                            cTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAt          <  6:137624/149‑1 (MQ=255)                                                                                                                                               cGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAtctc       >  5:149130/1‑149 (MQ=255)                                                                                                                                               cGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAtctc       <  4:36715/149‑1 (MQ=255)                                                                                                                                                gCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAtctc       >  1:253282/1‑148 (MQ=255)                                                                                                                                                  tGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCaa    <  6:222801/149‑1 (MQ=255)                                                                                                                                                    tCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAGt  <  1:142900/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       CCAGCACGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAG‑CCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAGT  >  NC_002947/421014‑421309

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑