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.100.0% 1134.8 / NA 247coding (611/1764 nt)dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (129/118);  total (129/118)

Only 100 of 267 total aligned reads displayed.

CCCAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA  >  NC_000913/3870538‑3870718
                                                                                            |                                                                                        
cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                          <  1:668356/93‑1 (MQ=255)
cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                          <  1:3454442/93‑1 (MQ=255)
cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                          <  1:3743442/93‑1 (MQ=255)
cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                          <  2:2020054/93‑1 (MQ=255)
  cAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActg                                                                                        <  2:302329/93‑3 (MQ=255)
  cAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActg                                                                                        >  2:2326659/1‑91 (MQ=255)
   aaCCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActgc                                                                                       >  2:3599883/1‑90 (MQ=255)
    aCCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCActgcc                                                                                      <  2:1377126/93‑5 (MQ=255)
     ccGCCCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGCTACGTGGTAATTTTGGTGATTTTCATTGCACTGCct                                                                                    >  2:4641516/1‑93 (MQ=255)
      cGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCctc                                                                                   >  2:3052790/1‑93 (MQ=255)
      cGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCctc                                                                                   <  1:4531548/93‑1 (MQ=255)
         ccACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCg                                                                                <  2:363556/93‑1 (MQ=255)
          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                               >  2:1604168/1‑93 (MQ=255)
          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                               <  2:3987931/93‑1 (MQ=255)
          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                               <  1:1394200/93‑1 (MQ=255)
            cGCCTTCATCGGTTTCAATTTTCAGGAGCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATa                                                                             >  2:3152254/1‑93 (MQ=255)
            cGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATa                                                                             <  1:4286330/93‑1 (MQ=255)
            cGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATa                                                                             <  2:4388373/93‑1 (MQ=255)
             gCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATAc                                                                            <  1:2995765/93‑1 (MQ=255)
              ccTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACg                                                                           >  1:3636896/1‑93 (MQ=255)
               cTTCATCGGTTTCAATTTTCAGGAACATGCATCGGGGAGGTTAATGATACGTGGTAACTTTGGTGATTTGTAATCCACTGCCTCTCGATACGc                                                                          >  2:89624/1‑93 (MQ=255)
               cTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGc                                                                          >  2:210420/1‑93 (MQ=255)
                 tCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCt                                                                        <  2:3798314/93‑1 (MQ=255)
                   aTCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                                      >  2:1361196/1‑93 (MQ=255)
                    tCGGTTTTAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTTATTTTCATTGCACTGCCTCTCGATACGCCTTaa                                                                     <  1:2767759/93‑1 (MQ=255)
                    tCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTaa                                                                     >  2:2358244/1‑93 (MQ=255)
                    tCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTaa                                                                     <  1:331009/93‑1 (MQ=255)
                     cGGTTTCAATTTTCAGGACCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAAc                                                                    <  1:842313/93‑1 (MQ=255)
                      ggTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACa                                                                   >  2:2457553/1‑93 (MQ=255)
                       gTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaa                                                                  >  1:1559739/1‑93 (MQ=255)
                        tttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaaa                                                                 >  1:2299077/1‑93 (MQ=255)
                        tttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaaa                                                                 <  2:2949983/93‑1 (MQ=255)
                         ttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                                                <  2:3867879/93‑1 (MQ=255)
                          tCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                                               >  2:1144332/1‑93 (MQ=255)
                          tCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                                               <  1:1986254/93‑1 (MQ=255)
                          tCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                                               <  2:2749345/93‑1 (MQ=255)
                           cAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgc                                                              >  1:1642522/1‑93 (MQ=255)
                            aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgct                                                             >  2:72415/1‑93 (MQ=255)
                             aTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctg                                                            >  1:3880213/1‑93 (MQ=255)
                              ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgc                                                           <  2:2429926/93‑1 (MQ=255)
                                 tCAGGAACATCCAGCGGGGAGGTAAACGATACGTTGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                                        <  2:4203524/93‑1 (MQ=255)
                                   aGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcct                                                      <  1:403360/93‑1 (MQ=255)
                                   aGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcct                                                      <  1:2538483/93‑1 (MQ=255)
                                    ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                                     >  1:733235/1‑93 (MQ=255)
                                      aaCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                                   <  2:2883133/93‑1 (MQ=255)
                                      aaCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                                   >  2:1411791/1‑93 (MQ=255)
                                       aCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctg                                                  >  1:1008981/1‑93 (MQ=255)
                                          tCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcgc                                               >  2:4326806/1‑93 (MQ=255)
                                           ccAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGcgcg                                              <  1:3277502/93‑1 (MQ=255)
                                            cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCgg                                             >  1:996053/1‑93 (MQ=255)
                                               cGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgc                                          <  1:1301086/93‑1 (MQ=255)
                                               cGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgc                                          <  2:628672/93‑1 (MQ=255)
                                                 gggAGGTAAGCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcgc                                        <  1:1875212/93‑1 (MQ=255)
                                                      gTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                                   <  1:2267354/93‑1 (MQ=255)
                                                      gTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                                   <  2:451055/93‑1 (MQ=255)
                                                      gTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTAc                                   <  1:2517028/93‑1 (MQ=255)
                                                         aaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGa                                <  2:2392856/93‑1 (MQ=255)
                                                         aaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGa                                <  2:4633522/93‑1 (MQ=255)
                                                          aCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAt                               <  1:427078/93‑1 (MQ=255)
                                                          aCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAt                               >  1:4195659/1‑93 (MQ=255)
                                                          aCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAt                               <  2:783852/93‑1 (MQ=255)
                                                           cGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAtt                              <  1:941008/93‑1 (MQ=255)
                                                           cGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAtt                              <  1:2951707/93‑1 (MQ=255)
                                                             aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGc                            >  1:2891387/1‑93 (MQ=255)
                                                             aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGc                            >  2:4383516/1‑93 (MQ=255)
                                                             aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGc                            >  1:3962640/1‑93 (MQ=255)
                                                             aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGc                            >  1:4096367/1‑93 (MQ=255)
                                                              tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGcc                           >  2:3916455/1‑93 (MQ=255)
                                                              tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGcc                           <  2:4074948/93‑1 (MQ=255)
                                                              tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGcc                           >  1:1101265/1‑93 (MQ=255)
                                                              tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGcc                           <  1:4587492/93‑1 (MQ=255)
                                                                cGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCg                         >  1:1961632/1‑93 (MQ=255)
                                                                 gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg                        <  1:1749564/93‑1 (MQ=255)
                                                                 gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg                        <  1:3061274/93‑1 (MQ=255)
                                                                  tGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgc                       >  2:936930/1‑93 (MQ=255)
                                                                  tGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgc                       <  2:2511259/93‑1 (MQ=255)
                                                                  tGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgc                       >  1:620091/1‑93 (MQ=255)
                                                                  tGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgc                       >  1:3875645/1‑93 (MQ=255)
                                                                    ggaaTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc                     >  2:3733028/3‑93 (MQ=255)
                                                                        ttttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATa                 <  1:3621730/93‑1 (MQ=255)
                                                                          ttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga               <  2:849134/93‑1 (MQ=255)
                                                                          ttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga               <  1:2847538/93‑1 (MQ=255)
                                                                          ttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga               <  2:1921445/93‑1 (MQ=255)
                                                                          ttGGTGATTTTCATTACACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga               <  2:667364/93‑1 (MQ=255)
                                                                           tGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATagag              <  2:2858717/93‑1 (MQ=255)
                                                                            ggTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAgaga             >  2:3315932/1‑93 (MQ=255)
                                                                             gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGAt            <  1:2969367/93‑1 (MQ=255)
                                                                               gATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCg          >  2:2974724/1‑93 (MQ=255)
                                                                                aTTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGc         >  1:4336686/1‑93 (MQ=255)
                                                                                aTTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGc         <  1:1269500/93‑1 (MQ=255)
                                                                                 ttttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCt        <  2:3451565/93‑1 (MQ=255)
                                                                                   ttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGc      >  2:2644951/1‑93 (MQ=255)
                                                                                   ttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGc      >  1:1469746/1‑93 (MQ=255)
                                                                                     cATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCt    >  2:4023535/1‑93 (MQ=255)
                                                                                      aTTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTa   >  2:2238304/1‑93 (MQ=255)
                                                                                      aTTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTa   >  1:2436100/1‑93 (MQ=255)
                                                                                      aTTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTa   >  2:3294398/1‑93 (MQ=255)
                                                                                       ttGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTaa  <  2:4262734/93‑1 (MQ=255)
                                                                                       ttGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTaa  <  2:1802369/93‑1 (MQ=255)
                                                                                       ttGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTaa  <  2:65224/93‑1 (MQ=255)
                                                                                            |                                                                                        
CCCAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA  >  NC_000913/3870538‑3870718

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