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

Only 100 of 263 total aligned reads displayed.

CCCAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA  >  NC_000913/3870538‑3870718
                                                                                            |                                                                                          
cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                            <  1:5960984/93‑1 (MQ=255)
cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                            >  1:5141073/1‑93 (MQ=255)
cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                            >  2:4619975/1‑93 (MQ=255)
 ccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAct                                                                                           >  2:3671370/1‑92 (MQ=255)
 ccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAct                                                                                           >  1:2047108/1‑92 (MQ=255)
      cGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCctc                                                                                     <  1:6175020/93‑1 (MQ=255)
      cGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCctc                                                                                     <  2:2616221/93‑1 (MQ=255)
        aCCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCtctc                                                                                   >  2:2052598/1‑93 (MQ=255)
         ccACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCg                                                                                  >  1:3903130/1‑93 (MQ=255)
          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                                 <  1:294290/93‑1 (MQ=255)
          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                                 <  1:2772132/93‑1 (MQ=255)
          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                                 <  2:4414235/93‑1 (MQ=255)
           aCGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGAt                                                                                <  1:1898606/93‑1 (MQ=255)
           aCGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGAt                                                                                >  2:1389725/1‑93 (MQ=255)
           aCGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGAt                                                                                <  2:281267/93‑1 (MQ=255)
             gCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATAc                                                                              <  1:2810638/93‑1 (MQ=255)
             gCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATAc                                                                              >  1:4662738/1‑93 (MQ=255)
              ccTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTTCACTGCCTCTCGATACg                                                                             >  2:3744993/1‑93 (MQ=255)
                 tCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCt                                                                          >  1:3236767/1‑93 (MQ=255)
                   aTCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                                        >  1:6269189/1‑93 (MQ=255)
                   aTCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                                        >  2:4278392/1‑93 (MQ=255)
                   aTCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                                        >  1:2415288/1‑93 (MQ=255)
                   aTCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                                        >  2:5767620/1‑93 (MQ=255)
                    tCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTaa                                                                       <  2:3485906/93‑1 (MQ=255)
                     cgctttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAAc                                                                      <  1:800144/90‑1 (MQ=255)
                       gTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaa                                                                    >  1:1430207/1‑93 (MQ=255)
                       gTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaa                                                                    >  1:4185102/1‑93 (MQ=255)
                        tttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaaa                                                                   <  2:5732980/93‑1 (MQ=255)
                        tttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaaa                                                                   >  2:2214303/1‑93 (MQ=255)
                        tttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACaaa                                                                   >  1:4252864/1‑93 (MQ=255)
                         ttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                                                  <  2:2911762/93‑1 (MQ=255)
                         ttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                                                  <  1:450681/93‑1 (MQ=255)
                          tCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                                                 >  1:2560945/1‑93 (MQ=255)
                          tCAATTTTCAGGAACATCCAGCGGGGAGGTAAACAATACGCGGGAATTTTGGTGATCTTCATTGCACTGCCTCTCGATACGCCCTAACaaaaa                                                                 >  2:4049480/1‑91 (MQ=255)
                            aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgct                                                               <  1:953977/93‑1 (MQ=255)
                            aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgct                                                               >  1:271878/1‑93 (MQ=255)
                              ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgc                                                            >  2:3848396/1‑93 (MQ=255)
                               tttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgct                                                           >  2:6186390/1‑93 (MQ=255)
                                 tCAGGAACATCCAGCGGGGCGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                                         >  2:6129820/1‑93 (MQ=255)
                                 tCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                                         >  1:3840513/1‑93 (MQ=255)
                                 tCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                                         >  1:653146/1‑93 (MQ=255)
                                 tCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                                         >  1:1238404/1‑93 (MQ=255)
                                 tCAGGAACATCCAGCGGGGAGGAAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGACTCTCGATACGCCTTAACAAAtgctgctgc                                                         >  2:982112/1‑93 (MQ=255)
                                  cAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcc                                                        >  1:2232796/1‑93 (MQ=255)
                                   aGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgcgt                                                       >  1:5298602/1‑91 (MQ=255)
                                    ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                                      >  1:1282251/1‑93 (MQ=255)
                                    ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                                      >  2:197953/1‑93 (MQ=255)
                                    ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                                      >  2:1299589/1‑93 (MQ=255)
                                    ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                                      >  1:2102995/1‑93 (MQ=255)
                                     gAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgc                                                     <  2:5615238/93‑1 (MQ=255)
                                      aaCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                                    <  1:3148117/93‑1 (MQ=255)
                                        cATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctgc                                                  <  2:643202/93‑1 (MQ=255)
                                        cATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctgc                                                  <  1:3779478/93‑1 (MQ=255)
                                        cATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctgc                                                  <  1:6069796/93‑1 (MQ=255)
                                          tCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcgc                                                >  2:5088405/1‑93 (MQ=255)
                                            cAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCgg                                              >  2:2558726/1‑93 (MQ=255)
                                              gCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTg                                            <  2:1717937/93‑1 (MQ=255)
                                                ggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcg                                          >  2:141959/1‑93 (MQ=255)
                                                   gAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGctc                                      >  1:4367017/1‑93 (MQ=255)
                                                    aGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGctct                                     >  1:3980484/1‑93 (MQ=255)
                                                     ggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTa                                    >  2:1572556/1‑93 (MQ=255)
                                                     ggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTa                                    >  2:5881747/1‑93 (MQ=255)
                                                     ggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTa                                    >  2:239901/1‑93 (MQ=255)
                                                       tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                                  >  1:5891626/1‑93 (MQ=255)
                                                       tAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACg                                  >  2:1494333/1‑93 (MQ=255)
                                                          aCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAt                               <  1:5412952/93‑1 (MQ=255)
                                                          aCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAt                               <  1:1577625/93‑1 (MQ=255)
                                                           cGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAtt                              <  1:1394343/93‑1 (MQ=255)
                                                           cGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAtt                              <  2:5184326/93‑1 (MQ=255)
                                                           cGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAtt                              >  2:780937/1‑93 (MQ=255)
                                                               aCGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGccc                          <  1:1805329/93‑1 (MQ=255)
                                                               aCGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGccc                          >  1:3665/1‑93 (MQ=255)
                                                               aCGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGccc                          <  2:848315/93‑1 (MQ=255)
                                                                 gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg                        >  2:5461948/1‑93 (MQ=255)
                                                                 gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAAAAAATCTGCCTGCCGCCGCGCCGTTCGGTCCTACGGATTGCCCgc                        >  2:4879400/1‑92 (MQ=255)
                                                                  tGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgc                       >  1:2336592/1‑93 (MQ=255)
                                                                  tGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgc                       <  2:3369442/93‑1 (MQ=255)
                                                                    gTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc                     <  2:5946380/93‑1 (MQ=255)
                                                                     tAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGcgcg                    <  1:1735789/93‑1 (MQ=255)
                                                                       aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGAt                  >  2:5353536/1‑93 (MQ=255)
                                                                        ttttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATa                 >  2:398339/1‑93 (MQ=255)
                                                                         tttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATag                <  1:3630583/93‑1 (MQ=255)
                                                                          ttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga               >  1:963870/1‑93 (MQ=255)
                                                                          ttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATaga               <  1:5591195/93‑1 (MQ=255)
                                                                           tGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATagag              >  1:5762443/1‑93 (MQ=255)
                                                                            ggTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAgaga             <  1:2166174/93‑1 (MQ=255)
                                                                             gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGAt            <  1:1122925/93‑1 (MQ=255)
                                                                               gATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCg          >  2:1660087/1‑93 (MQ=255)
                                                                                agtgtCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGc         <  2:4698631/89‑1 (MQ=255)
                                                                                aTTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGc         <  1:1000479/93‑1 (MQ=255)
                                                                                aTTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGc         <  1:1489671/93‑1 (MQ=255)
                                                                                aTTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGc         >  2:4540639/1‑93 (MQ=255)
                                                                                 ttttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCt        >  1:3559860/1‑93 (MQ=255)
                                                                                  tttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTg       >  2:18719/1‑93 (MQ=255)
                                                                                   ttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGc      <  1:3467558/93‑1 (MQ=255)
                                                                                   ttCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGc      <  2:2757827/93‑1 (MQ=255)
                                                                                    tCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGcc     >  2:600130/1‑93 (MQ=255)
                                                                                      aTTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTa   >  1:4929924/1‑93 (MQ=255)
                                                                                      aTTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTa   >  2:493896/1‑93 (MQ=255)
                                                                                       ttGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTaa  >  1:4839756/1‑93 (MQ=255)
                                                                                            |                                                                                          
CCCAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA  >  NC_000913/3870538‑3870718

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