Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,345 A→G 15.0% intergenic (‑141/+167) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,177,3450AG15.0% 337.2 / 46.0 139intergenic (‑141/+167)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base A (50/68);  new base G (15/6);  total (65/74)
Fisher's exact test for biased strand distribution p-value = 1.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 140 total aligned reads displayed.

ATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGC  >  NC_000913/2177248‑2177443
                                                                                                 |                                                                                                  
 ttattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                                                                                                 <  1:1665388/100‑1 (MQ=255)
    tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCTATCGGGAGGCCGt                                                                                               <  1:297520/99‑1 (MQ=255)
    tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                              <  1:1844783/100‑1 (MQ=255)
    tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                              <  1:1844512/100‑1 (MQ=255)
    tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                              <  2:3349920/100‑1 (MQ=255)
     ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGttt                                                                                             <  1:2202965/100‑3 (MQ=38)
     ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCTATCGGGAGGCCGtt                                                                                              <  1:1326770/99‑2 (MQ=255)
     ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                              <  2:1502849/99‑1 (MQ=255)
     ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                                                                             <  2:3151866/100‑1 (MQ=255)
     ttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                                                                             <  1:3577192/100‑1 (MQ=255)
      tGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTTAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                              <  1:2001615/98‑1 (MQ=255)
      tGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGc                                                                                            >  2:3384311/1‑100 (MQ=255)
       gTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCa                                                                                           >  1:3842488/1‑100 (MQ=255)
       gTTTTTAAAGATCGAGATAAAGATCACAAAAACGACAATATGAAAATTATTCGAGGAAAAAAAAACAAAAAAACGGCTTCCCTATCGGGAGGCCGtt                                                                                              <  1:1178793/97‑2 (MQ=19)
       gTTTTTAAAAATCGAGATAACGATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAATGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                                                                             <  1:739222/98‑1 (MQ=255)
        tttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCAGATAGGGAAGCCGTAGCaa                                                                                          >  1:1950877/1‑100 (MQ=255)
          tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttg                                                                                        <  2:125967/100‑8 (MQ=37)
          tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAg                                                                                        >  2:987140/1‑100 (MQ=255)
          tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGAAAGGGAAGACGTAGCAAAg                                                                                        >  1:995517/1‑100 (MQ=255)
           ttAAATATCGGGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttttgc                                                                                       <  2:3028176/100‑9 (MQ=25)
            tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                                                                                      <  2:3434068/100‑1 (MQ=255)
            tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                                                                                      >  1:279017/1‑100 (MQ=255)
            tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                                                                                      <  2:1609025/100‑1 (MQ=255)
            tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTg                                                                                      <  1:3384293/100‑1 (MQ=255)
              natatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcg                                                                                    >  2:1140142/2‑100 (MQ=255)
                  tCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                                                                                >  2:3896238/1‑100 (MQ=255)
                        tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                                          <  2:215474/100‑1 (MQ=255)
                        tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                                          >  2:2086523/1‑100 (MQ=255)
                        tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                                          <  1:2065825/100‑1 (MQ=255)
                         aaCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt                                                                         >  2:2538191/1‑100 (MQ=255)
                          aCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                                        <  2:3906445/100‑1 (MQ=255)
                           cGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                                                                       >  2:4077612/1‑100 (MQ=255)
                             aTCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCg                                                                     <  1:2654011/100‑1 (MQ=255)
                               cacaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa                                                                   >  2:554160/1‑100 (MQ=255)
                               cacaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa                                                                   <  2:3533465/100‑1 (MQ=255)
                                acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                                                                  >  2:1947408/1‑100 (MQ=255)
                                  aaaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAACACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                                 >  2:575678/1‑99 (MQ=255)
                                    aaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                              <  1:3910984/100‑1 (MQ=255)
                                    aaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg                                                              >  1:2953027/1‑100 (MQ=255)
                                      aCGACAATATTAAAATTATTCGAGGAGTGACAGGCAAAAGACCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtg                                                            <  1:3721117/100‑1 (MQ=255)
                                      aCGACAATATGCAAATTATTCGAGGAGAAAAAGAAAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgt                                                             <  1:3674589/99‑1 (MQ=38)
                                      aCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtg                                                            <  2:4135154/100‑1 (MQ=255)
                                       cGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtgt                                                           >  2:1851073/1‑100 (MQ=255)
                                          cAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                                        >  2:2108195/1‑100 (MQ=255)
                                          cAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                                        <  2:4317316/100‑1 (MQ=255)
                                          cAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                                        <  1:4327082/100‑1 (MQ=255)
                                           aaTATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc                                                       <  2:2012332/100‑1 (MQ=255)
                                           aaTATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc                                                       <  2:3727706/100‑1 (MQ=255)
                                              atGAAAATTATTCTAGGAGTGAAAGGCAAAAAAAGGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGGGTaaaga                                                      >  1:2066743/1‑96 (MQ=39)
                                              atGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa                                                    >  2:1667292/1‑100 (MQ=255)
                                              atGAAAATTATTCGAGGAGGGAAAGGAAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa                                                    <  1:313776/100‑1 (MQ=255)
                                               tGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACagat                                                   >  2:270130/1‑97 (MQ=255)
                                                  aaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                                >  2:3714904/1‑100 (MQ=255)
                                                  aaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                                <  2:3452396/100‑1 (MQ=255)
                                                    attattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGc                                              <  2:4398633/100‑1 (MQ=255)
                                                    attattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGc                                              >  2:4088554/1‑100 (MQ=255)
                                                      tattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACATGGCAg                                            >  1:1783267/1‑100 (MQ=255)
                                                        ttCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                                          <  2:415633/100‑1 (MQ=255)
                                                        ttCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                                          <  1:2499263/100‑1 (MQ=255)
                                                        ttCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                                          <  1:13752/100‑1 (MQ=255)
                                                          cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                                        >  1:2266226/1‑100 (MQ=255)
                                                          cGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                                        >  2:2495892/1‑100 (MQ=255)
                                                           gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTa                                       >  2:2370704/1‑100 (MQ=255)
                                                                gTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAg                                  >  2:3365252/1‑100 (MQ=255)
                                                                gTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAg                                  <  2:3252661/100‑1 (MQ=255)
                                                                 tGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAGGTGCGCGTGTTTTTATGCCGGAGGCGGTGTAAACATACGTGGCAGTGCTATGTAgg                                 <  1:4213518/100‑1 (MQ=255)
                                                                   acaGGCAAAGAACCGGTCTCCCGATAGGGAAGCCGTAGCAAAATGCGTGTGTTTTTGTGCCGGATGCGGTGTAAACATACGTGGCAGTGTTATGTa                                   <  1:39970/94‑1 (MQ=255)
                                                                      ggCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGa                            <  1:1997295/100‑1 (MQ=255)
                                                                       gCAAAAAAACGGCCTCTCGATAGGGAAGCTGTAGCAAAGTGAGCGTGTTGTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATg                             <  1:2749236/98‑1 (MQ=255)
                                                                       gCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCGGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGAt                           >  1:2813366/1‑100 (MQ=255)
                                                                       gCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGAt                           <  1:2217189/100‑1 (MQ=255)
                                                                        cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa                          <  2:763775/100‑1 (MQ=255)
                                                                            gctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc                      >  2:4308835/4‑100 (MQ=38)
                                                                            gctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc                      >  1:1570684/4‑100 (MQ=38)
                                                                            gctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc                      >  2:469838/4‑100 (MQ=38)
                                                                            aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAc                      <  2:619309/100‑1 (MQ=255)
                                                                             ctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCTCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                     >  2:3102735/3‑100 (MQ=38)
                                                                             ctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                     >  2:303171/3‑100 (MQ=38)
                                                                             ctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                     >  2:1411510/3‑100 (MQ=38)
                                                                             ctaCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                     >  2:594380/3‑100 (MQ=38)
                                                                              taCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTTCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCGGTGCTATGTAGTCGTGATAAGACg                     >  1:358037/2‑99 (MQ=37)
                                                                               aCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCt                   >  1:1188975/1‑100 (MQ=38)
                                                                               aCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCt                   >  2:1456753/1‑100 (MQ=38)
                                                                               aCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCt                   >  2:1985792/1‑100 (MQ=255)
                                                                                 ggCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTc                 >  2:880280/1‑100 (MQ=255)
                                                                                  gCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCa                >  2:2030678/1‑100 (MQ=255)
                                                                                  gCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCa                <  2:1842293/100‑1 (MQ=255)
                                                                                   ccTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAg               <  2:2176848/100‑1 (MQ=255)
                                                                                    cTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGc              >  1:4362312/1‑100 (MQ=255)
                                                                                     tCCCGCTAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGACGGATGCGGTGTAAACATACGTGGCCGTGCTATGTAGGCCTGATCAGACGATCCAg               >  1:3077531/1‑98 (MQ=255)
                                                                                     tCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCg             <  2:1534546/100‑1 (MQ=255)
                                                                                      cccGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGt            <  2:3029632/100‑1 (MQ=255)
                                                                                      cccGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGt            >  2:2669936/1‑100 (MQ=255)
                                                                                        cGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCg          >  2:3478211/1‑100 (MQ=255)
                                                                                        cGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCg          <  2:44332/100‑1 (MQ=255)
                                                                                          aTAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCa        <  2:4174188/100‑1 (MQ=255)
                                                                                            aGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTGAAGCGTcttat       <  2:3683325/99‑5 (MQ=38)
                                                                                             gggAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGTAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCa     <  2:1163354/100‑1 (MQ=255)
                                                                                             gggAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGTCGCTTAAGCGTCGCATCa     <  1:101350/100‑1 (MQ=255)
                                                                                             gggAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCa     >  1:522541/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
ATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGC  >  NC_000913/2177248‑2177443

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: