Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 583,846 C→T 6.2% T571M (ACG→ATG)  selB → selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947583,8460CT 6.2% 191.0 / 2.5 64T571M (ACG→ATG) selBselenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base T (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGACCGCATCGGCCTGACCCGGAGG  >  NC_002947/583705‑583994
                                                                                                                                             |                                                                                                                                                    
ggCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAg                                                                                                                                               >  7:278069/1‑149 (MQ=255)
   aaTTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGATGCTGCAgcg                                                                                                                                             >  4:18263/1‑148 (MQ=255)
    tcctGACCCAGGGTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCa                                                                                                                                           >  2:7567/4‑149 (MQ=255)
    tcctGACCCAGGGTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCAACAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCa                                                                                                                                           >  8:113272/4‑149 (MQ=255)
      ctGACCCAGGGTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATg                                                                                                                                         >  7:77332/2‑149 (MQ=255)
       tGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATg                                                                                                                                         <  4:301452/148‑1 (MQ=255)
              gggtGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGtt                                                                                                                                 <  1:7567/146‑1 (MQ=255)
               ccTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCAt                                                                                                                                          >  1:231760/1‑139 (MQ=255)
               ccTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCAt                                                                                                                                          <  2:231760/139‑1 (MQ=255)
                gtGGGTCAAGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGa                                                                                                                               <  1:214145/148‑1 (MQ=255)
                 tGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTg                                                                                                                                                  <  3:18512/129‑1 (MQ=255)
                 tGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTg                                                                                                                                                  >  4:18512/1‑129 (MQ=255)
                  gggTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATg                                                                                                                             >  6:11625/1‑149 (MQ=255)
                     tCAGGACGCTGGCGAGTGAGGAGAACTGGGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTg                                                                                                                          <  4:267881/149‑1 (MQ=255)
                       aGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTggg                                                                                                                        >  5:102492/1‑149 (MQ=255)
                       aGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCGGCAGCGCATGGCAGCGTTGATGCTggg                                                                                                                        <  5:12426/149‑1 (MQ=255)
                         gACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTgggg                                                                                                                       >  2:40681/1‑148 (MQ=255)
                           cGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAg                                                                                                                    >  1:71289/1‑149 (MQ=255)
                                gCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAg                                                                                                               <  8:77332/149‑1 (MQ=255)
                                gCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAg                                                                                                               >  7:244728/1‑149 (MQ=255)
                                     tGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCAACAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGaaa                                                                                                          <  7:113272/149‑1 (MQ=255)
                                         gagaACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGGTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACgga                                                                                                      <  8:244728/149‑3 (MQ=255)
                                         gagaACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTAGATGCTGGGGCAGGCCAGCGAAACg                                                                                                        >  4:286813/1‑147 (MQ=255)
                                           gaACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGa                                                                                                    <  5:11625/149‑1 (MQ=255)
                                                  gCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTg                                                                                               >  7:38231/1‑147 (MQ=255)
                                                   ccGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAg                                                                                            <  6:145958/149‑1 (MQ=255)
                                                   ccGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAg                                                                                            >  6:168772/1‑149 (MQ=255)
                                                    cGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGATGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAgg                                                                                           >  4:231973/1‑149 (MQ=255)
                                                    cGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAgg                                                                                           <  7:133704/149‑1 (MQ=255)
                                                      aaGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGtt                                                                                         >  5:159690/1‑149 (MQ=255)
                                                               tgCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGtt                                                                                 <  8:278069/148‑1 (MQ=255)
                                                                   ccTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTcgc                                                                             >  2:94623/1‑148 (MQ=255)
                                                                     tGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTcgc                                                                             >  2:283503/1‑146 (MQ=255)
                                                                     tGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTcgc                                                                             <  1:283503/146‑1 (MQ=255)
                                                                         aCTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTg                                                                                                                          >  8:293018/1‑97 (MQ=255)
                                                                         aCTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTg                                                                                                                          <  7:293018/97‑1 (MQ=255)
                                                                         aCTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGc                                                                                                                      <  4:94617/101‑1 (MQ=255)
                                                                         aCTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGc                                                                                                                      >  3:94617/1‑101 (MQ=255)
                                                                          cTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCg                                                                                                                                   <  8:272614/87‑1 (MQ=255)
                                                                          cTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCg                                                                                                                                   >  7:272614/1‑87 (MQ=255)
                                                                             cgcAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAg                                                                                                                    <  6:247436/99‑1 (MQ=255)
                                                                             cgcAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAg                                                                                                                    >  5:247436/1‑99 (MQ=255)
                                                                             cgcAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTggg                                                                  >  2:76309/1‑149 (MQ=255)
                                                                             cgcAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTgg                                                                   <  5:168772/148‑1 (MQ=255)
                                                                              gcAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTggg                                                                  >  5:260984/1‑148 (MQ=255)
                                                                               cAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATg                                                                       >  3:83426/1‑142 (MQ=255)
                                                                               cAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATg                                                                       <  4:83426/142‑1 (MQ=255)
                                                                                      gCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGATGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTc                                                          <  3:18263/148‑1 (MQ=255)
                                                                                      gCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTc                                                          <  2:8133/148‑1 (MQ=255)
                                                                                               ggCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTc                                                          <  3:277941/139‑1 (MQ=255)
                                                                                               ggCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTc                                                          >  4:277941/1‑139 (MQ=255)
                                                                                                     gTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGtgtg                                            >  7:210993/1‑147 (MQ=255)
                                                                                                          tCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTc                                     >  4:244391/1‑149 (MQ=255)
                                                                                                                    cgcgACTTGTTCTACCCCGAGGCGATGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACt                            <  3:231973/148‑1 (MQ=255)
                                                                                                                    cgcgACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACtt                           >  7:2857/1‑149 (MQ=255)
                                                                                                                       gACTTGTTCTACACCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGAGAAACGCAGATAGTGCAGCTAGCCGCGTTTCGCGATATGTTGGGCATCGGTCCCAAGCGCAGTGTGCAGATTCTGGAATACTTAGa                        >  7:14160/1‑149 (MQ=255)
                                                                                                                        aCTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGAc                       <  1:76309/149‑1 (MQ=255)
                                                                                                                        aCTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGAc                       <  2:71289/149‑1 (MQ=255)
                                                                                                                         cTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGAcc                      <  6:102492/149‑1 (MQ=255)
                                                                                                                             ttCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGACCGCAt                  <  8:2857/149‑1 (MQ=255)
                                                                                                                                tACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATATTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGACCGCATCgg               >  6:316061/1‑149 (MQ=255)
                                                                                                                                 aCCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGACCGCATCGGc              >  5:36102/1‑149 (MQ=255)
                                                                                                                                           gACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGACCGCATCGGCCTGACCCg      >  8:121961/1‑147 (MQ=255)
                                                                                                                                           gACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGACCGCATCGGCCTGACCCg      >  8:252797/1‑147 (MQ=255)
                                                                                                                                             cGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGACCGCATCGGCCTGACCCGGAgg  <  8:38231/149‑1 (MQ=255)
                                                                                                                                             cGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGACCGCATCGGCCTGACCCGGAgg  <  6:159690/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GGCAATTTGACCCACCCTGGGTCAGGACGCTGGCGAGTGAGGAGAACTGTGCCGAAGCCGATGTGCGCCTGTTACTGCGCAAGTTGGCCCGGCTGGGCCTGGTGCATCAGGTGGTGCGCGACTTGTTCTACCCCGAGGCGACGCTGCAGCGCATGGCAGCGTTGATGCTGGGGCAGGCCAGCGAAACGCAGATAGTGCAGGTTGCCGCGTTTCGCGATATGTTGGGCATCGGTCGCAAGCGCAGTGTGCAGATTCTGGAATACTTCGACCGCATCGGCCTGACCCGGAGG  >  NC_002947/583705‑583994

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: