Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,067,548 T→C H137R (CAC→CGC)  b3879 ← putative aldose‑1‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,067,5480TC100.0% 586.9 / NA 186H137R (CAC→CGC) b3879putative aldose‑1‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (92/94);  total (92/94)

Only 100 of 220 total aligned reads displayed.

TTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTG  >  NC_000913/4067458‑4067637
                                                                                          |                                                                                         
ttGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCgtgcgt                                                                                         <  1:4524978/93‑1 (MQ=255)
ttGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCgtgcgt                                                                                         <  1:1477486/93‑1 (MQ=255)
  gCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCgtgcgtgc                                                                                       <  1:29320/93‑1 (MQ=255)
  gCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCgtgcgtgc                                                                                       >  2:2610461/1‑93 (MQ=255)
  gCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCgtgcgtgc                                                                                       <  1:4626443/93‑1 (MQ=255)
    aGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTgcgtgcgt                                                                                     >  1:2411089/1‑93 (MQ=255)
     ggTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGGCCCGTGCGTGCGtt                                                                                    >  2:504077/1‑93 (MQ=255)
     ggTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGtt                                                                                    >  1:4559839/1‑93 (MQ=255)
     ggTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGtt                                                                                    <  1:4658488/93‑1 (MQ=255)
     ggTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGtt                                                                                    >  1:1332110/1‑93 (MQ=255)
     ggTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGtt                                                                                    >  2:2471516/1‑93 (MQ=255)
      gTTAAATAAGGATGAATACCCGCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTg                                                                                   <  2:4526427/93‑1 (MQ=255)
      gTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTg                                                                                   >  2:698435/1‑93 (MQ=255)
      gTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTg                                                                                   <  1:2050656/93‑1 (MQ=255)
         aaaTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGc                                                                                >  2:1832462/1‑93 (MQ=255)
         aaaTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGGGCGTTGGGc                                                                                >  1:365891/1‑93 (MQ=255)
          aaTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCg                                                                               >  1:1941734/1‑93 (MQ=255)
           aTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGa                                                                              >  2:4007487/1‑93 (MQ=255)
           aTAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGa                                                                              >  1:2174757/1‑93 (MQ=255)
            tAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGaa                                                                             >  1:2583802/1‑93 (MQ=255)
                gATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaaa                                                                         <  1:3573346/93‑1 (MQ=255)
                gATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATaaa                                                                         <  1:2133182/93‑1 (MQ=255)
                  tGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaaca                                                                       >  2:3952571/1‑93 (MQ=255)
                  tGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAaaca                                                                       >  1:1601972/1‑93 (MQ=255)
                      tACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACtt                                                                   >  2:2563808/1‑93 (MQ=255)
                      tACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACtt                                                                   >  2:3574491/1‑93 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 <  1:3157504/93‑1 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 <  2:353533/93‑1 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 <  1:3411589/93‑1 (MQ=255)
                        cccACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTtgt                                                                 >  1:1020592/1‑93 (MQ=255)
                          cACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGa                                                               <  2:2715661/93‑1 (MQ=255)
                           ccGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              <  1:497481/92‑1 (MQ=255)
                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              <  1:3513195/93‑1 (MQ=255)
                           aCGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGaa                                                              >  1:972772/1‑93 (MQ=255)
                            cGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAg                                                             <  1:2916569/93‑1 (MQ=255)
                             gCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGc                                                            >  2:3906631/1‑93 (MQ=255)
                             gCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGc                                                            >  1:4302827/1‑93 (MQ=255)
                             gCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGc                                                            >  2:3489427/1‑93 (MQ=255)
                              ccATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGCCCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGcc                                                           <  2:2682013/93‑1 (MQ=255)
                              ccATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGcc                                                           <  2:4254755/93‑1 (MQ=255)
                              ccATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGcc                                                           >  1:3252035/1‑93 (MQ=255)
                              ccATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACCTGTGAAGcc                                                           >  2:4379708/1‑93 (MQ=255)
                               cATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCa                                                          <  2:2025800/93‑1 (MQ=255)
                                 tAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGc                                                        >  2:479533/1‑93 (MQ=255)
                                 tAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGc                                                        >  2:4100119/1‑93 (MQ=255)
                                   aGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCAt                                                      <  2:3356834/93‑1 (MQ=255)
                                    ggCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATa                                                     >  1:3648965/1‑93 (MQ=255)
                                     gcgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaa                                                    <  1:63258/93‑1 (MQ=255)
                                      cgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaa                                                   <  1:751249/93‑1 (MQ=255)
                                      cgcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATaaa                                                   <  2:3026374/93‑1 (MQ=255)
                                       gcTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAAc                                                  <  1:3622697/93‑1 (MQ=255)
                                        cTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACg                                                 <  2:3165359/93‑1 (MQ=255)
                                         tGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACgg                                                <  1:1178971/93‑1 (MQ=255)
                                            gACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATa                                             <  1:1790940/93‑1 (MQ=255)
                                             aCAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAg                                            >  1:208848/1‑93 (MQ=255)
                                             aCAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAg                                            >  2:1732756/1‑93 (MQ=255)
                                              cAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGc                                           <  1:2530557/93‑1 (MQ=255)
                                               aGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGcc                                          <  1:3401223/93‑1 (MQ=255)
                                               aGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGcc                                          <  1:4005527/93‑1 (MQ=255)
                                                    cGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAAc                                     <  1:78495/93‑1 (MQ=255)
                                                       tatTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTgg                                  <  1:1679817/93‑1 (MQ=255)
                                                       tatTCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTgg                                  >  2:1096431/1‑93 (MQ=255)
                                                         ttCTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGCAGCCAGCATAAACGGATAGCCGTAACTgggg                                >  1:3783942/1‑92 (MQ=255)
                                                           cTGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggcgg                              <  2:4542822/93‑1 (MQ=255)
                                                            tGGCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGggcggc                             <  2:3827674/93‑1 (MQ=255)
                                                             ggCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCa                            <  1:4137374/93‑1 (MQ=255)
                                                              gCTTGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTTAAGCCAGCATAAACGGATAACCGTAACTTGTCGTCAg                           >  2:4484775/1‑93 (MQ=255)
                                                                 tGCGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGaaa                        <  1:2707540/93‑1 (MQ=255)
                                                                   cGATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAg                      >  2:2792764/1‑93 (MQ=255)
                                                                    gATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGc                     <  2:2103465/93‑1 (MQ=255)
                                                                    gATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGc                     <  1:3017592/93‑1 (MQ=255)
                                                                    gATCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGc                     <  1:1759713/93‑1 (MQ=255)
                                                                     aTCTCAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAACGCt                    >  2:3374367/1‑93 (MQ=255)
                                                                      tctcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGTGCGGCAGAAAAGCTg                   >  2:2839226/1‑93 (MQ=255)
                                                                        tcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTc                 <  2:4135638/93‑1 (MQ=255)
                                                                        tcAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTc                 <  1:401902/93‑1 (MQ=255)
                                                                         cAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCa                >  2:85691/1‑93 (MQ=255)
                                                                         cAACGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCa                <  1:1988319/93‑1 (MQ=255)
                                                                          aaCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAg               <  2:3704342/93‑1 (MQ=255)
                                                                          aaCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAg               <  1:2223334/93‑1 (MQ=255)
                                                                          aaCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAg               <  2:535831/93‑1 (MQ=255)
                                                                           aCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              <  1:1259415/93‑1 (MQ=255)
                                                                           aCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              >  1:3948986/1‑93 (MQ=255)
                                                                           aCGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAgg              >  2:4193539/1‑93 (MQ=255)
                                                                            cGGAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAggg             <  1:1867594/93‑1 (MQ=255)
                                                                              gAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTa           <  2:4663075/93‑1 (MQ=255)
                                                                              gAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTa           <  1:984307/93‑1 (MQ=255)
                                                                              gAAAGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTa           >  2:1402176/1‑93 (MQ=255)
                                                                                aaGACCCGTTCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAAc         <  2:1526214/93‑1 (MQ=255)
                                                                                aaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAAc         <  2:2400137/93‑1 (MQ=255)
                                                                                aaGACCCGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGAAAc         >  1:4228265/1‑93 (MQ=255)
                                                                                  gccccGGGGGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGc       <  2:1230386/91‑1 (MQ=255)
                                                                                  gACCCGTGCGTGCGTTGAGCAAATAAGCAACTGGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGc       <  2:738534/93‑1 (MQ=255)
                                                                                     ccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtgg    >  2:4343940/1‑93 (MQ=255)
                                                                                     ccGTGCGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtgg    >  1:4070528/1‑93 (MQ=255)
                                                                                      cgtgcgtgCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggt   <  1:2769863/93‑1 (MQ=255)
                                                                                      cgtgcgtgCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggt   <  1:908271/93‑1 (MQ=255)
                                                                                      cgtgcgtgCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggt   >  2:1710333/1‑93 (MQ=255)
                                                                                      cgtgcgtgCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggt   >  2:1974732/1‑93 (MQ=255)
                                                                                       gtgcgtgcGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCtggtg  <  2:1653277/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTGCAGGTTAAATAAGGATGAATACCCACGCCATAAGGCGCTGCGACAGTGCCGATATTCTGGCTTGCGATCTCAACGGAAAGACCCGTGTGTGCGTTGAGCGAATAAACAACTTGTGAAGCCAGCATAAACGGATAGCCGTAACTGGGCGGCAGAAAAGCTGTCAGGGTAACGCTGGTG  >  NC_000913/4067458‑4067637

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