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.5% 710.9 / ‑3.8 217coding (477/1764 nt)dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (113/103);  total (114/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

Only 100 of 242 total aligned reads displayed.

GCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATG  >  NC_000913/3870672‑3870851
                                                                                              |                                                                                       
gcgGTGCGCTCTACGGATTTCCGGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCt                                                                                          <  2:3345831/93‑1 (MQ=255)
gcgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCt                                                                                          >  1:3319356/1‑93 (MQ=255)
 cgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTg                                                                                        <  2:363598/93‑1 (MQ=255)
 cgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGAACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTg                                                                                        <  1:1943518/93‑1 (MQ=255)
  ggTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTgg                                                                                       >  1:4094075/1‑92 (MQ=255)
    ggcgcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggcg                                                                                     <  2:524380/92‑4 (MQ=255)
        ctctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAc                                                                                  <  1:4320285/93‑1 (MQ=255)
         tctACGTATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGCTTTCTGGCGTCAcg                                                                                 >  2:2143354/1‑93 (MQ=255)
          ctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                >  1:231680/1‑93 (MQ=255)
           tACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgcc                                                                               >  2:167213/1‑93 (MQ=255)
             cGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgccgc                                                                             >  1:2044719/1‑93 (MQ=255)
              ggATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACgccgcc                                                                            >  1:1524938/1‑93 (MQ=255)
               gATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCa                                                                           >  2:1458883/1‑93 (MQ=255)
                 ttGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAAc                                                                         <  1:3530963/93‑1 (MQ=255)
                  tGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACg                                                                        <  2:1454322/93‑1 (MQ=255)
                  tGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACg                                                                        >  2:2263878/1‑93 (MQ=255)
                   gCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACgg                                                                       <  2:1695674/93‑1 (MQ=255)
                    cccGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                      <  1:1276844/93‑1 (MQ=255)
                      cGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaa                                                                    >  1:643731/1‑93 (MQ=255)
                          gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGTGCCAGGTTTTCTGGCGTCACGCCGCCATCGGCAAAGAc                                                                >  2:4366808/1‑93 (MQ=255)
                          gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                <  2:1992343/93‑1 (MQ=255)
                          gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                <  1:2982092/93‑1 (MQ=255)
                          gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                <  2:3558713/93‑1 (MQ=255)
                          gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                <  1:212244/93‑1 (MQ=255)
                            gATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTg                                                              <  1:2147419/93‑1 (MQ=255)
                            gATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTg                                                              <  1:4711703/93‑1 (MQ=255)
                            gATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTg                                                              >  2:4149410/1‑93 (MQ=255)
                                  gaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                                        >  2:458892/1‑93 (MQ=255)
                                  gaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                                        >  1:1403236/1‑93 (MQ=255)
                                  gaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                                        >  2:3902794/1‑93 (MQ=255)
                                  gaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGt                                                        >  1:3926547/1‑93 (MQ=255)
                                   aTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTc                                                       <  2:1251622/93‑1 (MQ=255)
                                   aTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTc                                                       >  2:1082006/1‑93 (MQ=255)
                                     cGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCg                                                     >  2:4736271/1‑93 (MQ=255)
                                      gCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                    <  2:4105759/93‑1 (MQ=255)
                                       cTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGAt                                                   >  1:1110819/1‑93 (MQ=255)
                                       cTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGAt                                                   <  1:2650321/93‑1 (MQ=255)
                                       cTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGAt                                                   <  1:3016399/93‑1 (MQ=255)
                                       cTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGAt                                                   <  1:3132497/93‑1 (MQ=255)
                                        tGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                  >  1:396131/1‑93 (MQ=255)
                                             aaGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAAt                                             >  1:552924/1‑93 (MQ=255)
                                             aaGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAAt                                             >  1:613498/1‑93 (MQ=255)
                                              aGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCCACGACTGCGATGTCCGCTGGCAATa                                            >  1:258675/1‑93 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           >  1:2437937/1‑93 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  2:3090030/93‑1 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  2:2284717/93‑1 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  2:193703/93‑1 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTATATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  2:2590859/93‑1 (MQ=255)
                                                cccgcccCTGCACAACCTGCGTCTATCCATAGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAcc                                          <  2:3485253/93‑1 (MQ=255)
                                                  cgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGc                                        <  2:2860162/93‑1 (MQ=255)
                                                  cgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGc                                        <  1:3356872/93‑1 (MQ=255)
                                                   gcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                                       <  2:917338/93‑1 (MQ=255)
                                                   gcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                                       >  2:1219192/1‑93 (MQ=255)
                                                   cccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                                       <  1:37299/92‑1 (MQ=255)
                                                   cccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                                       <  1:429293/92‑1 (MQ=255)
                                                    ccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtt                                      >  2:1336377/1‑93 (MQ=255)
                                                    ccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtt                                      >  2:4377959/1‑93 (MQ=255)
                                                    ccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtt                                      >  1:3082217/1‑93 (MQ=255)
                                                     cccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCttt                                     >  2:446486/1‑93 (MQ=255)
                                                      ccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCGACGGCAAAGACTGCGATGTCCGATGGCAATACTGCtttt                                    >  1:2550600/1‑93 (MQ=255)
                                                      ccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtttt                                    >  2:4037101/1‑93 (MQ=255)
                                                      ccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtttt                                    >  1:3635970/1‑93 (MQ=255)
                                                       ccgCAAAACCTGCGTCTAGCCACTTCGCCAGGTTTTCTGGCGTCACCCCGCCAACGGCAAAGACTTCGATGTCCGATTGCAATACCGCTTTTa                                   <  1:1671766/91‑1 (MQ=255)
                                                        tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAGTGCGATGTCCGATGGCAATACCGCTTTTaa                                  >  2:4637896/1‑93 (MQ=255)
                                                           acacCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGc                               <  2:1579574/93‑1 (MQ=255)
                                                           acaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGc                               >  1:3193597/1‑93 (MQ=255)
                                                           acaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGc                               <  2:2323955/93‑1 (MQ=255)
                                                            caACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                              >  1:1808783/1‑93 (MQ=255)
                                                             aaCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                             <  1:1651817/93‑1 (MQ=255)
                                                              aCCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCttg                            >  1:3712953/1‑92 (MQ=255)
                                                                cTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCTCGCCGCTAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGa                          <  2:453295/93‑1 (MQ=255)
                                                                 tGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                         >  1:2769585/1‑93 (MQ=255)
                                                                   cGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGTCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGt                       <  1:4005634/93‑1 (MQ=255)
                                                                     tCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAt                     >  2:2984552/1‑93 (MQ=255)
                                                                      cTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  1:154770/93‑1 (MQ=255)
                                                                      cTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  1:3443098/93‑1 (MQ=255)
                                                                      cTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAtt                    <  2:1827924/93‑1 (MQ=255)
                                                                       tATCCTCTGCGCCAGATTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTg                   >  2:393583/1‑93 (MQ=255)
                                                                        aTCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGc                  >  2:3646787/1‑93 (MQ=255)
                                                                         tCCACTGCGCCAGGTTTTCTTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg                 >  1:444754/1‑93 (MQ=255)
                                                                         tCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg                 >  1:1028743/1‑93 (MQ=255)
                                                                         tCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg                 <  2:815599/93‑1 (MQ=255)
                                                                          ccACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCgg                >  1:135813/1‑93 (MQ=255)
                                                                           cACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGa               <  1:307894/93‑1 (MQ=255)
                                                                            aCTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAc              >  1:2570195/1‑93 (MQ=255)
                                                                            aCTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAc              >  1:4179338/1‑93 (MQ=255)
                                                                             cTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAcc             <  2:2863808/93‑1 (MQ=255)
                                                                               gcgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaa           >  2:2427503/1‑93 (MQ=255)
                                                                                cgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaa          <  1:3247814/93‑1 (MQ=255)
                                                                                 gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  2:3290102/93‑1 (MQ=255)
                                                                                  ccAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAg        <  2:474924/93‑1 (MQ=255)
                                                                                    aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTACCGCTTTGATGTATTGCGGACCAAAAGcc      >  1:4368963/1‑93 (MQ=255)
                                                                                    aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      <  2:3065011/93‑1 (MQ=255)
                                                                                    aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      <  2:2617479/93‑1 (MQ=255)
                                                                                    aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      >  1:464875/1‑93 (MQ=255)
                                                                                    aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      <  2:1724667/93‑1 (MQ=255)
                                                                                     ggTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCg     >  2:1657426/1‑93 (MQ=255)
                                                                                      gTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    <  1:2281473/93‑1 (MQ=255)
                                                                                      gTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    >  2:3426300/1‑93 (MQ=255)
                                                                                        tttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATg  <  2:138598/93‑1 (MQ=255)
                                                                                              |                                                                                       
GCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATG  >  NC_000913/3870672‑3870851

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