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.7% 1235.0 / ‑4.5 381coding (477/1764 nt)dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (176/204);  total (177/204)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 423 total aligned reads displayed.

GCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATG  >  NC_000913/3870672‑3870851
                                                                                             |                                                                                       
gcgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCt                                                                                          <  2:1214802/93‑1 (MQ=255)
gcgGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCt                                                                                          <  1:498315/93‑1 (MQ=255)
    tGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggcg                                                                                     >  2:4318146/1‑90 (MQ=255)
    tGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggcg                                                                                     >  1:6798234/1‑90 (MQ=255)
     gcgcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTggcgt                                                                                    >  1:9705973/1‑89 (MQ=255)
       gcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCa                                                                                   <  1:195721/93‑1 (MQ=255)
          ctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGGCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                >  1:5034865/1‑93 (MQ=255)
          ctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                >  1:1339606/1‑93 (MQ=255)
          ctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                >  1:8701605/1‑93 (MQ=255)
          ctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgc                                                                                <  1:6320118/93‑1 (MQ=255)
           tACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgcc                                                                               <  1:9897756/93‑1 (MQ=255)
             cGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCAcgccgc                                                                             <  2:9643513/93‑1 (MQ=255)
              ggATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACgccgcc                                                                            >  2:6518470/1‑93 (MQ=255)
               gATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCa                                                                           >  2:8909983/1‑93 (MQ=255)
                aTTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCaa                                                                          <  1:8225498/93‑1 (MQ=255)
                       ggCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaaa                                                                   >  1:5883011/1‑93 (MQ=255)
                       ggCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaaa                                                                   <  2:7768204/93‑1 (MQ=255)
                        gcgcGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAg                                                                  <  1:2696626/93‑1 (MQ=255)
                        gcgcGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAg                                                                  >  1:4326208/1‑93 (MQ=255)
                         cgcgATAGAGATCGCTGCCTAAGCCCGCTCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGa                                                                 <  2:940752/93‑1 (MQ=255)
                         cgcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGa                                                                 <  1:7631280/93‑1 (MQ=255)
                         cgcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGa                                                                 <  1:5751108/93‑1 (MQ=255)
                          gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                <  1:4127723/93‑1 (MQ=255)
                          gcgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGAc                                                                >  1:7361929/1‑93 (MQ=255)
                           ggATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACt                                                               <  2:7452936/92‑1 (MQ=255)
                           cgATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACt                                                               <  2:6734773/93‑1 (MQ=255)
                             aTAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                                             <  1:6146796/93‑1 (MQ=255)
                             aTAGAGATCGCTGCCTAAGCCCGCCCATGCACAACCTGCGTCTATCCACGGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGc                                                             <  2:3280175/93‑1 (MQ=255)
                              tAGAGATCGCTGCCTAAGCCCGCCCCTGCCCAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCg                                                            <  2:2549481/93‑1 (MQ=255)
                              tAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCg                                                            >  2:1522392/1‑93 (MQ=255)
                              tAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCg                                                            >  2:265911/1‑93 (MQ=255)
                               agagATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGa                                                           <  2:4689038/93‑1 (MQ=255)
                                gagaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGAt                                                          >  2:8142684/1‑93 (MQ=255)
                                gagaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGAt                                                          >  2:2194241/1‑93 (MQ=255)
                                 agaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATg                                                         <  2:1410053/93‑1 (MQ=255)
                                 agaTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGACTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATg                                                         <  1:7024358/93‑1 (MQ=255)
                                  gaTCGCTGCCTAAGCCCACCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGTCAAAGACTGCGATGt                                                        <  1:5400940/93‑1 (MQ=255)
                                   aTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTc                                                       >  2:7792236/1‑93 (MQ=255)
                                   aTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTc                                                       <  2:1688088/93‑1 (MQ=255)
                                   aTCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTc                                                       >  1:9412565/1‑93 (MQ=255)
                                      gCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                    <  1:9812879/93‑1 (MQ=255)
                                      gCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                    <  1:4753937/93‑1 (MQ=255)
                                      gCTGCCGATGGCCGCCCCTGCACAGGCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                    <  1:5796629/93‑1 (MQ=255)
                                       cTGCCTCAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCAAACGGCAAAGACTGCGATGTCCGAt                                                   <  2:3944312/93‑1 (MQ=255)
                                        tGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                  >  1:3693974/1‑93 (MQ=255)
                                        tGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATg                                                  >  2:7300646/1‑93 (MQ=255)
                                           cTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCa                                               <  2:2855775/93‑1 (MQ=255)
                                            tAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCaa                                              >  1:9583125/1‑93 (MQ=255)
                                            tAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCaa                                              <  1:232198/93‑1 (MQ=255)
                                            tAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCaa                                              >  1:1768703/1‑93 (MQ=255)
                                            tAACCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCaa                                              <  2:2027934/93‑1 (MQ=255)
                                             aaGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAAt                                             <  1:8534773/93‑1 (MQ=255)
                                             aaGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAAt                                             >  1:6009298/1‑93 (MQ=255)
                                              aGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATa                                            >  1:9475178/1‑93 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  2:9482950/93‑1 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  1:4410840/93‑1 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           >  2:885140/1‑93 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  2:8253116/93‑1 (MQ=255)
                                               gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAc                                           <  1:2950854/93‑1 (MQ=255)
                                                cccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAcc                                          <  1:2433696/93‑1 (MQ=255)
                                                cccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATAcc                                          <  2:7035338/93‑1 (MQ=255)
                                                 ccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCg                                         >  1:3339139/1‑93 (MQ=255)
                                                   cccccTGCACAACCTGCTTTTATCCACTGCGCAAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCt                                       <  2:3708064/92‑1 (MQ=255)
                                                    ccccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtt                                      <  1:2574653/93‑1 (MQ=255)
                                                     cccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCttt                                     >  2:5483841/1‑93 (MQ=255)
                                                      ccTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCtttt                                    <  2:7377470/93‑1 (MQ=255)
                                                       cTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                                   <  2:3479873/93‑1 (MQ=255)
                                                        tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  >  2:3894926/1‑93 (MQ=255)
                                                        tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  <  1:8090665/93‑1 (MQ=255)
                                                        tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  <  2:400617/93‑1 (MQ=255)
                                                        tGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTaa                                  <  1:2231551/93‑1 (MQ=255)
                                                         gCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAAc                                 <  1:86147/93‑1 (MQ=255)
                                                           acaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGc                               <  2:9212503/93‑1 (MQ=255)
                                                           acaACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGc                               <  2:6394742/93‑1 (MQ=255)
                                                            caACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                              >  2:4496941/1‑93 (MQ=255)
                                                                cTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGa                          <  2:4983136/93‑1 (MQ=255)
                                                                cTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGa                          <  1:6558957/93‑1 (MQ=255)
                                                                  gCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATg                        <  2:5921887/93‑1 (MQ=255)
                                                                    gTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTa                      <  2:1845495/93‑1 (MQ=255)
                                                                     tCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTAt                     >  2:1329890/1‑93 (MQ=255)
                                                                       tATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTg                   <  2:2473493/93‑1 (MQ=255)
                                                                        atacacaGCGCAAGGTTTTCTGGCTTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACGGCTTTTAACGCTTTGATGTATTGc                  <  2:5169878/90‑1 (MQ=255)
                                                                         tCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCg                 >  1:1119072/1‑93 (MQ=255)
                                                                          ccACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCgg                >  1:834089/1‑93 (MQ=255)
                                                                            ccTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGAc              <  2:7204350/92‑1 (MQ=255)
                                                                              tGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCa            <  2:5639588/93‑1 (MQ=255)
                                                                                cgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaa          <  2:313526/93‑1 (MQ=255)
                                                                                cgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaa          >  1:3579223/1‑93 (MQ=255)
                                                                                cgcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaa          <  2:8587954/93‑1 (MQ=255)
                                                                                 gcCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCaaaa         <  2:6182876/93‑1 (MQ=255)
                                                                                  ccAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAg        <  1:10003811/93‑1 (MQ=255)
                                                                                  ccAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAg        >  2:2991321/1‑93 (MQ=255)
                                                                                  ccAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAg        <  2:4144510/93‑1 (MQ=255)
                                                                                   cAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGc       <  1:8996904/93‑1 (MQ=255)
                                                                                    aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      >  1:8785660/1‑93 (MQ=255)
                                                                                    aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      >  1:5224932/1‑93 (MQ=255)
                                                                                    aGGTTTCCTGGCATCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGcc      <  1:3955616/93‑1 (MQ=255)
                                                                                      gTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    <  2:3101788/93‑1 (MQ=255)
                                                                                      gTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGa    <  1:7233013/93‑1 (MQ=255)
                                                                                        tttCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATg  <  1:4589577/93‑1 (MQ=255)
                                                                                             |                                                                                       
GCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATG  >  NC_000913/3870672‑3870851

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