Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,758,656 A→T 8.8% G402G (GGT→GGA) ‡ selB ← selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,758,6560AT 8.8% 299.3 / 15.6 102G402G (GGT→GGA) ‡selBselenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base A (45/48);  new base T (6/3);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

Only 100 of 102 total aligned reads displayed.

TTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCATTGC  >  NC_000913/3758508‑3758803
                                                                                                                                                    |                                                                                                                                                   
ttGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAAc                                                                                                                                                    <  1:1059325/150‑1 (MQ=255)
ttGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAAc                                                                                                                                                    <  1:162807/150‑1 (MQ=255)
  ggCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCa                                                                                                                                                  >  2:345947/1‑150 (MQ=255)
  ggCAACGCCATACGTCGCAGACGTTCGCGCCCAGAGCCAGGTTAATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCa                                                                                                                                                  >  2:1014677/1‑150 (MQ=255)
     aaCGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGc                                                                                                                                               <  1:63210/150‑1 (MQ=255)
     aaCGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGc                                                                                                                                               <  2:1206460/150‑1 (MQ=255)
     aaCGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGc                                                                                                                                               <  2:755606/150‑1 (MQ=255)
        gCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAAGCAGGCTGt                                                                                                                                            <  1:1635418/150‑1 (MQ=255)
        gCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGt                                                                                                                                            >  2:1205224/1‑150 (MQ=255)
         ccATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGGTTATATTCAGGCTggt                                                                                                                                           <  2:143536/150‑3 (MQ=255)
           aTACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGGTTATATTCAGGCTggtta                                                                                                                                         <  2:169196/150‑5 (MQ=255)
           aTACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTgc                                                                                                                                         <  2:1625141/150‑1 (MQ=255)
           aTACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTgc                                                                                                                                         <  2:659030/150‑1 (MQ=255)
           aTACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTgc                                                                                                                                         >  1:1576637/1‑150 (MQ=255)
           aTACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTgc                                                                                                                                         <  2:771912/150‑1 (MQ=255)
            tACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTgca                                                                                                                                        >  1:691661/1‑150 (MQ=255)
              cGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTgcagc                                                                                                                                      <  1:1199918/150‑1 (MQ=255)
              cGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTgcagc                                                                                                                                      >  1:631408/1‑150 (MQ=255)
                   cAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTc                                                                                                                                 >  2:1050328/1‑150 (MQ=255)
                           cgcgCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCAAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTCTTGCAGCAAATCGCGCATcc                                                                                                                         >  2:563396/1‑150 (MQ=255)
                               cccAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAAAGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTc                                                                                                                     >  1:1111430/1‑150 (MQ=255)
                                ccAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCg                                                                                                                    <  2:1006506/150‑1 (MQ=255)
                                 cAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGc                                                                                                                   <  2:1576637/150‑1 (MQ=255)
                                 cAGGGCCAGGTTAATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGc                                                                                                                   >  2:784291/1‑150 (MQ=255)
                                    ggCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGt                                                                                                                <  2:116096/150‑1 (MQ=255)
                                     gCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGtt                                                                                                               >  2:183292/1‑150 (MQ=255)
                                      ccAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTg                                                                                                              >  2:870874/1‑150 (MQ=255)
                                        aGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAg                                                                                                            <  1:1205224/150‑1 (MQ=255)
                                         ggTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGc                                                                                                           <  2:631408/150‑1 (MQ=255)
                                              aTCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACGGGCGCATTCAACAAGCTATAACCAGCCTGAATATTACCAGGCTTTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgcg                                                                                                      <  1:563396/150‑1 (MQ=255)
                                              aTCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgcg                                                                                                      >  1:1478384/1‑150 (MQ=255)
                                              aTCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgcg                                                                                                      >  1:1560338/1‑150 (MQ=255)
                                               tCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgcgc                                                                                                     >  2:1623403/1‑150 (MQ=255)
                                                        tGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCg                                                                                            <  2:691661/150‑1 (MQ=255)
                                                          cTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGaa                                                                                          >  2:128992/1‑150 (MQ=255)
                                                          cTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGaa                                                                                          >  1:858881/1‑150 (MQ=255)
                                                          cTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGaa                                                                                          >  1:218270/1‑150 (MQ=255)
                                                           tCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGaaa                                                                                         <  1:1014677/150‑1 (MQ=255)
                                                           tCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGaaa                                                                                         >  1:331504/1‑150 (MQ=255)
                                                            cATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAAt                                                                                        >  1:454962/1‑150 (MQ=255)
                                                              tGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATcc                                                                                      >  2:1047386/1‑150 (MQ=255)
                                                               gATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgc                                                                                                       <  1:289651/132‑1 (MQ=255)
                                                               gATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgc                                                                                                       >  2:289651/1‑132 (MQ=255)
                                                               gATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCg                                                                                     >  2:142873/1‑150 (MQ=255)
                                                                  aaGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCaa                                                                                  <  2:229654/150‑1 (MQ=255)
                                                                  aaGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCAGGCGAAATCCGCAAg                                                                                 <  1:805221/150‑1 (MQ=255)
                                                                   aGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAg                                                                                 >  1:207580/1‑150 (MQ=255)
                                                                    gTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAgg                                                                                >  1:350707/1‑150 (MQ=255)
                                                                      cGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGtt                                                                              >  2:378390/1‑150 (MQ=255)
                                                                      cGCTAATGTGACGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGtt                                                                              >  2:1083089/1‑150 (MQ=255)
                                                                         tAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAc                                                                           <  2:331504/150‑1 (MQ=255)
                                                                          aaTGTGGCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAcc                                                                          <  2:789006/150‑1 (MQ=255)
                                                                             gtgtCGAGAATTTTCCGCTGCCAGCGGGCGGCAACAGGCGCATTCAACAAGCTATAACCAGCATGAATATAACAAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcg                                                                       >  2:662810/1‑150 (MQ=255)
                                                                                tCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccg                                                                    <  2:362326/150‑1 (MQ=255)
                                                                                 cGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                                   >  1:1249038/1‑150 (MQ=255)
                                                                                  gagaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgcg                                                                  >  1:372715/1‑150 (MQ=255)
                                                                                    gaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGtt                                                                <  1:784291/150‑1 (MQ=255)
                                                                                      aTTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGGTTATATTCAGGCTggtt                                                                                                                                          <  2:1146087/74‑4 (MQ=37)
                                                                                      aTTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGGTTATATTCAGGCTggtt                                                                                                                                          >  1:1146087/1‑71 (MQ=37)
                                                                                      aTTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTcc                                                              <  1:1083089/150‑1 (MQ=255)
                                                                                       ttttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACTGCGCCGCGTTCCa                                                             <  2:858881/150‑1 (MQ=255)
                                                                                       ttttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCa                                                             <  2:570784/150‑1 (MQ=255)
                                                                                       ttttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCa                                                             <  2:1560338/150‑1 (MQ=255)
                                                                                         ttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGa                                                           <  1:1669345/150‑1 (MQ=255)
                                                                                         ttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGa                                                           <  1:1047386/150‑1 (MQ=255)
                                                                                         ttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGa                                                           <  1:378390/150‑1 (MQ=255)
                                                                                           ccGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATg                                                         <  2:1477684/150‑1 (MQ=255)
                                                                                            cGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAAGAAGGCTGTTGCAGCAATTAGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCTAAATGCGCAAGGTTAACCGCGCCGCGTTCCAGATGa                                                        <  1:144535/150‑1 (MQ=255)
                                                                                              cTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAAc                                                      >  1:170793/1‑150 (MQ=255)
                                                                                              cTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAAc                                                      <  1:345947/150‑1 (MQ=255)
                                                                                               tGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACa                                                     >  2:153715/1‑150 (MQ=255)
                                                                                                 ccAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGa                                                   <  1:128992/150‑1 (MQ=255)
                                                                                                  cAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGAt                                                  >  2:88232/1‑150 (MQ=255)
                                                                                                   aGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATa                                                 <  2:193582/150‑1 (MQ=255)
                                                                                                     cgggcggCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAAc                                               <  1:492634/150‑1 (MQ=255)
                                                                                                              aCCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACgcatcggca                                      >  1:43863/1‑150 (MQ=255)
                                                                                                                cGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACgcatcggtatc                                    <  1:531447/150‑1 (MQ=255)
                                                                                                                   cgcATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGct                                 >  2:169655/1‑150 (MQ=255)
                                                                                                                       ttCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgc                             <  1:358260/150‑1 (MQ=255)
                                                                                                                         caacaaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgcgc                           >  1:359758/1‑150 (MQ=255)
                                                                                                                         caacaaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgcgc                           >  1:452084/1‑150 (MQ=255)
                                                                                                                           acaaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGccc                         >  1:1035835/1‑150 (MQ=255)
                                                                                                                           acaaGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGccc                         >  1:1239522/1‑150 (MQ=255)
                                                                                                                              aGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTg                      <  2:1338624/150‑1 (MQ=255)
                                                                                                                              aGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTg                      <  2:43863/150‑1 (MQ=255)
                                                                                                                                gctgcaacAGCCTGGTTATATTCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCa                    >  1:1022493/8‑150 (MQ=255)
                                                                                                                                gctgcaacAGCCTGGTTATATTCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCa                    >  2:581418/8‑150 (MQ=255)
                                                                                                                                    aaCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAagag                <  2:113492/150‑1 (MQ=255)
                                                                                                                                     aacAGCCTGGTTATATTCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAgaga               >  2:286724/3‑150 (MQ=255)
                                                                                                                                     aCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAgaga               >  1:1351668/1‑150 (MQ=255)
                                                                                                                                     aCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAgaga               >  1:1346094/1‑150 (MQ=255)
                                                                                                                                      ccAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGAc              <  1:1050328/150‑1 (MQ=255)
                                                                                                                                      ccAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCACATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGAc              >  2:1073499/1‑150 (MQ=255)
                                                                                                                                      ccAGCCTCAAAATAACCATGCTCTTCCAGCAATTCGCGCATCCCTTCGCCGTTGACCTGGCGCGCCCAGGCGAAATCCGCAAGCTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGAc              >  2:1446663/1‑150 (MQ=255)
                                                                                                                                       cAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACg             >  2:469095/1‑150 (MQ=255)
                                                                                                                                        aGCCTGGTTATATTCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACgc            >  1:1219342/1‑150 (MQ=255)
                                                                                                                                         gCCTGGTTATATTCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACgcc           >  1:610774/1‑150 (MQ=255)
                                                                                                                                          ccTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACgcca          <  1:469095/150‑1 (MQ=255)
                                                                                                                                                tataACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCAtt    <  2:1478384/150‑1 (MQ=255)
                                                                                                                                                  taACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCATTGc  <  1:579307/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCATTGC  >  NC_000913/3758508‑3758803

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