Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,980,918 2 bp→CG coding (348‑349/1236 nt) aslB → putative arylsulfatase regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,9180GC99.6% 883.3 / ‑5.9 273W116C (TGG→TGCaslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base C (109/162);  minor base T (0/1);  total (109/164)
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.94e-01
*NC_0009133,980,9190CG100.0% 895.6 / NA 278P117A (CCG→GCG) aslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (109/169);  total (109/169)

Only 100 of 283 total aligned reads displayed.

CGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGT  >  NC_000913/3980831‑3981005
                                                                                       ||                                                                                      
cgtAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGtt                                                                                    <  1:2371973/93‑1 (MQ=255)
 gtAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGttt                                                                                   <  1:163933/93‑1 (MQ=255)
  tAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtc                                                                                  <  2:3303286/93‑1 (MQ=255)
   aGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtct                                                                                 >  1:1922489/1‑93 (MQ=255)
    gcgcTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                                                <  2:1658911/93‑1 (MQ=255)
    gcgcTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                                                >  1:4948230/1‑93 (MQ=255)
        tACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTACCAGCTTCCAGAGTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCgg                                                                            <  2:6513994/93‑1 (MQ=255)
        tACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCgg                                                                            >  2:4879861/1‑93 (MQ=255)
         aCAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                                           <  2:752216/93‑1 (MQ=255)
          cAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAg                                                                          >  2:3437628/1‑93 (MQ=255)
           aGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGc                                                                         >  1:3178978/1‑93 (MQ=255)
           aGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGc                                                                         >  1:3061840/1‑93 (MQ=255)
              cGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatc                                                                      <  1:3792692/93‑1 (MQ=255)
               gAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatca                                                                     >  2:360545/1‑93 (MQ=255)
                aaaTATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatcat                                                                    >  1:2138845/1‑93 (MQ=255)
                aaaTATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatcat                                                                    >  2:1003681/1‑93 (MQ=255)
                  atatGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAttt                                                                  <  1:3549488/93‑1 (MQ=255)
                     tGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCt                                                               <  2:2898343/93‑1 (MQ=255)
                     tGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCt                                                               >  2:4286387/1‑93 (MQ=255)
                     tGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCt                                                               <  1:2453432/93‑1 (MQ=255)
                      ggTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                                              >  2:6250658/1‑93 (MQ=255)
                        tGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgt                                                            <  1:6299631/93‑1 (MQ=255)
                         gCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgtt                                                           <  2:3087515/93‑1 (MQ=255)
                         gCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgtt                                                           <  1:4184626/93‑1 (MQ=255)
                         gCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgtt                                                           <  2:3380583/93‑1 (MQ=255)
                         gCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgtt                                                           <  2:5470669/93‑1 (MQ=255)
                            ggCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGt                                                        <  1:3289982/93‑1 (MQ=255)
                            ggCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGt                                                        <  1:4054986/93‑1 (MQ=255)
                             gCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtt                                                       <  2:6375372/93‑1 (MQ=255)
                             gCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtt                                                       <  2:5647561/93‑1 (MQ=255)
                              cAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGttt                                                      <  2:534394/93‑1 (MQ=255)
                              cAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGttt                                                      <  2:6550571/93‑1 (MQ=255)
                               aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtttg                                                     >  2:5845053/1‑92 (MQ=255)
                               aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtttg                                                     >  1:857261/1‑92 (MQ=255)
                                ggAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtttgt                                                    >  2:6072288/1‑91 (MQ=255)
                                ggAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtttgt                                                    >  2:3953035/1‑91 (MQ=255)
                                 gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                   <  2:6531245/93‑1 (MQ=255)
                                 gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                   >  1:615319/1‑93 (MQ=255)
                                  aaGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCg                                                  <  1:2713477/93‑1 (MQ=255)
                                  aaGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCg                                                  <  2:251705/93‑1 (MQ=255)
                                   aGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGc                                                 <  1:1205802/93‑1 (MQ=255)
                                    gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCt                                                >  1:2042896/1‑93 (MQ=255)
                                      tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTgg                                              <  1:2581001/93‑1 (MQ=255)
                                      tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTgg                                              <  1:290522/93‑1 (MQ=255)
                                      tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTgg                                              <  2:1246069/93‑1 (MQ=255)
                                        aGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGAt                                            <  1:5219457/93‑1 (MQ=255)
                                         gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                           >  1:60389/1‑93 (MQ=255)
                                         gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                           >  1:3902873/1‑93 (MQ=255)
                                          tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATgg                                          >  2:6629144/1‑93 (MQ=255)
                                          tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATgg                                          >  2:6099625/1‑93 (MQ=255)
                                           aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGc                                         <  1:31445/93‑1 (MQ=255)
                                           aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGc                                         >  1:5415338/1‑93 (MQ=255)
                                           aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGc                                         >  1:5199463/1‑93 (MQ=255)
                                            aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGcc                                        >  2:2037091/1‑93 (MQ=255)
                                             cAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGccc                                       <  1:1722206/93‑1 (MQ=255)
                                             cAGCCTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGccc                                       <  2:5135122/93‑1 (MQ=255)
                                              aGCTTCCAGACTAACGGCGTGCTGCTGGATGCCGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCg                                      >  1:3396409/1‑93 (MQ=255)
                                              aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCg                                      <  1:1694883/93‑1 (MQ=255)
                                              aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCg                                      >  2:185406/1‑93 (MQ=255)
                                               gCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGc                                     >  1:812064/1‑93 (MQ=255)
                                               gCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGc                                     >  1:3647077/1‑93 (MQ=255)
                                                cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGcc                                    >  2:2214423/1‑93 (MQ=255)
                                                cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGcc                                    <  2:1328663/93‑1 (MQ=255)
                                                cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGcc                                    <  2:4470735/93‑1 (MQ=255)
                                                 ttCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCt                                   >  1:1397278/1‑93 (MQ=255)
                                                   ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                 <  2:1699511/93‑1 (MQ=255)
                                                   ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                 <  2:2241434/93‑1 (MQ=255)
                                                   ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                 >  2:6303222/1‑93 (MQ=255)
                                                   ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                 <  1:4593389/93‑1 (MQ=255)
                                                            cGGCGTGCTGCTGGATGACGAATGGTGCGCATTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaa                        <  2:4063727/93‑1 (MQ=255)
                                                             ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaat                       >  1:2240477/1‑93 (MQ=255)
                                                              gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                      >  2:3903621/1‑93 (MQ=255)
                                                              gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                      <  1:5252257/93‑1 (MQ=255)
                                                              gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                      >  2:4302169/1‑93 (MQ=255)
                                                              gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                      <  1:3737826/93‑1 (MQ=255)
                                                                gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaa                    <  1:655310/93‑1 (MQ=255)
                                                                 tgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCGcaatcaat                   >  2:4562152/1‑93 (MQ=255)
                                                                   ctgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAatat                 >  1:5960832/1‑93 (MQ=255)
                                                                    tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                >  2:1755546/1‑93 (MQ=255)
                                                                    tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                <  2:530579/93‑1 (MQ=255)
                                                                     gctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcg               <  1:453009/93‑1 (MQ=255)
                                                                      ctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgc              >  1:5104627/1‑93 (MQ=255)
                                                                       tgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgcg             >  1:4376546/1‑93 (MQ=255)
                                                                        ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGt            <  1:736650/93‑1 (MQ=255)
                                                                        ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGt            <  1:4779667/93‑1 (MQ=255)
                                                                         gATGACGAATGGTGCGCGTTTCTCGCTGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTg           <  1:5614649/93‑1 (MQ=255)
                                                                           tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGAc         <  1:4711557/93‑1 (MQ=255)
                                                                           tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGAc         <  1:5079353/93‑1 (MQ=255)
                                                                            gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACt        <  1:5715785/93‑1 (MQ=255)
                                                                              cGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaa      <  1:6203392/93‑1 (MQ=255)
                                                                               gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     <  2:2712095/93‑1 (MQ=255)
                                                                               gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     <  1:1536562/93‑1 (MQ=255)
                                                                               gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     <  1:1061491/93‑1 (MQ=255)
                                                                                aaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAg    <  2:979994/93‑1 (MQ=255)
                                                                                 aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg   <  2:5182628/93‑1 (MQ=255)
                                                                                 aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg   <  1:675628/93‑1 (MQ=255)
                                                                                 aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg   <  1:556422/93‑1 (MQ=255)
                                                                                  tggtgcGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGt  <  2:6161420/93‑1 (MQ=255)
                                                                                  tggtgcGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGt  <  2:2556801/93‑1 (MQ=255)
                                                                                  tggtgcGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGt  <  1:5318334/93‑1 (MQ=255)
                                                                                       ||                                                                                      
CGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGT  >  NC_000913/3980831‑3981005

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: