Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,693 C→G 5.6% P268R (CCG→CGG)  cca → fused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,202,6930CG 5.6% 368.9 / 10.5 120P268R (CCG→CGG) ccafused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base C (63/49);  new base G (3/4);  total (67/53)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

Only 100 of 129 total aligned reads displayed.

CCCTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCAAAA  >  NC_000913/3202544‑3202841
                                                                                                                                                     |                                                                                                                                                    
cccTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAGCTCTGGcc                                                                                                                                                      >  1:335515/1‑150 (MQ=255)
cccTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGcc                                                                                                                                                      <  1:1518826/150‑1 (MQ=255)
 ccTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcg                                                                                                                                                     <  2:852252/150‑1 (MQ=255)
  cTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgc                                                                                                                                                    >  2:1496720/1‑150 (MQ=255)
   tGCCAAGTGGCATCCGGAAATCGATACGGGTCTTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgcg                                                                                                                                                   <  1:1673219/150‑1 (MQ=255)
   tGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGTCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgcg                                                                                                                                                   <  2:1231924/150‑1 (MQ=255)
     ccAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtc                                                                                                                                                 <  2:1478194/150‑1 (MQ=255)
      cAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtca                                                                                                                                                >  2:1467920/1‑150 (MQ=255)
        aGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtcagc                                                                                                                                              <  1:69266/150‑3 (MQ=255)
         gTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCGCTCAATGGCGGCGATGCTGAGTCCGCAGGGCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatca                                                                                                                                             <  2:1611241/150‑1 (MQ=255)
          tGGCATCCTGATATCGATACTGTTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCAGTTTCGCAAGTTTAAGCCACGATCTCGGTAAAGGGATGACGCGGCCAGAGTTCTGTCGGCGtcagccc                                                                                                                                            <  2:1537012/150‑5 (MQ=255)
           ggCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGATGACGCGGCCAGAGTTCTGGCGGCGtcagccct                                                                                                                                           <  1:1259387/150‑6 (MQ=255)
            gCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACATTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatgg                                                                                                                                          <  2:1183977/150‑1 (MQ=255)
             cATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCGGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggt                                                                                                                                         >  1:167696/1‑150 (MQ=255)
             cATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggt                                                                                                                                         >  2:497351/1‑150 (MQ=255)
              aTCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                                                                                        >  1:384516/1‑150 (MQ=255)
                  ggAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtcatgg                                                                                                                                    >  2:1052567/1‑150 (MQ=255)
                  ggAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtcatgg                                                                                                                                    >  2:1857152/1‑150 (MQ=255)
                  ggAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTAGCAACTTTATGCCACGAGCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtcatgg                                                                                                                                    >  2:609106/1‑150 (MQ=255)
                   gAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGc                                                                                                                                   <  2:1038695/150‑1 (MQ=255)
                      aTCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcg                                                                                                                                >  1:1034150/1‑150 (MQ=255)
                       tCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcgg                                                                                                                               >  1:1493103/1‑150 (MQ=255)
                         gATACGGGTATTCATACGTTAATGACGCTCGCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGCCGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggcg                                                                                                                             >  2:1281141/1‑150 (MQ=255)
                                gTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCTTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                                                                      <  2:1358520/150‑1 (MQ=255)
                                gTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                                                                      <  1:609673/150‑1 (MQ=255)
                                 tATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTaa                                                                                                                     <  2:357305/150‑1 (MQ=255)
                                    tCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGtt                                                                                                                  >  2:1624967/1‑150 (MQ=255)
                                       tACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                                                                               <  2:808873/150‑1 (MQ=255)
                                        aCCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTg                                                                                                              >  1:1544313/1‑150 (MQ=255)
                                        aCCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTg                                                                                                              <  2:1776344/150‑1 (MQ=255)
                                          cTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGa                                                                                                            >  1:1333000/1‑150 (MQ=255)
                                              aTGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGGTTCGCAACTTTATGCCACGCTCTCGGTTAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGGCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                                                                        >  2:100506/1‑150 (MQ=255)
                                               tGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCTACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGCTCATGGCCCGGCGGGTGTTTAGTTAGTGGaacaac                                                                                                       >  1:450773/1‑150 (MQ=255)
                                               tGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaac                                                                                                       >  1:1328827/1‑150 (MQ=255)
                                               tGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaac                                                                                                       <  1:544519/150‑1 (MQ=255)
                                                gACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACt                                                                                                      >  1:480273/1‑150 (MQ=255)
                                                caCGCGCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTCTGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACt                                                                                                      <  2:138714/149‑1 (MQ=255)
                                                 aCGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTa                                                                                                     >  1:1057743/1‑150 (MQ=255)
                                                  cGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTAt                                                                                                    >  2:1829784/1‑150 (MQ=255)
                                                  cGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTAt                                                                                                    <  1:1842001/150‑1 (MQ=255)
                                                  cGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTAt                                                                                                    >  1:1244812/1‑150 (MQ=255)
                                                     tctcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGcc                                                                                                 >  2:632650/1‑150 (MQ=255)
                                                        cAATGGCGGCGATGCTGGGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGc                                                                                              >  2:898224/1‑150 (MQ=255)
                                                                    tGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGc                                                                                  >  1:1653672/1‑150 (MQ=255)
                                                                    tGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGAGCGGGc                                                                                  >  1:147627/1‑150 (MQ=255)
                                                                     gCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGcc                                                                                 <  2:1333000/150‑1 (MQ=255)
                                                                      cTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCGCGGCGGGTGTTAAGTTAGTGGAACACCTATGCCAGCGTCTGCGCGTGCCa                                                                                >  2:110601/1‑150 (MQ=255)
                                                                       tGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCaa                                                                               >  1:807203/1‑150 (MQ=255)
                                                                               cAGGTCGATGTCCGTTTCGCAACTTTATGCCGCGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAAtt                                                                       >  1:197357/1‑150 (MQ=255)
                                                                                 gggcGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCTTCTGCGCGTGCCAAATGAAATTcg                                                                     <  1:1235742/147‑1 (MQ=255)
                                                                                   gcGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgcg                                                                   <  2:1797163/149‑1 (MQ=255)
                                                                                   gcGATGTCCGTGTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGGCATGGCCCGGCGGGTGTTAAGTTAGGGGAACAACTATGCCAGCGGCTGCGCGTTCCAAATGAAATTcgcg                                                                   <  1:457158/149‑1 (MQ=255)
                                                                                     gATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAATTGAAATTCGCGAt                                                                 >  2:884078/1‑150 (MQ=255)
                                                                                        ggccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTa                                                              <  2:1114548/148‑1 (MQ=255)
                                                                                        gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGGCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTa                                                              <  2:1140916/150‑1 (MQ=255)
                                                                                          ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGc                                                            <  1:1091997/150‑1 (MQ=255)
                                                                                          cagTTTCGCAACTTTATGCCACGATCGCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGGCATGGCCCGGCGGGTGTTAAGGTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGc                                                            <  1:1158340/148‑1 (MQ=255)
                                                                                           cGTATCGCAACTTTATGCCCCGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGcc                                                           <  2:477776/150‑1 (MQ=255)
                                                                                               tCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGAc                                                       >  1:7970/1‑150 (MQ=255)
                                                                                                cGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGATAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACt                                                      >  1:955843/1‑150 (MQ=255)
                                                                                                   aaCTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATGTAGCCAGCCtggt                                                   >  2:1356220/1‑150 (MQ=255)
                                                                                                    catttATGCCACGATCTCGGTAACGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGGCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACtggtg                                                  <  2:132125/148‑1 (MQ=255)
                                                                                                          tGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAg                                            >  2:1221182/1‑150 (MQ=255)
                                                                                                          tGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAg                                            >  2:858083/1‑150 (MQ=255)
                                                                                                          tGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAg                                            <  1:872064/150‑1 (MQ=255)
                                                                                                            ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGCGtt                                          >  1:697046/1‑150 (MQ=255)
                                                                                                                   ctcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGag                                   >  1:1124743/1‑149 (MQ=255)
                                                                                                                   ctcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAt                                   >  2:536102/1‑150 (MQ=255)
                                                                                                                    tcGGTTAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAtc                                  <  2:852375/150‑1 (MQ=255)
                                                                                                                    tcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAtc                                  <  1:707646/150‑1 (MQ=255)
                                                                                                                    tcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAtc                                  <  2:411010/150‑1 (MQ=255)
                                                                                                                     cGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAtct                                 >  2:506037/1‑150 (MQ=255)
                                                                                                                     cGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAtcg                                 >  1:949948/1‑149 (MQ=255)
                                                                                                                      ggTAAAGGGCTGCCGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAtctc                                <  2:733172/150‑1 (MQ=255)
                                                                                                                        tAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCAt                              >  2:72031/1‑150 (MQ=255)
                                                                                                                        gaaaGGGCGCACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCAt                              <  2:1624888/149‑1 (MQ=255)
                                                                                                                          aaGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATcc                            <  1:1304424/150‑1 (MQ=255)
                                                                                                                              gCTGCCGCCGCCAGAACTCGGGCAGCGGCATCATGGTCAGGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCGGCGCGTTCACAATTAAATACGCGATTTCGCCCGACTGGTGGCTGAGGTTCACGATCGCATCCgcac                        >  2:104250/1‑150 (MQ=255)
                                                                                                                              gCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCacac                        >  2:1056180/1‑150 (MQ=255)
                                                                                                                                 gACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGTACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCGCATCCACACCtg                     >  1:329992/1‑149 (MQ=255)
                                                                                                                                   cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTcc                   >  1:396193/1‑150 (MQ=255)
                                                                                                                                   cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTcc                   >  2:599495/1‑150 (MQ=255)
                                                                                                                                   cgccgcCAGAACTCTGGCCGCGTCATCATGGGCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTcc                   <  1:1433475/150‑1 (MQ=255)
                                                                                                                                    gccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTccc                  <  1:1521829/150‑1 (MQ=255)
                                                                                                                                        ccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATg              >  2:681390/1‑150 (MQ=255)
                                                                                                                                         cAGAACTCTGGCCTCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCAGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGc             <  1:609106/150‑1 (MQ=255)
                                                                                                                                         cAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGc             <  1:1045870/150‑1 (MQ=255)
                                                                                                                                             aCTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGaa         <  1:858083/150‑1 (MQ=255)
                                                                                                                                              ctctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGCTCTCATCCACACCTTCCCAATGCTGAAc        >  2:1215247/1‑150 (MQ=255)
                                                                                                                                              ctctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAAc        >  1:789619/1‑150 (MQ=255)
                                                                                                                                              ctctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAAc        >  1:868879/1‑150 (MQ=255)
                                                                                                                                              ctctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGGTAAGTTTGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAAc        <  1:731678/150‑1 (MQ=255)
                                                                                                                                                 tGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCa     <  1:861563/150‑1 (MQ=255)
                                                                                                                                                  ggCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTCAACCCaa    >  1:213100/1‑150 (MQ=255)
                                                                                                                                                  ggCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGAGTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCaa    >  1:1637118/1‑150 (MQ=255)
                                                                                                                                                   gCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCaaa   <  2:4254/150‑1 (MQ=255)
                                                                                                                                                   gCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCaaa   >  1:1051595/1‑150 (MQ=255)
                                                                                                                                                    ccGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGCGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCataa  >  1:1698610/1‑150 (MQ=255)
                                                                                                                                                    ccGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCaaaa  <  1:1624967/150‑1 (MQ=255)
                                                                                                                                                    ccGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTGAGTGGACCAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGACCCCaaaa  >  2:726094/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCAAAA  >  NC_000913/3202544‑3202841

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: