Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,039,960 Δ1 bp intergenic (+42/‑36) yihD → / → yihE hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,039,9590A.100.0% 2709.2 / NA 592intergenic (+41/‑37)yihD/yihEhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (262/330);  total (262/330)

Only 100 of 639 total aligned reads displayed.

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047
                                                                                            |                                                                                        
agaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                                          <  1:7805848/93‑1 (MQ=255)
agaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                                          <  2:2099222/93‑1 (MQ=255)
 gaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                                         <  2:8179754/93‑2 (MQ=255)
   agaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcc                                                                                       >  1:9301370/1‑90 (MQ=255)
     aCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCtt                                                                                    <  1:1954813/93‑1 (MQ=255)
        tCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGt                                                                                 >  1:9675535/1‑93 (MQ=255)
         cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTc                                                                                <  1:1157829/93‑1 (MQ=255)
         cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTc                                                                                >  2:7894397/1‑93 (MQ=255)
         cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTc                                                                                <  1:537373/93‑1 (MQ=255)
         cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTc                                                                                >  2:5528692/1‑93 (MQ=255)
          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCa                                                                               <  2:7106762/93‑1 (MQ=255)
          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCa                                                                               <  2:7410009/93‑1 (MQ=255)
           aaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCAt                                                                              <  2:2661925/93‑1 (MQ=255)
           aaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCATTCGTCAt                                                                              >  1:341656/1‑93 (MQ=255)
            aaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCAtt                                                                             >  2:6857447/1‑93 (MQ=255)
             aCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGCGAAATGATATCCTTCGTCATTc                                                                            >  1:4918071/1‑93 (MQ=255)
              cGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCg                                                                           >  2:5680552/1‑93 (MQ=255)
                 cgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGCAATGATATCCTTCGTCATTCGTaa                                                                        >  2:9305575/1‑93 (MQ=255)
                 cgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTaa                                                                        <  1:6884165/93‑1 (MQ=255)
                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAAt                                                                       <  1:5785356/93‑1 (MQ=255)
                   cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATg                                                                      >  2:8381730/1‑93 (MQ=255)
                    tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGt                                                                     <  1:669671/93‑1 (MQ=255)
                    tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGt                                                                     >  2:8738049/1‑93 (MQ=255)
                     gTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtt                                                                    <  2:4531592/93‑1 (MQ=255)
                      ttACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGttt                                                                   >  1:3137112/1‑93 (MQ=255)
                        aCGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTc                                                                 <  1:1728924/93‑1 (MQ=255)
                         cgcgCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTcc                                                                <  1:8770641/93‑1 (MQ=255)
                          gcgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCg                                                               <  2:2769720/93‑1 (MQ=255)
                           cgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCGTCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCgg                                                              <  2:7534354/93‑1 (MQ=255)
                           cgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCgg                                                              >  2:7617155/1‑93 (MQ=255)
                            gcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGga                                                             <  1:3968317/93‑1 (MQ=255)
                             cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            >  1:7447502/1‑93 (MQ=255)
                              tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatg                                                           >  2:7671564/1‑93 (MQ=255)
                              tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatg                                                           <  2:118492/93‑1 (MQ=255)
                               cgcgGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatga                                                          <  1:7169166/93‑1 (MQ=255)
                               cgcgGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatga                                                          <  2:905474/93‑1 (MQ=255)
                               cgcgGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatga                                                          >  1:10000811/1‑93 (MQ=255)
                                gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatgat                                                         <  2:3350284/93‑1 (MQ=255)
                                gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatgat                                                         >  1:5656656/1‑93 (MQ=255)
                                gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatgat                                                         <  2:1133830/93‑1 (MQ=255)
                                 cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGatgatg                                                        <  2:8031194/93‑1 (MQ=255)
                                 cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGatgatg                                                        >  1:8532834/1‑93 (MQ=255)
                                  ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATgg                                                       <  1:2906839/93‑1 (MQ=255)
                                   gcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                                      >  2:9428183/1‑93 (MQ=255)
                                     gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                                    <  1:1351303/93‑1 (MQ=255)
                                      tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                                   >  2:6431688/1‑93 (MQ=255)
                                       aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                                  <  1:9061892/93‑1 (MQ=255)
                                       aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                                  <  2:6311101/93‑1 (MQ=255)
                                        aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                 >  1:358783/1‑93 (MQ=255)
                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                <  1:822018/93‑1 (MQ=255)
                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                >  2:5387359/1‑93 (MQ=255)
                                          tAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                               <  2:4388301/93‑1 (MQ=255)
                                           aaaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                                              <  1:2083972/93‑1 (MQ=255)
                                           aaaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                                              >  2:3445712/1‑93 (MQ=255)
                                           aaaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                                              >  2:8572447/1‑93 (MQ=255)
                                            aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAc                                             <  2:2342223/93‑1 (MQ=255)
                                            aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAc                                             >  1:8968481/1‑93 (MQ=255)
                                              gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAg                                           <  1:2624966/93‑1 (MQ=255)
                                               agTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgc                                          >  1:5096404/1‑93 (MQ=255)
                                               agTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgc                                          <  1:6226790/93‑1 (MQ=255)
                                                 tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                        <  2:10007971/93‑1 (MQ=255)
                                                  aaaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCt                                       <  2:6072361/93‑1 (MQ=255)
                                                  aaaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCt                                       <  2:5214834/93‑1 (MQ=255)
                                                   aaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtt                                      <  2:6518958/93‑1 (MQ=255)
                                                    aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                                     <  2:7006122/93‑1 (MQ=255)
                                                    aaGCTTCTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACACCGCttt                                     >  2:9661415/1‑93 (MQ=255)
                                                     aGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtttt                                    >  2:3149203/1‑93 (MQ=255)
                                                      gCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCTTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTa                                   <  1:3827959/93‑1 (MQ=255)
                                                      gCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTa                                   >  2:3973483/1‑93 (MQ=255)
                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 >  1:7943683/1‑93 (MQ=255)
                                                         tGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                <  1:9442097/93‑1 (MQ=255)
                                                            aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTcc                             >  1:4827280/1‑93 (MQ=255)
                                                            aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTcc                             <  1:5423074/93‑1 (MQ=255)
                                                             aGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCa                            <  1:8302773/93‑1 (MQ=255)
                                                             aGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCa                            <  1:2337858/93‑1 (MQ=255)
                                                                 gcgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca                        <  2:289502/93‑1 (MQ=255)
                                                                 gcgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca                        >  1:4381329/1‑93 (MQ=255)
                                                                   gcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACt                      >  2:1596194/1‑93 (MQ=255)
                                                                    ccaccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTa                     >  1:4560354/1‑93 (MQ=255)
                                                                      accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                   >  2:4900658/1‑93 (MQ=255)
                                                                       ccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTacac                  >  1:9931204/1‑93 (MQ=255)
                                                                       ccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTacac                  <  2:2237510/93‑1 (MQ=255)
                                                                       ccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTacac                  >  2:3915823/1‑93 (MQ=255)
                                                                       ccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTacac                  >  1:6357076/1‑93 (MQ=255)
                                                                        caAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAcc                 >  2:9824335/1‑93 (MQ=255)
                                                                         aaaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                >  1:4160425/1‑93 (MQ=255)
                                                                         aaaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                <  1:251964/93‑1 (MQ=255)
                                                                          aaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAATAGCGCTTTTACTTTCCAGACACTACACCCg               <  2:8852514/93‑1 (MQ=255)
                                                                          aaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCg               <  1:6496750/93‑1 (MQ=255)
                                                                            atcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGa             <  2:7261099/93‑1 (MQ=255)
                                                                             tcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGAt            <  1:4624837/93‑1 (MQ=255)
                                                                              catcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGAtc           >  1:1522310/1‑92 (MQ=255)
                                                                              catcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATa           <  1:3034740/93‑1 (MQ=255)
                                                                              catcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATa           >  2:1378336/1‑93 (MQ=255)
                                                                                tcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCTGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCGGGATagc         >  1:6067983/1‑91 (MQ=255)
                                                                                 cGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACca        <  2:2559869/93‑1 (MQ=255)
                                                                                  gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcat       >  2:5138849/1‑93 (MQ=255)
                                                                                  gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcat       <  2:4170923/93‑1 (MQ=255)
                                                                                    gAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatca     <  2:1888587/93‑1 (MQ=255)
                                                                                     aaaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatcat    <  2:3737790/93‑1 (MQ=255)
                                                                                            |                                                                                        
AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

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