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

Only 100 of 297 total aligned reads displayed.

CCCAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACCTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAA  >  NC_000913/3870538‑3870718
                                                                                            |                                                                                        
cccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAc                                                                                          >  2:1103860/1‑93 (MQ=255)
 ccAACCGACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCAct                                                                                         <  2:1637216/93‑2 (MQ=255)
     ccGACCACGCCTGCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCct                                                                                    <  1:3565339/93‑1 (MQ=255)
       gtccATGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCctct                                                                                  <  1:4188403/91‑1 (MQ=255)
       gACCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTTATTGCACTGCctct                                                                                  >  2:2968114/1‑93 (MQ=255)
        aCCACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCtctc                                                                                 <  1:5468115/93‑1 (MQ=255)
         ccACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCg                                                                                >  1:3697438/1‑93 (MQ=255)
          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                               <  1:5543824/93‑1 (MQ=255)
          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                               <  1:4329756/93‑1 (MQ=255)
          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                               >  1:192171/1‑93 (MQ=255)
          cACGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGa                                                                               <  1:2827375/93‑1 (MQ=255)
           aCGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGAt                                                                              >  2:1693218/1‑93 (MQ=255)
            cGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATa                                                                             <  1:4960177/93‑1 (MQ=255)
            cGCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATa                                                                             <  2:1289952/93‑1 (MQ=255)
             gCCTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATAc                                                                            >  1:3795353/1‑93 (MQ=255)
              ccTTCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACg                                                                           >  1:2997261/1‑93 (MQ=255)
                ttCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGcc                                                                         <  1:4807093/93‑1 (MQ=255)
                 tCATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCt                                                                        <  2:3906342/93‑1 (MQ=255)
                  cATCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCtt                                                                       >  1:262411/1‑93 (MQ=255)
                   aTCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTa                                                                      >  2:384186/1‑93 (MQ=255)
                    tCGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTaa                                                                     >  1:5597036/1‑93 (MQ=255)
                     cGGTTTCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAAc                                                                    >  1:771880/1‑93 (MQ=255)
                         ttCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAt                                                                <  1:2566863/93‑1 (MQ=255)
                          tCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                                               <  2:5348993/93‑1 (MQ=255)
                          tCAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtg                                                               <  2:2616986/93‑1 (MQ=255)
                           cAATTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgc                                                              <  1:3933722/93‑1 (MQ=255)
                            aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgct                                                             >  2:5790051/1‑93 (MQ=255)
                            aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgct                                                             <  2:3212190/93‑1 (MQ=255)
                            aaTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGGAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTCACAAAtgct                                                             >  1:346786/1‑93 (MQ=255)
                             aTTTTCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctg                                                            <  1:2092249/93‑1 (MQ=255)
                              ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgc                                                           <  2:5740956/93‑1 (MQ=255)
                              ttttCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgc                                                           >  2:2477163/1‑93 (MQ=255)
                                 tCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                                        <  2:4110417/93‑1 (MQ=255)
                                 tCAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAAtgctgctgc                                                        <  2:6021965/93‑1 (MQ=255)
                                  cAGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcc                                                       <  1:4896037/93‑1 (MQ=255)
                                   aGGACCATCCAGCGTGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcct                                                      <  1:6017101/93‑1 (MQ=255)
                                   aGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcct                                                      <  2:5043007/93‑1 (MQ=255)
                                   aGGAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcct                                                      >  2:944513/1‑93 (MQ=255)
                                    ggAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctg                                                     >  1:1555693/1‑93 (MQ=255)
                                     gAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgc                                                    <  2:1255692/93‑1 (MQ=255)
                                     gAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgc                                                    <  1:3332028/93‑1 (MQ=255)
                                     gAACATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgc                                                    <  1:6413209/93‑1 (MQ=255)
                                      aaCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                                   >  2:2262953/1‑93 (MQ=255)
                                      aaCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                                   >  1:2747366/1‑93 (MQ=255)
                                      aaCATCCAGCGGGGAGGGAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgct                                                   >  1:1749525/1‑93 (MQ=255)
                                       aCATCTAGCGGGAAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctg                                                  <  1:944232/93‑1 (MQ=255)
                                       aCATCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGctgctgcctgctg                                                  >  2:5690331/1‑93 (MQ=255)
                                         aTCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcg                                                <  1:4474772/93‑1 (MQ=255)
                                          tCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcgc                                               <  1:2419229/93‑1 (MQ=255)
                                          tCCAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTgcgc                                               <  1:145457/93‑1 (MQ=255)
                                           ccAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGcgcg                                              >  2:2631407/1‑93 (MQ=255)
                                           ccAGCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGcgcg                                              >  2:559597/1‑93 (MQ=255)
                                              gCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTg                                           >  2:2512001/1‑93 (MQ=255)
                                              gCGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTg                                           <  2:6345585/93‑1 (MQ=255)
                                               cGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgc                                          <  2:4321087/93‑1 (MQ=255)
                                               cGGGGAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgc                                          <  2:64541/93‑1 (MQ=255)
                                                ggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcg                                         >  1:2946199/1‑93 (MQ=255)
                                                ggggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcg                                         >  2:4439491/1‑93 (MQ=255)
                                                 gggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcgc                                        >  1:1234193/1‑93 (MQ=255)
                                                 gggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTgcgc                                        >  2:209731/1‑93 (MQ=255)
                                                  ggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGct                                       <  2:4532416/93‑1 (MQ=255)
                                                  ggAGGTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCCGCTGCGAGGTGCGct                                       >  1:6148782/1‑93 (MQ=255)
                                                     ggTAAACGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTa                                    >  2:5144717/1‑93 (MQ=255)
                                                        aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACgg                                 >  1:2192210/1‑93 (MQ=255)
                                                        aaaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACgg                                 >  1:5104669/1‑93 (MQ=255)
                                                         aaCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGa                                >  1:1068928/1‑93 (MQ=255)
                                                          aCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAt                               <  2:5436169/93‑1 (MQ=255)
                                                          aCGATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGAt                               <  2:4790281/93‑1 (MQ=255)
                                                            gATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTg                             >  2:2669907/1‑93 (MQ=255)
                                                            gATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTg                             >  1:409727/1‑93 (MQ=255)
                                                            gATACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTg                             >  1:319679/1‑93 (MQ=255)
                                                             aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGc                            >  1:5039242/1‑93 (MQ=255)
                                                             aTACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGc                            >  1:4641309/1‑93 (MQ=255)
                                                              tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGcc                           <  1:2668658/93‑1 (MQ=255)
                                                              tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGcc                           >  1:99871/1‑93 (MQ=255)
                                                              tACGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGcc                           <  2:1539362/93‑1 (MQ=255)
                                                               aCGTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGccc                          <  1:6338388/93‑1 (MQ=255)
                                                                 gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg                        >  2:339214/1‑93 (MQ=255)
                                                                 gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg                        <  2:3541799/93‑1 (MQ=255)
                                                                 gTGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCgg                        >  2:362761/1‑93 (MQ=255)
                                                                  tgttAATTTGGGTGATTTTCATTGCACTGCCGCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgc                       <  1:3107443/90‑1 (MQ=255)
                                                                  tGGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgc                       >  2:5909141/1‑93 (MQ=255)
                                                                    gTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc                     >  2:4906946/1‑93 (MQ=255)
                                                                    gTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc                     >  1:5823360/1‑93 (MQ=255)
                                                                    gTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGgcgc                     <  2:2784243/93‑1 (MQ=255)
                                                                     tAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGcgcg                    >  1:3037604/1‑93 (MQ=255)
                                                                       aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGAt                  <  2:181162/93‑1 (MQ=255)
                                                                       aTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGAt                  <  2:3683216/93‑1 (MQ=255)
                                                                         tttGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATag                <  2:1970173/93‑1 (MQ=255)
                                                                            ggTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAgaga             >  1:4850553/1‑93 (MQ=255)
                                                                             gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGAt            >  2:1455638/1‑93 (MQ=255)
                                                                             gTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGAt            <  1:5953037/93‑1 (MQ=255)
                                                                              tGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATc           <  1:2369613/93‑1 (MQ=255)
                                                                               gATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCg          >  1:3243192/1‑93 (MQ=255)
                                                                                aTTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGc         <  1:215327/93‑1 (MQ=255)
                                                                                  tttCATTGCACTGCCTCTCGATACGCTTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTg       <  1:5201855/93‑1 (MQ=255)
                                                                                    tCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGcc     <  1:5354375/93‑1 (MQ=255)
                                                                                    tCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGcc     >  2:692053/1‑93 (MQ=255)
                                                                                       ttGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTaa  >  2:179672/1‑93 (MQ=255)
                                                                                       ttGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTaa  <  2:311762/93‑1 (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 < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: