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,9180GC100.0% 669.8 / NA 202W116C (TGG→TGCaslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (93/109);  total (93/109)
*NC_0009133,980,9190CG99.5% 660.8 / ‑5.1 202P117A (CCG→GCG) aslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (93/108);  minor base T (0/1);  total (93/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 205 total aligned reads displayed.

GTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGT  >  NC_000913/3980832‑3981005
                                                                                      ||                                                                                      
 tAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtc                                                                                  >  2:2352643/1‑93 (MQ=255)
 tAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtc                                                                                  >  2:585865/1‑93 (MQ=255)
   gcgcTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                                                >  2:1989036/1‑93 (MQ=255)
   gcgcTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                                                <  2:3836629/93‑1 (MQ=255)
     gcTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTcgc                                                                              >  1:713833/1‑93 (MQ=255)
       tACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCgg                                                                            <  2:3067098/93‑1 (MQ=255)
        aCAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                                           <  1:1659535/93‑1 (MQ=255)
        aCAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                                           >  1:4700899/1‑93 (MQ=255)
        aCAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                                           <  2:1951621/93‑1 (MQ=255)
         cAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAg                                                                          <  1:2543050/93‑1 (MQ=255)
          aGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGc                                                                         <  1:3847232/93‑1 (MQ=255)
            gCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcat                                                                       <  2:2304782/93‑1 (MQ=255)
             cGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatc                                                                      <  1:1467944/93‑1 (MQ=255)
             cGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGGGAGGGTTTCTCGCTGAGcatc                                                                      >  1:2690904/1‑93 (MQ=255)
              gAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatca                                                                     >  1:1007260/1‑93 (MQ=255)
                aaTATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAtt                                                                   <  2:1728585/93‑1 (MQ=255)
                   atGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTc                                                                <  2:3945219/93‑1 (MQ=255)
                   atGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTc                                                                >  1:4291792/1‑93 (MQ=255)
                     ggTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                                              <  1:1952267/93‑1 (MQ=255)
                     ggTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                                              <  2:121691/93‑1 (MQ=255)
                     ggTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                                              >  2:4017205/1‑93 (MQ=255)
                     ggTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                                              <  2:569051/93‑1 (MQ=255)
                       tGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgt                                                            <  1:2650329/93‑1 (MQ=255)
                       tGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgt                                                            <  1:3203863/93‑1 (MQ=255)
                       tGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgt                                                            <  1:3381861/93‑1 (MQ=255)
                        gCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgtt                                                           <  2:3126302/93‑1 (MQ=255)
                          tGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTgg                                                         >  1:3227248/1‑93 (MQ=255)
                          tGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTgg                                                         <  1:1419631/93‑1 (MQ=255)
                          tGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTgg                                                         <  2:2851600/93‑1 (MQ=255)
                          tGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTgg                                                         <  2:2888519/93‑1 (MQ=255)
                           ggCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGt                                                        <  2:4289439/93‑1 (MQ=255)
                           ggCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGt                                                        >  2:4315789/1‑93 (MQ=255)
                             cAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGttt                                                      >  2:1178869/1‑93 (MQ=255)
                              aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtttg                                                     <  2:794120/93‑2 (MQ=255)
                              aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtttg                                                     <  1:844817/93‑2 (MQ=255)
                              aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtttg                                                     >  1:108721/1‑92 (MQ=255)
                               ggAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtttgt                                                    <  1:2959008/93‑3 (MQ=255)
                               ggAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtttgt                                                    >  2:192218/1‑91 (MQ=255)
                               ggAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtttgt                                                    >  2:3763587/1‑91 (MQ=255)
                                gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                   <  2:114857/93‑1 (MQ=255)
                                gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                   <  2:3941972/93‑1 (MQ=255)
                                gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                   <  1:3394463/93‑1 (MQ=255)
                                gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                   >  1:360038/1‑93 (MQ=255)
                                   gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCt                                                >  1:787473/1‑93 (MQ=255)
                                   gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCt                                                <  2:275009/93‑1 (MQ=255)
                                   gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCt                                                >  2:2479194/1‑93 (MQ=255)
                                    aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGCGCATCATTTTCTTGTTGGGTTGTCGCTg                                               >  2:3152855/1‑93 (MQ=255)
                                    aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTg                                               >  1:3587202/1‑93 (MQ=255)
                                    aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTg                                               <  2:3550403/93‑1 (MQ=255)
                                      aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGa                                             <  2:2181289/93‑1 (MQ=255)
                                      aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGa                                             >  2:4355172/1‑93 (MQ=255)
                                       aGTAACAGCTTCCAGACTAACGGCGTGCTGCTTGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGAt                                            <  1:153737/93‑1 (MQ=255)
                                       aGTAACAGCTTCCAGACAAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGAt                                            <  2:1430899/93‑1 (MQ=255)
                                         tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATgg                                          >  1:4425232/1‑93 (MQ=255)
                                          aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                                                <  1:169399/54‑1 (MQ=255)
                                          aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                                                >  2:169399/1‑54 (MQ=255)
                                          aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGc                                         >  2:102317/1‑93 (MQ=255)
                                          aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGc                                         <  2:2075786/93‑1 (MQ=255)
                                           aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGcc                                        <  2:2240038/93‑1 (MQ=255)
                                              gCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGc                                     >  1:1771500/1‑93 (MQ=255)
                                              gCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGc                                     >  2:2688343/1‑93 (MQ=255)
                                                ttCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCt                                   >  1:4378108/1‑93 (MQ=255)
                                                ttCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCt                                   <  1:519372/93‑1 (MQ=255)
                                                 tCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGATTTGTCGCTGGATGGCCCGCCTg                                  >  2:2581464/1‑93 (MQ=255)
                                                  ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                 <  1:1700408/93‑1 (MQ=255)
                                                   cAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTgag                                <  1:114033/93‑1 (MQ=255)
                                                   cAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTgag                                >  2:3296321/1‑93 (MQ=255)
                                                    cgACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTgaga                               <  2:1671886/92‑1 (MQ=255)
                                                     gACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGAt                              >  1:4015703/1‑93 (MQ=255)
                                                     gACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGAt                              >  2:1413619/1‑93 (MQ=255)
                                                     gACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGAt                              >  1:794730/1‑93 (MQ=255)
                                                     gACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGAt                              <  1:207565/93‑1 (MQ=255)
                                                       cTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATcc                            >  2:3735407/1‑93 (MQ=255)
                                                        tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGCTCca                           >  1:985802/1‑93 (MQ=255)
                                                        tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCca                           >  2:703207/1‑93 (MQ=255)
                                                         aaCGGCGTGCTGCTGGATTACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcac                          >  1:2391080/1‑93 (MQ=255)
                                                         aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcac                          >  1:1926592/1‑93 (MQ=255)
                                                           cGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaa                        <  1:942409/93‑1 (MQ=255)
                                                             gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                      <  1:49170/93‑1 (MQ=255)
                                                              cGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatca                     <  1:2486482/93‑1 (MQ=255)
                                                               gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaa                    <  2:896219/93‑1 (MQ=255)
                                                                tgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaat                   <  2:4057778/93‑1 (MQ=255)
                                                                 gcagctgGATGACGAAGGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCGGGATGGCCCGCCTGAGATCCACAATCAata                  <  2:4031759/90‑1 (MQ=255)
                                                                  ctgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAatat                 <  1:4062159/93‑1 (MQ=255)
                                                                  ctgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAatat                 <  2:1348381/93‑1 (MQ=255)
                                                                     ctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgc              <  2:2081172/93‑1 (MQ=255)
                                                                      tgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgcg             >  2:179889/1‑93 (MQ=255)
                                                                       ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGGCCGCCTTAGATCCACAATCAATATCGCGt            >  1:2092947/1‑93 (MQ=255)
                                                                        gATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTg           >  2:422750/1‑93 (MQ=255)
                                                                         aTGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGa          <  1:4196132/93‑1 (MQ=255)
                                                                          tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGAc         >  1:2725856/1‑93 (MQ=255)
                                                                           gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACt        <  1:39224/93‑1 (MQ=255)
                                                                           gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACt        >  2:3464422/1‑93 (MQ=255)
                                                                           gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACt        >  1:4351267/1‑93 (MQ=255)
                                                                           gACGAATGGTGCGCGTTTCTAGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCCATCTCGCGTGACt        >  1:2809519/1‑93 (MQ=255)
                                                                              gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     <  2:1498213/93‑1 (MQ=255)
                                                                              gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     <  2:1072683/93‑1 (MQ=255)
                                                                              gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     <  1:4595483/93‑1 (MQ=255)
                                                                              gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     <  2:2792999/93‑1 (MQ=255)
                                                                               aaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAg    <  1:4476491/93‑1 (MQ=255)
                                                                                      ||                                                                                      
GTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGT  >  NC_000913/3980832‑3981005

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