Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	583,798	T→C	7.1%	L555P (CTG→CCG)	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,798	0	T	C	7.1%	200.5 / 3.8	70	L555P (CTG→CCG)	selB	selenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base T (40/25);  new base C (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 = 9.80e-01

AGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGC  >  NC_002947/583653‑583944                                                                                                                                                    |                                                                                                                                                   aGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTgg                                                                                                                                                  <  3:215168/148‑1 (MQ=255)    gAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCt                                                                                                                                              >  8:229768/1‑149 (MQ=255)     aaGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTg                                                                                                                                             >  4:178625/1‑149 (MQ=255)     aaGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTg                                                                                                                                             <  6:210517/149‑1 (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‑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                                                                                                                      <  4:301897/137‑1 (MQ=255)                                        aaGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACtt                                                                                                                      >  3:301897/1‑137 (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                                                                                                                       <  8:347504/115‑1 (MQ=255)                                                              gACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACt                                                                                                                       >  7:347504/1‑115 (MQ=255)                                                              gACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATgg                                                                                      <  3:381518/148‑1 (MQ=255)                                                              gACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATgg                                                                                      >  1:158852/1‑148 (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                                          >  6:386006/1‑149 (MQ=255)                                                                                                         cc‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAg                                          <  3:378675/149‑1 (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)                                                                                                                                                  gCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAgc  >  4:158292/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   AGCGAAGCCGACAGTGCATTGTGGGCACGTCTGCAGCCGAAGTTGATCGCAGGGC‑AATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCC‑GAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGC  >  NC_002947/583653‑583944

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