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% 870.6 / NA 248W116C (TGG→TGCaslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (111/137);  total (111/137)
*NC_0009133,980,9190CG99.6% 851.6 / NA 248P117A (CCG→GCG) aslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (110/137);  minor base T (0/1);  total (110/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

CGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGT  >  NC_000913/3980831‑3981005
                                                                                       ||                                                                                      
cgtAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGtt                                                                                    >  1:5678459/1‑93 (MQ=255)
 gtAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGttt                                                                                   <  2:5066239/93‑1 (MQ=255)
 gtAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGttt                                                                                   >  2:1974467/1‑93 (MQ=255)
  tAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtc                                                                                  >  2:2376718/1‑93 (MQ=255)
  tAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtc                                                                                  <  1:2413437/93‑1 (MQ=255)
   aGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtct                                                                                 >  1:805343/1‑93 (MQ=255)
   aGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtct                                                                                 <  1:1551221/93‑1 (MQ=255)
   aGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtct                                                                                 <  1:2873731/93‑1 (MQ=255)
    gcgcTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGttacac                                                                                <  1:411059/93‑5 (MQ=255)
    gcgcTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                                                >  1:506085/1‑93 (MQ=255)
    gcgcTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                                                >  1:4957653/1‑93 (MQ=255)
      gcTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTcgc                                                                              <  2:5224676/93‑1 (MQ=255)
       cTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTcgcg                                                                             >  2:4782626/1‑93 (MQ=255)
         aCAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                                           <  2:220150/93‑1 (MQ=255)
           aGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGc                                                                         <  1:3693100/93‑1 (MQ=255)
            ggCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGca                                                                        >  2:806323/1‑93 (MQ=255)
             gCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcat                                                                       <  2:5850180/93‑1 (MQ=255)
             gCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcat                                                                       <  2:5437741/93‑1 (MQ=255)
                aaaTATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatcat                                                                    >  1:5853365/1‑93 (MQ=255)
                aaaTATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatcat                                                                    <  2:3989315/93‑1 (MQ=255)
                 aaTATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAtt                                                                   >  1:4579435/1‑93 (MQ=255)
                 aaTATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAtt                                                                   <  2:1757197/93‑1 (MQ=255)
                  atatGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAttt                                                                  >  1:6159776/1‑93 (MQ=255)
                   tatGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAtttt                                                                 >  2:1328106/1‑93 (MQ=255)
                     tGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCt                                                               <  1:3886106/93‑1 (MQ=255)
                     tGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCt                                                               <  2:5624343/93‑1 (MQ=255)
                     tGGTGCTGGCAGGAAGATAAGTAACAGCATCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCt                                                               <  2:6208651/93‑1 (MQ=255)
                      ggTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                                              <  2:5705136/93‑1 (MQ=255)
                       gTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttg                                                             >  1:4736709/1‑93 (MQ=255)
                         gCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgtt                                                           >  1:1397949/1‑93 (MQ=255)
                          cTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgttg                                                          <  2:4685263/93‑1 (MQ=255)
                           tGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCCTTTTCTTGTTgg                                                         <  1:5152211/93‑1 (MQ=255)
                            ggCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGt                                                        <  1:5920458/93‑1 (MQ=255)
                             gCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtt                                                       <  2:1618141/93‑1 (MQ=255)
                             gCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtt                                                       <  2:6451922/93‑1 (MQ=255)
                               atgAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTa                                                     <  1:6259249/91‑1 (MQ=255)
                               aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTa                                                     <  2:2484455/93‑1 (MQ=255)
                                 gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATc                                                   <  1:5642807/93‑1 (MQ=255)
                                  aaGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCg                                                  >  1:3205981/1‑93 (MQ=255)
                                   aGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGc                                                 >  1:4335085/1‑93 (MQ=255)
                                    gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCt                                                <  2:2942971/93‑1 (MQ=255)
                                    gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCt                                                >  2:3121339/1‑93 (MQ=255)
                                    gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCt                                                >  1:3717222/1‑93 (MQ=255)
                                     aaaaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTg                                               <  1:3502277/91‑1 (MQ=255)
                                      tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTgg                                              >  2:6352542/1‑93 (MQ=255)
                                      tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTgg                                              >  2:86927/1‑93 (MQ=255)
                                      tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTgg                                              >  1:4192645/1‑93 (MQ=255)
                                       aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGa                                             >  1:3976162/1‑93 (MQ=255)
                                        aGTAACAGCTTCCGGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGAt                                            >  1:3124079/1‑93 (MQ=255)
                                        aGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGAt                                            <  2:1473823/93‑1 (MQ=255)
                                        aGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATAGCTGGAt                                            <  2:1639905/93‑1 (MQ=255)
                                         gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATg                                           <  1:1153907/93‑1 (MQ=255)
                                         gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCACTTTCTTGTTGGGTTATCGCTGGATg                                           >  2:4426176/1‑93 (MQ=255)
                                           aaCAGCTTCCAGACTAACGGCGTGCAGCTGGATGACGAATGGTGCGCGGTTCTCGCGGAGCATCATATTCTTGTTGGTTTATCGCTGGATgac                                         >  1:5201511/1‑91 (MQ=255)
                                            aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                                        >  2:6129941/1‑93 (MQ=255)
                                              aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcccc                                      >  2:4985388/1‑92 (MQ=255)
                                               gCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGc                                     >  2:2816101/1‑93 (MQ=255)
                                                cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGcc                                    <  1:5487016/93‑1 (MQ=255)
                                                 ttCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCt                                   >  2:1136627/1‑93 (MQ=255)
                                                 ttCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCt                                   <  2:1570087/93‑1 (MQ=255)
                                                 ttCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCt                                   <  1:3421457/93‑1 (MQ=255)
                                                 ttCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCt                                   >  2:1921797/1‑93 (MQ=255)
                                                 ttCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCt                                   <  2:428658/93‑1 (MQ=255)
                                                   ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTga                                 <  1:890287/93‑1 (MQ=255)
                                                   ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTga                                 >  2:3892728/1‑93 (MQ=255)
                                                   ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTga                                 <  1:3910132/93‑1 (MQ=255)
                                                    cAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTgag                                >  1:6228827/1‑93 (MQ=255)
                                                    cAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTgag                                >  1:2725582/1‑93 (MQ=255)
                                                       aCTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATc                             <  2:1293925/93‑1 (MQ=255)
                                                       aCTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATc                             >  2:776090/1‑93 (MQ=255)
                                                          aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcac                          >  1:1036597/1‑93 (MQ=255)
                                                          aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcac                          >  1:923729/1‑93 (MQ=255)
                                                           aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcaca                         >  2:3388245/1‑93 (MQ=255)
                                                           aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcaca                         <  1:1454402/93‑1 (MQ=255)
                                                             ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaat                       <  1:2126745/93‑1 (MQ=255)
                                                             ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaat                       <  1:2348542/93‑1 (MQ=255)
                                                             ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaat                       <  1:5742080/93‑1 (MQ=255)
                                                              gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatc                      >  1:672670/1‑93 (MQ=255)
                                                              gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatc                      >  1:1582769/1‑93 (MQ=255)
                                                              gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatc                      <  2:6425782/93‑1 (MQ=255)
                                                              gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatc                      <  2:4562531/93‑1 (MQ=255)
                                                                gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatcaa                    >  1:2553270/1‑93 (MQ=255)
                                                                  gctgctGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAata                  <  1:1696327/93‑1 (MQ=255)
                                                                    tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGGCCGCCTGAGATCCACAATCAACATc                >  2:3400820/1‑93 (MQ=255)
                                                                     gctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcg               >  2:6040740/1‑93 (MQ=255)
                                                                       tgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgcg             <  2:1193690/93‑1 (MQ=255)
                                                                        ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGt            <  1:4853609/93‑1 (MQ=255)
                                                                        ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGt            <  1:1876659/93‑1 (MQ=255)
                                                                         gATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTg           <  2:2096741/93‑1 (MQ=255)
                                                                             aCGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTa       <  2:2049929/93‑1 (MQ=255)
                                                                             aCGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTa       >  1:5386311/1‑93 (MQ=255)
                                                                             aCGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTa       <  2:1063561/93‑1 (MQ=255)
                                                                              cGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaa      >  1:4417799/1‑93 (MQ=255)
                                                                               gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     <  2:5955532/93‑1 (MQ=255)
                                                                               gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     <  2:2587451/93‑1 (MQ=255)
                                                                               gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTCGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     <  2:4528277/93‑1 (MQ=255)
                                                                               gAATGGGGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGCATCGTCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa     >  1:4896099/1‑93 (MQ=255)
                                                                                aaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAg    <  2:4341122/93‑1 (MQ=255)
                                                                                 aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg   <  2:3087844/93‑1 (MQ=255)
                                                                                 aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg   <  1:5214133/93‑1 (MQ=255)
                                                                                       ||                                                                                      
CGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGT  >  NC_000913/3980831‑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: