Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	583,795	G→A	7.1%	R554Q (CGG→CAG)	selB →	selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	583,795	0	G	A	7.1%	203.8 / 4.6	70	R554Q (CGG→CAG)	selB	selenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base G (40/25);  new base A (2/3);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

GTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCG  >  NC_002947/583649‑583939                                                                                                                                                     |                                                                                                                                                 gTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGc                                                                                                                                                >  4:381518/1‑149 (MQ=255)     aGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTgg                                                                                                                                             <  3:215168/148‑1 (MQ=255)        gAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCt                                                                                                                                         >  8:229768/1‑149 (MQ=255)         aaGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTg                                                                                                                                        <  6:210517/149‑1 (MQ=255)         aaGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTg                                                                                                                                        >  4:178625/1‑149 (MQ=255)          aGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCAGCCGGGccaac                                                                                                                                       <  5:19077/149‑4 (MQ=255)           gCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGt                                                                                                                                      <  3:178625/149‑1 (MQ=255)            ccGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGCAAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGt                                                                                                                                      >  1:96904/1‑149 (MQ=255)                cAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATc                                                                                                                                 <  5:115306/149‑1 (MQ=255)                  gTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCAGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCa                                                                                                                                <  1:112016/149‑1 (MQ=255)                    gCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggt                                                                                                                             >  7:119572/1‑149 (MQ=255)                      aTTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCgg                                                                                                                                                 <  1:18410/127‑1 (MQ=255)                      aTTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCgg                                                                                                                                                 >  2:18410/1‑127 (MQ=255)                      aTTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggtgg                                                                                                                           >  3:325197/1‑149 (MQ=255)                      aTTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCAGCCGGGccaactt                                                                                                                                     <  3:36333/139‑6 (MQ=255)                      aTTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCAGCCGGGccaactt                                                                                                                                     >  4:36333/1‑134 (MQ=255)                          tgGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTgc                                                                                                                        <  7:229768/148‑1 (MQ=255)                             gCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggtg                                                                                                                            >  7:169056/1‑141 (MQ=255)                             gCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggtg                                                                                                                            <  8:169056/141‑1 (MQ=255)                              cACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGAc                                                                                                                   <  8:119572/149‑1 (MQ=255)                                cGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACtt                                                                                                                 <  1:369837/149‑1 (MQ=255)                                    tGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTc                                                                                                             >  7:124031/1‑149 (MQ=255)                                        gCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGtt                                                                                                              >  4:312040/1‑144 (MQ=255)                                        gCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGtt                                                                                                              <  3:312040/144‑1 (MQ=255)                                           gAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTAcccc                                                                                                       >  1:226127/1‑148 (MQ=255)                                            aaGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACtt                                                                                                                 >  3:301897/1‑137 (MQ=255)                                            aaGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACtt                                                                                                                 <  4:301897/137‑1 (MQ=255)                                              gTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTgg                                                                                                                                             <  2:237708/107‑1 (MQ=255)                                              gTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTgg                                                                                                                                             >  1:237708/1‑107 (MQ=255)                                                 gATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGa                                                                                                >  1:159708/1‑149 (MQ=255)                                                  aTCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGAc                                                                                               >  3:51275/1‑149 (MQ=255)                                                       aGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTg                                                                                           >  3:181396/1‑148 (MQ=255)                                                         ggC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggt                                                                                                                             <  8:358841/112‑1 (MQ=255)                                                         ggC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAggt                                                                                                                             >  7:358841/1‑112 (MQ=255)                                                         ggC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGAc                                                                                               >  6:232678/1‑142 (MQ=255)                                                         ggC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGAc                                                                                               <  5:232678/142‑1 (MQ=255)                                                           c‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAgcg                                                                                      >  7:155900/1‑149 (MQ=255)                                                               tttGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCAt                                                                                   >  2:123890/1‑149 (MQ=255)                                                                 tGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGACTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATg                                                                                  <  4:325197/148‑1 (MQ=255)                                                                  gACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACt                                                                                                                  >  7:347504/1‑115 (MQ=255)                                                                  gACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACt                                                                                                                  <  8:347504/115‑1 (MQ=255)                                                                  gACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATgg                                                                                 >  1:158852/1‑148 (MQ=255)                                                                  gACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATgg                                                                                 <  3:381518/148‑1 (MQ=255)                                                                    cccacccTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGAC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAg                                                                              >  7:80432/1‑149 (MQ=255)                                                                       acccTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGt                                                                           >  8:254553/1‑149 (MQ=255)                                                                                         ggCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGc                                                           >  4:378675/1‑147 (MQ=255)                                                                                            gAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCg                                                      <  2:96904/149‑1 (MQ=255)                                                                                              gTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGa                                                     >  8:226656/1‑148 (MQ=255)                                                                                                 aggagAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACg                                                 >  8:2516/1‑149 (MQ=255)                                                                                                   gagaACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGGGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCa                                               <  2:317795/149‑1 (MQ=255)                                                                                                      aaCTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGa                                             >  4:368325/1‑148 (MQ=255)                                                                                                           tgCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGt                                         <  7:2516/147‑1 (MQ=255)                                                                                                           tgCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTg                                        >  2:208622/1‑148 (MQ=255)                                                                                                            gCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTg                                        >  1:133961/1‑147 (MQ=255)                                                                                                             cc‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAg                                     <  3:378675/149‑1 (MQ=255)                                                                                                             cc‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAg                                     >  6:386006/1‑149 (MQ=255)                                                                                                             cc‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAg                                     >  4:250433/1‑149 (MQ=255)                                                                                                                   gCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGc                                <  1:123890/149‑1 (MQ=255)                                                                                                                      gATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCcgc                             >  8:143164/1‑149 (MQ=255)                                                                                                                           gcgcCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTc                        >  1:294688/1‑149 (MQ=255)                                                                                                                           gcgcCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTc                        <  2:158852/149‑1 (MQ=255)                                                                                                                                tGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGa                    >  2:196051/1‑148 (MQ=255)                                                                                                                                   tACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGatat                 <  2:133961/148‑1 (MQ=255)                                                                                                                                    aCTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGt               >  7:144401/1‑149 (MQ=255)                                                                                                                                      tGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCcgc                             <  7:124237/133‑1 (MQ=255)                                                                                                                                      tGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCcgc                             >  8:124237/1‑133 (MQ=255)                                                                                                                                       gcgcAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTgg            >  2:259218/1‑149 (MQ=255)                                                                                                                                         gcAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTggg           >  4:205931/1‑148 (MQ=255)                                                                                                                                           aaGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCAt        >  5:232657/1‑149 (MQ=255)                                                                                                                                            accaggCCCAGCCGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTgg            <  1:235834/140‑1 (MQ=255)                                                                                                                                            accaggCCCAGCCGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTgg            >  2:235834/5‑144 (MQ=255)                                                                                                                                                 gCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCg  >  6:267396/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GTTGAGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCG  >  NC_002947/583649‑583939

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