Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,870,631 Δ1 bp coding (610/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,6300C.99.4% 2412.7 / 1.7 506coding (611/1764 nt)dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (276/226);  total (277/229)
Fisher's exact test for biased strand distribution p-value = 9.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 529 total aligned reads displayed.

GCCCCAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCC  >  NC_000913/3870536‑3870721
                                                                                              |                                                                                           
gCCCCAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGGGATTTTCATTGCAct                                                                                            >  1:5657660/1‑95 (MQ=255)
 ccccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActg                                                                                           >  2:6808329/1‑94 (MQ=255)
    cAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCct                                                                                       >  1:9098223/1‑96 (MQ=255)
     aaCCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCctc                                                                                      >  1:6893352/1‑96 (MQ=255)
     aaCCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGAAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCctc                                                                                      >  2:12786907/1‑96 (MQ=255)
       ccGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCtctc                                                                                    >  2:12986133/1‑96 (MQ=255)
        cGACCCCGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCg                                                                                   <  2:882728/96‑1 (MQ=255)
         gACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                                  <  1:9479279/96‑1 (MQ=255)
         gACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                                  <  2:3860484/96‑1 (MQ=255)
          aCCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGAt                                                                                 <  1:13550294/96‑1 (MQ=255)
          aCCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGGGGTAATTTTGGGGATTTTCATTGCACTGCCTCTCGAt                                                                                 >  2:5555860/1‑96 (MQ=255)
             aCGCCTTCATCGTTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGTTTTTCACTGCACTGCCTCTCGATACg                                                                              <  1:13223645/96‑1 (MQ=255)
             aCGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACg                                                                              <  1:7364595/96‑1 (MQ=255)
              cGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGc                                                                             <  2:1802950/96‑1 (MQ=255)
                ccTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCt                                                                           >  2:11621098/1‑96 (MQ=255)
                 cctCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCtt                                                                          >  1:11874846/3‑96 (MQ=255)
                 cTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCtt                                                                          >  1:2918442/1‑96 (MQ=255)
                 cTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCtt                                                                          <  2:7360694/96‑1 (MQ=255)
                 cTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCtt                                                                          <  2:4752540/96‑1 (MQ=255)
                  ttCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                                         <  1:12024882/96‑1 (MQ=255)
                  ttCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                                         >  2:8621715/1‑96 (MQ=255)
                    cATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAAc                                                                       >  2:10770505/1‑96 (MQ=255)
                       cGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGGTACGTGGGAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaaa                                                                    >  1:4599234/1‑96 (MQ=255)
                        ggTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                                                   <  2:13529785/96‑1 (MQ=255)
                        ggTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                                                   >  1:1711869/1‑96 (MQ=255)
                         gTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                                                  <  2:8269472/96‑1 (MQ=255)
                         gTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGAAATTTTGGTGATTTTCATTGCACTGCATCTCGATACGCCTTAACAAAtg                                                                  >  2:9382702/1‑96 (MQ=255)
                           ttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgct                                                                <  1:10369555/96‑1 (MQ=255)
                           atCAATCTTCAGGACCATCAAGCCGGGAGGTAACCGATACGTTGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtttt                                                                <  1:3524667/95‑4 (MQ=255)
                             cAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgc                                                              >  1:4763132/1‑96 (MQ=255)
                              aaTTTTCAGGAACTTCCAGCGATGATGTAAACGATACGTTGTAATTTTGGTGATTTTCATTGCACTGTCTCTCGATACTCCTTAACAAAtgctgct                                                             <  2:13905784/96‑1 (MQ=255)
                              aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgct                                                             <  1:1956541/96‑1 (MQ=255)
                              aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgct                                                             >  1:2212263/1‑96 (MQ=255)
                               aTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctg                                                            >  2:6109887/1‑96 (MQ=255)
                               aTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctg                                                            >  1:9984541/1‑96 (MQ=255)
                               aTTTTCAGGAACATCCAGCGGGGAGGTAAACGATAAGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctg                                                            <  1:3238239/96‑1 (MQ=255)
                                ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                                           <  2:2456986/96‑1 (MQ=255)
                                ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                                           >  1:5317374/1‑96 (MQ=255)
                                 tttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcc                                                          <  1:8725548/96‑1 (MQ=255)
                                 tttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcc                                                          >  1:13979096/1‑96 (MQ=255)
                                  ttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTCATTTTCATTGCACTGCCTCTCGATACGCCTTAACTAATGctgctgcct                                                         <  1:11004103/96‑1 (MQ=255)
                                   tCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                                        <  1:9698316/96‑1 (MQ=255)
                                   tCAGGAACATCCAGCGGGGAGGCAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGAt                                                                                 <  2:10285941/71‑1 (MQ=255)
                                     aGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                                      <  1:8907855/96‑1 (MQ=255)
                                     aGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                                      >  1:10016127/1‑96 (MQ=255)
                                      ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctg                                                     <  1:4099501/96‑1 (MQ=255)
                                        aaCATCCAGCTGTGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcg                                                   <  1:13385592/96‑1 (MQ=255)
                                        aaCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcg                                                   >  2:3580890/1‑96 (MQ=255)
                                        aaCATCCAGCGGGGAGGCAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcg                                                   >  1:6176949/1‑96 (MQ=255)
                                         aCATCCAGCGGGGAGGAAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcgc                                                  >  2:13233635/1‑96 (MQ=255)
                                            tCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGt                                               >  2:13741645/1‑96 (MQ=255)
                                            tCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGt                                               <  1:12567611/96‑1 (MQ=255)
                                              cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgc                                             <  1:8386083/96‑1 (MQ=255)
                                              cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgc                                             <  1:4982254/96‑1 (MQ=255)
                                              cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgc                                             <  2:3163049/96‑1 (MQ=255)
                                                gtggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcgc                                           >  1:10587055/3‑96 (MQ=255)
                                                    ggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTa                                       <  1:708547/96‑1 (MQ=255)
                                                     ggggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                                      >  1:13016852/3‑96 (MQ=255)
                                                     gAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                                      >  2:1651327/1‑96 (MQ=255)
                                                      gggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                                     >  1:11392669/2‑96 (MQ=255)
                                                      aGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                                     >  2:5075138/1‑96 (MQ=255)
                                                       ggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACgg                                    <  2:13191326/96‑1 (MQ=255)
                                                         tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAt                                  >  2:12061665/1‑96 (MQ=255)
                                                          aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAtt                                 >  2:12559101/1‑96 (MQ=255)
                                                          aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAtt                                 >  2:8350317/1‑96 (MQ=255)
                                                          aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAtt                                 <  1:8633547/96‑1 (MQ=255)
                                                            aCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGc                               <  2:11459710/96‑1 (MQ=255)
                                                              gATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTCAGCTCTACGGATTGccc                             >  1:12773851/1‑96 (MQ=255)
                                                               aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCg                            >  1:13665633/1‑96 (MQ=255)
                                                               aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCg                            <  1:6629059/96‑1 (MQ=255)
                                                                tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg                           >  2:11019491/1‑96 (MQ=255)
                                                                tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg                           >  1:2045823/1‑96 (MQ=255)
                                                                 aCGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgcc                          <  1:7543569/96‑3 (MQ=255)
                                                                  cGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcg                         <  2:1050690/96‑1 (MQ=255)
                                                                  cGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcg                         >  1:8196377/1‑96 (MQ=255)
                                                                   gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc                        <  2:6195702/96‑1 (MQ=255)
                                                                   gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCAGCGCGGTGCGCTCTACGGATTGCCCGgcgc                        >  2:14125301/1‑96 (MQ=255)
                                                                       tAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATa                    >  2:4375893/1‑96 (MQ=255)
                                                                       tAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATa                    <  1:10770505/96‑1 (MQ=255)
                                                                         aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga                  >  1:256573/1‑96 (MQ=255)
                                                                         aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga                  >  2:14302655/1‑96 (MQ=255)
                                                                         aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga                  >  2:6477024/1‑96 (MQ=255)
                                                                          ttttGGTGATTTTCATTGCACTTCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATagag                 <  2:11274507/96‑1 (MQ=255)
                                                                          ttttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCAGTGCGCTCTACGGATTGCCCGGCGCGATagag                 <  2:2210028/96‑1 (MQ=255)
                                                                           tttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAgaga                >  1:14061122/1‑96 (MQ=255)
                                                                            ttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGAt               >  2:2674376/1‑96 (MQ=255)
                                                                             tGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGGGCGATAGAGATc              >  2:2930630/1‑96 (MQ=255)
                                                                             tGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATc              <  1:1012597/96‑1 (MQ=255)
                                                                               gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGc            >  1:3875558/1‑96 (MQ=255)
                                                                                tGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGTTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCt           >  2:12251626/1‑96 (MQ=255)
                                                                                tGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCt           <  2:12694559/96‑1 (MQ=255)
                                                                                tGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCt           <  2:9222295/96‑1 (MQ=255)
                                                                                 gATTTTCATTGCACTGCCTCTCGATACGTCATAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTg          <  1:10863522/96‑1 (MQ=255)
                                                                                 gATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTg          <  2:11906834/96‑1 (MQ=255)
                                                                                 gATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTg          >  2:9710984/1‑96 (MQ=255)
                                                                                  aTTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGTCGCGATATAGATCGCTGc         >  2:5320371/1‑96 (MQ=255)
                                                                                    tttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCt       <  2:700340/96‑1 (MQ=255)
                                                                                     ttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTa      <  2:2044115/96‑1 (MQ=255)
                                                                                       cATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGCTAGAGATCGCTGCCTAAg    >  1:1211535/1‑96 (MQ=255)
                                                                                       cATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAg    >  2:4105768/1‑96 (MQ=255)
                                                                                              |                                                                                           
GCCCCAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCC  >  NC_000913/3870536‑3870721

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