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.8% 1807.1 / ‑4.1 462coding (477/1764 nt)dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (198/263);  total (199/263)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 508 total aligned reads displayed.

TGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACG  >  NC_000913/3870669‑3870854
                                                                                                |                                                                                          
tgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCGAAGCCCGCCCCTGCACAACCTGCGTCTACCCACTGCGCCAGGTTTTCt                                                                                             <  2:12773851/96‑1 (MQ=255)
   gcgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGGCTATCGACTGCGCCAGGTTTTCTgtc                                                                                         >  1:6957703/1‑94 (MQ=255)
   gcgGTGCGATCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggc                                                                                         <  1:1267051/96‑3 (MQ=255)
    cgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggcg                                                                                        <  2:6143679/96‑4 (MQ=255)
     ggTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggcgt                                                                                       <  1:9237998/96‑5 (MQ=255)
     ggTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggcgt                                                                                       <  2:4639210/96‑5 (MQ=255)
        gctctctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAc                                                                                     <  1:6818025/93‑1 (MQ=255)
        gcgcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTTTCCACTGCGCCAGGGTTTCTGGCGTCAc                                                                                     >  2:7574573/1‑96 (MQ=255)
        gcgcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAc                                                                                     <  1:8961959/96‑1 (MQ=255)
        gcgcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAc                                                                                     <  2:13153094/96‑1 (MQ=255)
         cgcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcg                                                                                    >  1:4606824/1‑96 (MQ=255)
          gcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                   >  2:1195380/1‑96 (MQ=255)
          gcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                   >  1:5952740/1‑96 (MQ=255)
          gcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                   <  1:5001826/96‑1 (MQ=255)
          gcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                   <  2:4763132/96‑1 (MQ=255)
            tctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgccg                                                                                 >  1:5662075/1‑96 (MQ=255)
             ctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgccgc                                                                                >  2:9103/1‑96 (MQ=255)
                cGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCaa                                                                             >  1:13330512/1‑96 (MQ=255)
                  gATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACg                                                                           <  1:773603/96‑1 (MQ=255)
                  gATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACg                                                                           >  2:10533864/1‑96 (MQ=255)
                       cccGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaaa                                                                      <  2:10690515/96‑1 (MQ=255)
                       cccGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaaa                                                                      >  2:8647540/1‑96 (MQ=255)
                       cccGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaaa                                                                      <  2:8341926/96‑1 (MQ=255)
                        cgtgcgcGATAGAGACCGCGGCATAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAg                                                                     <  1:9901144/93‑1 (MQ=255)
                        ccGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAg                                                                     >  2:10129892/1‑96 (MQ=255)
                         cGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGa                                                                    >  2:7864357/1‑96 (MQ=255)
                          ggCGCTCTAGAGATCGCTCCCTAAGCCCTCCCCTGCACAACCTGCGTCTATCCCCTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                   <  2:9680827/96‑1 (MQ=255)
                          ggCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                   >  1:6343041/1‑96 (MQ=255)
                          ggCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                   <  2:7255258/96‑1 (MQ=255)
                           gcgcGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCCCGCCGCCAACGGCAAAGACt                                                                  <  1:11928412/96‑1 (MQ=255)
                            cgcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTg                                                                 <  1:4038079/96‑1 (MQ=255)
                            cgcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGAGACAGGTTTTCTGGCGTCACGCCGCTAACAGCAAAGACTg                                                                 >  1:784972/1‑96 (MQ=255)
                             gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                                                <  2:10845364/96‑1 (MQ=255)
                             gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                                                >  1:13219445/1‑96 (MQ=255)
                               gATAGAGAGCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGa                                                              <  1:12155976/96‑1 (MQ=255)
                                  agagATCGCCGCTAATGCCCGCCCCTGCACAGCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                                           <  2:5153404/96‑1 (MQ=255)
                                    agaTCGCTGCCTAGCCCCGCCCCTGCCCAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGCCTGCGATGTcc                                                         <  1:7355355/96‑1 (MQ=255)
                                    agaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                         >  1:14252585/1‑96 (MQ=255)
                                    agaACGCTGCCTAAGCCCGCCCCTGCACAACCTCCGTCTATACACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                         <  2:6478857/96‑1 (MQ=255)
                                       ccGCTGCCTAAGCCCGCCCCTGCAAAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGAt                                                      <  1:10033790/95‑1 (MQ=255)
                                        cGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                     <  2:1479548/96‑1 (MQ=255)
                                        cGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                     >  2:3181395/1‑96 (MQ=255)
                                        cGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                     >  1:337686/1‑96 (MQ=255)
                                        cGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                     <  1:4274647/96‑1 (MQ=255)
                                          cTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGc                                                   >  2:11649483/1‑96 (MQ=255)
                                          cTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGc                                                   <  1:13552483/96‑1 (MQ=255)
                                            gCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCaa                                                 >  2:11491198/1‑96 (MQ=255)
                                             cttAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAAt                                                <  2:1784050/94‑1 (MQ=255)
                                             ccTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAAt                                                >  2:1974226/1‑96 (MQ=255)
                                               tAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                              >  1:833039/1‑96 (MQ=255)
                                                aaGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAcc                                             <  2:3503550/96‑1 (MQ=255)
                                                aaGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAcc                                             >  1:526340/1‑96 (MQ=255)
                                                 ggcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCg                                            <  2:2582550/95‑1 (MQ=255)
                                                  tcccgcccTTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGc                                           <  1:2930630/95‑1 (MQ=255)
                                                  gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGCCTGCGATGTCCGATGGCAATACCGc                                           >  1:1365393/1‑96 (MQ=255)
                                                   cccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                                          >  2:1336758/1‑96 (MQ=255)
                                                     cgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCttt                                        <  1:10549291/96‑1 (MQ=255)
                                                      gcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATCCCGCtttt                                       <  1:2162859/96‑1 (MQ=255)
                                                      gcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtttt                                       <  1:1146938/96‑1 (MQ=255)
                                                       ccccTGCACAACCTGCGTCTATCCACTGGGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                                      >  2:13918066/1‑96 (MQ=255)
                                                       ccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                                      <  2:9478073/96‑1 (MQ=255)
                                                       ccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                                      >  2:11348543/1‑96 (MQ=255)
                                                       ccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                                      <  2:337686/96‑1 (MQ=255)
                                                         cccgCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAAc                                    <  1:5445331/93‑1 (MQ=255)
                                                         ccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAAc                                    <  2:12094721/96‑1 (MQ=255)
                                                         ccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAAc                                    <  1:13942893/96‑1 (MQ=255)
                                                            gCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAATACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                                 <  1:12940728/96‑1 (MQ=255)
                                                            gCACAACCTGAGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                                 <  1:9473998/96‑1 (MQ=255)
                                                             cacaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                                <  1:12319026/96‑1 (MQ=255)
                                                              acaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                               <  2:2284361/96‑1 (MQ=255)
                                                              acaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                               <  1:9673155/96‑1 (MQ=255)
                                                              acaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttt                               >  1:9946263/1‑96 (MQ=255)
                                                                aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGa                             >  1:13153702/1‑96 (MQ=255)
                                                                aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGa                             >  1:280620/1‑96 (MQ=255)
                                                                 aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                            <  2:14037290/96‑1 (MQ=255)
                                                                 aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                            >  2:870698/1‑96 (MQ=255)
                                                                  ccTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATg                           >  2:4440014/1‑96 (MQ=255)
                                                                   cTGCGTCTATCCACTGCGCGAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGt                          <  2:13604598/96‑1 (MQ=255)
                                                                   cTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGt                          <  2:5853603/96‑1 (MQ=255)
                                                                     gCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAt                        <  2:8914968/96‑1 (MQ=255)
                                                                      cGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                       <  2:3870687/96‑1 (MQ=255)
                                                                       gTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTg                      >  2:3730465/1‑96 (MQ=255)
                                                                         cTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg                    <  1:11268176/96‑1 (MQ=255)
                                                                             ccACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAcc                <  2:5522965/96‑1 (MQ=255)
                                                                               aCTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaa              <  1:8637863/96‑1 (MQ=255)
                                                                                cTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaa             >  2:808107/1‑96 (MQ=255)
                                                                                cTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaa             <  2:4148739/96‑1 (MQ=255)
                                                                                  gcgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAg           <  2:11040946/96‑1 (MQ=255)
                                                                                   cgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGc          <  2:11874846/96‑1 (MQ=255)
                                                                                   cgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGc          >  1:2639018/1‑96 (MQ=255)
                                                                                    gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc         <  1:1080387/96‑1 (MQ=255)
                                                                                    gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc         <  2:2784739/96‑1 (MQ=255)
                                                                                    gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc         >  1:14136658/1‑96 (MQ=255)
                                                                                    gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc         <  1:1651327/96‑1 (MQ=255)
                                                                                    gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc         >  2:4957838/1‑96 (MQ=255)
                                                                                     ccAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCg        <  1:3674759/96‑1 (MQ=255)
                                                                                       aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt      <  2:12333702/96‑1 (MQ=255)
                                                                                         gTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGa    >  1:4918731/1‑96 (MQ=255)
                                                                                         gTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGa    >  2:69417/1‑96 (MQ=255)
                                                                                          ttttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGAc   <  2:9958319/96‑1 (MQ=255)
                                                                                                |                                                                                          
TGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACG  >  NC_000913/3870669‑3870854

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: