Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,870,765 +G coding (476/1764 nt) dgoA ← 2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,870,7641.G99.3% 1354.4 / 0.1 431coding (477/1764 nt)dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (209/218);  total (211/220)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.58e-01

Only 100 of 477 total aligned reads displayed.

GCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATG  >  NC_000913/3870672‑3870851
                                                                                               |                                                                                       
gcgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCt                                                                                          <  2:5169157/93‑1 (MQ=255)
  ggTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCCTGCACAACCTGCGTCTCTCCACTGCGCCAGGTTTTCTg                                                                                        >  2:405476/1‑93 (MQ=255)
   gTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggc                                                                                      <  1:6044183/93‑3 (MQ=255)
    tGCCCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggcg                                                                                     <  1:4367093/93‑4 (MQ=255)
       gcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCa                                                                                   >  1:5404004/1‑93 (MQ=255)
       gcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCa                                                                                   <  2:1838649/93‑1 (MQ=255)
        ctctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAc                                                                                  <  1:1894681/93‑1 (MQ=255)
          ctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                >  2:2163707/1‑93 (MQ=255)
           tACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTGTTCTGGCGGCAcgcc                                                                               >  2:7026542/1‑93 (MQ=255)
             cGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgccgc                                                                             <  1:7620488/93‑1 (MQ=255)
               gATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCa                                                                           >  2:1983580/1‑93 (MQ=255)
               gATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCa                                                                           <  2:5043410/93‑1 (MQ=255)
                 ttGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAAc                                                                         >  2:2949035/1‑93 (MQ=255)
                 ttGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAAc                                                                         >  2:6905982/1‑93 (MQ=255)
                     ccGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCa                                                                     >  2:4457825/1‑93 (MQ=255)
                     ccGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCa                                                                     <  2:1241437/93‑1 (MQ=255)
                      cGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaa                                                                    >  1:936369/1‑93 (MQ=255)
                       ggCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaaa                                                                   >  2:3849012/1‑93 (MQ=255)
                          gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                >  2:4657976/1‑93 (MQ=255)
                          gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                <  2:2733091/93‑1 (MQ=255)
                          gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                <  1:6220129/93‑1 (MQ=255)
                             aTAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                                             >  2:3017337/1‑93 (MQ=255)
                              tAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCg                                                            >  2:1316630/1‑93 (MQ=255)
                              tAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCg                                                            >  1:6929630/1‑93 (MQ=255)
                               agagATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGa                                                           <  1:4490400/93‑1 (MQ=255)
                               agagATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGa                                                           <  1:6711781/93‑1 (MQ=255)
                                gagaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGAt                                                          <  1:2754441/93‑1 (MQ=255)
                                gagaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGAt                                                          <  1:2494304/93‑1 (MQ=255)
                                aagaTCGCTACCTAAGCCCGCCCATGCACAACCTGCGTCTATCCACTGCCACAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGAt                                                          <  2:7933727/92‑1 (MQ=255)
                                  gaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                                        >  1:2119585/1‑93 (MQ=255)
                                   aTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTc                                                       >  2:2367493/1‑93 (MQ=255)
                                   aTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTc                                                       <  1:7883766/93‑1 (MQ=255)
                                   aTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTc                                                       <  1:7077741/93‑1 (MQ=255)
                                    tCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                      <  2:4207164/93‑1 (MQ=255)
                                    tCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                      >  1:5891608/1‑93 (MQ=255)
                                      gCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                    <  1:718104/93‑1 (MQ=255)
                                         gCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATgg                                                 <  2:2282497/93‑1 (MQ=255)
                                           cTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCa                                               >  2:5781315/1‑93 (MQ=255)
                                           cTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCa                                               >  1:3856422/1‑93 (MQ=255)
                                            tAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCTAAGACTGCGATGTCCGATGGCat                                              >  2:5415857/1‑92 (MQ=255)
                                            tAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCaa                                              >  2:4852033/1‑93 (MQ=255)
                                             aaGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAAt                                             <  1:7267898/93‑1 (MQ=255)
                                              aGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATa                                            <  2:6242743/93‑1 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  1:1313871/93‑1 (MQ=255)
                                                cccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAcc                                          <  1:1383306/93‑1 (MQ=255)
                                                  cgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGc                                        <  2:6548145/93‑1 (MQ=255)
                                                   gcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                                       >  2:7470232/1‑93 (MQ=255)
                                                      cccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCttt                                     >  2:1940458/1‑93 (MQ=255)
                                                      cccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCttt                                     <  1:1609069/93‑1 (MQ=255)
                                                        cTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                                   >  2:4366827/1‑93 (MQ=255)
                                                         tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  >  1:4807564/1‑93 (MQ=255)
                                                         tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  >  1:5309445/1‑93 (MQ=255)
                                                         tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  >  1:3316990/1‑93 (MQ=255)
                                                         tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  >  1:1502191/1‑93 (MQ=255)
                                                         cgtgcaACCGGCGGTTATTCACGGTGCCATGGTTTGTGGCGGCACGCCGCCAACGGCAAAGACTGCGATGGCCGATGGCAATACCGCTTTTaa                                  <  1:2000532/89‑1 (MQ=255)
                                                          gCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAAc                                 <  1:3417245/93‑1 (MQ=255)
                                                             caACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                              >  1:5026513/1‑93 (MQ=255)
                                                              aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                             >  1:5175204/1‑93 (MQ=255)
                                                              aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                             <  2:1741035/93‑1 (MQ=255)
                                                              aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                             <  2:7347027/93‑1 (MQ=255)
                                                               aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                            <  2:3362282/93‑1 (MQ=255)
                                                               aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                            <  1:4644255/93‑1 (MQ=255)
                                                                  tGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                         >  1:6654653/1‑93 (MQ=255)
                                                                  tGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                         >  1:5471312/1‑93 (MQ=255)
                                                                   gCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATg                        <  2:1086184/93‑1 (MQ=255)
                                                                    cGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGt                       <  1:4169359/93‑1 (MQ=255)
                                                                     gTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTa                      <  2:6457621/93‑1 (MQ=255)
                                                                     gTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTa                      >  2:1513529/1‑93 (MQ=255)
                                                                      tCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAt                     <  2:6813602/93‑1 (MQ=255)
                                                                      tCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAt                     >  2:877360/1‑93 (MQ=255)
                                                                      tCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAt                     <  2:3539294/93‑1 (MQ=255)
                                                                      tCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAt                     >  1:7492640/1‑93 (MQ=255)
                                                                       cTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  2:5996787/93‑1 (MQ=255)
                                                                       cTAACCACGGAGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  2:1431327/93‑1 (MQ=255)
                                                                        tATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTg                   <  1:6444102/93‑1 (MQ=255)
                                                                         aTCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGc                  <  1:6528931/93‑1 (MQ=255)
                                                                         aTCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGc                  >  2:5502308/1‑93 (MQ=255)
                                                                          tCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg                 >  1:3549747/1‑93 (MQ=255)
                                                                          tCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg                 >  1:5574545/1‑93 (MQ=255)
                                                                            cACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCggc               >  2:2476452/1‑92 (MQ=255)
                                                                            cACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGa               <  1:768890/93‑1 (MQ=255)
                                                                            cACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGa               >  1:4242808/1‑93 (MQ=255)
                                                                             aCTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAc              >  2:7206379/1‑93 (MQ=255)
                                                                              cTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAcc             >  2:3193708/1‑93 (MQ=255)
                                                                                  gcaaGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  2:3278430/90‑1 (MQ=255)
                                                                                  gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  1:1224124/93‑1 (MQ=255)
                                                                                  gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  2:3627658/93‑1 (MQ=255)
                                                                                    cAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGc       >  1:4072325/1‑93 (MQ=255)
                                                                                     aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      >  2:7804647/1‑93 (MQ=255)
                                                                                     aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      <  1:873741/93‑1 (MQ=255)
                                                                                       ggTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCg     >  2:7599837/1‑93 (MQ=255)
                                                                                       ggTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCg     >  1:1004271/1‑93 (MQ=255)
                                                                                       ggTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCg     <  2:6149023/93‑1 (MQ=255)
                                                                                        gTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    >  1:1766137/1‑93 (MQ=255)
                                                                                        gTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    >  1:979390/1‑93 (MQ=255)
                                                                                         ttttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt   <  2:5558087/93‑1 (MQ=255)
                                                                                         ttttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt   <  2:5351814/93‑1 (MQ=255)
                                                                                         ttttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt   >  2:4964136/1‑93 (MQ=255)
                                                                                          tttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATg  >  1:3634925/1‑93 (MQ=255)
                                                                                          tttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATg  <  1:3190211/93‑1 (MQ=255)
                                                                                               |                                                                                       
GCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATG  >  NC_000913/3870672‑3870851

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