Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,117,238 G→A 80.6% M144I (ATG→ATA yceA → putative rhodanese‑related sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,117,2380GA80.6% 242.5 / 54.1 123M144I (ATG→ATAyceAputative rhodanese‑related sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base A (52/47);  total (63/60)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

Only 100 of 124 total aligned reads displayed.

GAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATAC  >  NC_000913/1117139‑1117333
                                                                                                   |                                                                                               
gAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATg                                                                                                 <  1:3468128/100‑1 (MQ=255)
   ggTACGCGATCGCATTGTTGCCGACGGTCTTGACGATCCTCACTTTGATGCCAGCAATGTTGGTAAGTATCTGCAAGCGGCGGAAGTGAACGCCATACtt                                                                                              >  1:1861597/1‑100 (MQ=255)
   ggTACGCGATCGCATTGTTGCCGACGGTATCGCCGATGCCCCGTTTGCTGCCACCCCTGTCGGCGAGCCTCTGCAAGCGGCCGCAGTGAACGCCATACtt                                                                                              >  2:3673999/1‑100 (MQ=255)
   gctACGCGATCGCTTTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACtt                                                                                              >  1:2831455/3‑100 (MQ=255)
    gTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTg                                                                                             >  1:2289199/1‑100 (MQ=255)
      aCGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGAc                                                                                           <  2:2892814/100‑1 (MQ=255)
      aCGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGAc                                                                                           <  1:3415957/100‑1 (MQ=255)
       cgcgATCGCATTGTTGCCGACGTTATTGACGATCCTCTCTTCGATGCCAGCAGTGTTGGTGTGTATCTGCAAACGGCGGAAGTGAACGCCATACTTGa                                                                                            >  2:2212212/1‑98 (MQ=255)
        gcgATCGCATTGTTTGCGACGGTATTGACGATCCTCACTTTGACGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAATGAACGCCATACTTAAc                                                                                           <  2:1729421/98‑1 (MQ=255)
         ctaTCGCTTTGGTGTTGATGGTATGGACGATCCGGACTTTGATGCCAGTAATGTTGGTGAGGATTTGCAATCGGCGGAAGTGAACGCCAGGCt                                                                                               <  2:882143/91‑1 (MQ=255)
         cgATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGAt                                                                                        <  1:421002/100‑1 (MQ=255)
          ggtggtaTTGTTGCTGACGGTATTGACGAGCCTCACTTTGATGCCAGCCAGGTTGGTGCGTGTCTGCAAGCGGCGGAAGTGAACGCCATGCTTGAc                                                                                           <  2:3444430/90‑1 (MQ=255)
          gATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGAt                                                                                        >  2:3949542/1‑99 (MQ=255)
           aTCGTATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATcc                                                                                      <  2:3935282/100‑1 (MQ=255)
           aTCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATcc                                                                                      <  1:4147555/100‑1 (MQ=255)
             cGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCg                                                                                    <  2:880309/100‑1 (MQ=255)
              gCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGa                                                                                   <  2:3681805/100‑1 (MQ=255)
              gCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGa                                                                                   >  1:3871699/1‑100 (MQ=255)
              gCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGa                                                                                   <  1:1558206/100‑1 (MQ=255)
              gCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGa                                                                                   >  1:1940917/1‑100 (MQ=255)
               cATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGAt                                                                                  >  1:3318472/1‑100 (MQ=255)
               cATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGAt                                                                                  <  2:4099625/100‑1 (MQ=255)
               cATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGAt                                                                                  <  2:790850/100‑1 (MQ=255)
               cATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGAt                                                                                  >  2:1374555/1‑100 (MQ=255)
                 ttgttgCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGc                                                                                >  2:1093571/1‑100 (MQ=255)
                   gttgCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCAc                                                                              >  1:719573/1‑100 (MQ=255)
                   gttgCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCAc                                                                              >  2:301748/1‑100 (MQ=255)
                   gttgCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCAc                                                                              >  2:805114/1‑100 (MQ=255)
                      gCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTAt                                                                           >  1:968055/1‑100 (MQ=255)
                       ccGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTAtt                                                                          <  2:298111/100‑1 (MQ=255)
                        cGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTAttt                                                                         >  2:2540656/1‑100 (MQ=255)
                            ggTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCg                                                                     <  1:1093574/100‑1 (MQ=255)
                            ggTATTGACGATCCTCACTTTGATGCCAGAAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATc                                                                      <  2:155408/99‑1 (MQ=255)
                             gTATTGACGATCCTCACTTTGATGCCAGCATTGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGa                                                                    >  1:83635/1‑100 (MQ=255)
                             gTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCg                                                                     >  2:899371/1‑99 (MQ=255)
                             gTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGa                                                                    >  1:2123656/1‑100 (MQ=255)
                              tATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGAc                                                                   >  2:2365729/1‑100 (MQ=255)
                                ttGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACAt                                                                 >  1:2376953/1‑100 (MQ=255)
                                 tGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATg                                                                >  2:2233746/1‑100 (MQ=255)
                                 tGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATg                                                                >  1:1643063/1‑100 (MQ=255)
                                   aCGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCg                                                              <  2:4010791/100‑1 (MQ=255)
                                   aCGATCCTCACTTTGATGCCAGCAATGGTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCg                                                              <  1:197322/100‑1 (MQ=255)
                                    cGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGt                                                             >  1:255293/1‑100 (MQ=255)
                                     gATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGtn                                                            <  1:1157526/100‑2 (MQ=255)
                                       tCCTCACTTTGATGCCAGCAATGTGGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCTATGCTTGACGATCCCGATGCATTGTTTATCGACATGCGTAAc                                                          <  1:2131921/100‑1 (MQ=255)
                                        ccTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAAcc                                                         <  1:2986691/100‑1 (MQ=255)
                                         cTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCa                                                        >  1:600624/1‑100 (MQ=255)
                                           cACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCACGATGCACTATTTATCGACATGCGTAACCACt                                                      >  2:933931/1‑100 (MQ=255)
                                            aCTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACTTGCGTAACCACTa                                                     <  1:2803063/100‑1 (MQ=255)
                                            aCTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTa                                                     >  1:2321590/1‑100 (MQ=255)
                                             cTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTAt                                                    >  1:3220906/1‑100 (MQ=255)
                                             cTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTAt                                                    >  1:371585/1‑100 (MQ=255)
                                             cTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTAt                                                    <  1:1739789/100‑1 (MQ=255)
                                              tttGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATg                                                   <  2:1973803/100‑1 (MQ=255)
                                               ttGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATg                                                   <  2:465101/99‑1 (MQ=255)
                                                tGATGCCAGCACTGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAg                                                 <  2:1861633/100‑1 (MQ=255)
                                                tGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAg                                                 <  1:164651/100‑1 (MQ=255)
                                                  aTGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTa                                               <  1:1123297/100‑1 (MQ=255)
                                                  aTGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTa                                               <  2:126592/100‑1 (MQ=255)
                                                   tGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTAt                                              <  2:2915165/100‑1 (MQ=255)
                                                     ccAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGa                                            <  2:2169220/100‑1 (MQ=255)
                                                      cAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGaa                                           >  1:3940054/1‑100 (MQ=255)
                                                          aaTGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTg                                        >  2:1064496/1‑99 (MQ=255)
                                                           aTGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTggg                                      <  1:2133211/100‑1 (MQ=255)
                                                             gTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGc                                    <  2:1250644/100‑1 (MQ=255)
                                                              ttGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTaggg                                     >  2:3974351/1‑98 (MQ=255)
                                                               tGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTAGGGCAc                                  >  2:3904938/1‑100 (MQ=255)
                                                                 gTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACtt                                >  1:1530511/1‑100 (MQ=255)
                                                                  tGAGTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACttt                               >  1:3535500/1‑100 (MQ=255)
                                                                     gTATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGaa                            <  2:3698519/100‑1 (MQ=255)
                                                                      tATCTGCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGaaa                           <  1:2015622/100‑1 (MQ=255)
                                                                          tGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAAcgc                       >  1:1583725/1‑100 (MQ=255)
                                                                          tGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAAcgc                       <  1:3066234/100‑1 (MQ=255)
                                                                          tGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAAcgc                       >  2:3462613/1‑100 (MQ=255)
                                                                           gCAAGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAAcgcg                      >  1:79848/1‑100 (MQ=255)
                                                                           gCAAGCGGCGGAAGTGAACGCAATACTTGACGATCACGATGCACTATTTATCGACATGCGTAACAACTATGAGTATGAAGTGGGGGACTTTGAAAAcg                        >  2:3965696/1‑98 (MQ=255)
                                                                             aaGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCt                    >  1:1234015/1‑100 (MQ=255)
                                                                              aGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTc                   >  1:3251490/1‑100 (MQ=255)
                                                                              aGCGGCGGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTc                   <  2:2311312/100‑1 (MQ=255)
                                                                               gcggcgGAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCg                  <  2:2289184/100‑1 (MQ=255)
                                                                                  gcggAAGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGaaa               >  1:1631962/1‑100 (MQ=255)
                                                                                  gcggAAATGAACGCCATACTTGACGATCCCGATGCACTATTTATCGGCATGCGTAACCACTATGAGAATGAGGGGCGGCGCTTTGAATACGTGCCCGaa                >  2:3276599/1‑99 (MQ=255)
                                                                                       aGTGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTTAAAACGCGCTCGAAATTCCg          <  2:2007428/100‑1 (MQ=255)
                                                                                         tGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGc        <  1:2173599/100‑1 (MQ=255)
                                                                                         tGAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGc        <  2:79846/100‑1 (MQ=255)
                                                                                         tGAACGCCATACTTGACGATCCCGATGCACGATTTATCGACATTCGTAACCACTATGAGTATGAAGGGGGGCACTTTGAAAACGCGCTCGAAATTCCGGc        <  2:3436197/100‑1 (MQ=255)
                                                                                          gAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa       <  2:1026946/100‑1 (MQ=255)
                                                                                          gAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTCACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa       >  2:1725857/1‑100 (MQ=255)
                                                                                          gAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa       >  1:2074396/1‑100 (MQ=255)
                                                                                          gAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa       <  2:2210120/100‑1 (MQ=255)
                                                                                          gAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa       >  2:1680979/1‑100 (MQ=255)
                                                                                          gAACGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCa       >  1:1477651/1‑100 (MQ=255)
                                                                                           aaCGCCATGCTTGACGATCCCGATGCACTATGTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGACAACGCGCTCGAAATTCCGGCAg      >  2:3408225/1‑100 (MQ=255)
                                                                                           aaCGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAg      <  2:742856/100‑1 (MQ=255)
                                                                                           aaCGCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAg      <  2:1664728/100‑1 (MQ=255)
                                                                                            aCGCCATGCTTGACGATCCCGATGCACTATAAATCGACATCCCCAAACCATATAAGGATGAATTGGGGCAGTTTGACCAAGCGAACGAAAtgc            >  2:2057857/1‑91 (MQ=255)
                                                                                             cGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGAt    >  2:4132853/1‑100 (MQ=255)
                                                                                              gCCATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATa   <  1:3745295/100‑1 (MQ=255)
                                                                                               ccATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATAc  >  2:1339405/1‑100 (MQ=255)
                                                                                               ccATACTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATAc  >  1:1023793/1‑100 (MQ=255)
                                                                                                   |                                                                                               
GAAGGTACGCGATCGCATTGTTGCCGACGGTATTGACGATCCTCACTTTGATGCCAGCAATGTTGGTGAGTATCTGCAAGCGGCGGAAGTGAACGCCATGCTTGACGATCCCGATGCACTATTTATCGACATGCGTAACCACTATGAGTATGAAGTGGGGCACTTTGAAAACGCGCTCGAAATTCCGGCAGATAC  >  NC_000913/1117139‑1117333

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