Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,113,979 (C)5→4 coding (374/810 nt) pppA ← prepilin peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,113,9750C.100.0% 1514.3 / NA 289coding (378/810 nt)pppAprepilin peptidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (154/135);  total (154/135)

Only 100 of 311 total aligned reads displayed.

TCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGA  >  NC_000913/3113875‑3114071
                                                                                                    |                                                                                                
tCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGc                                                                                                  <  2:5248759/101‑1 (MQ=255)
 ccACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGcc                                                                                                 >  2:5249478/1‑101 (MQ=255)
  cacaATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGccc                                                                                                <  1:5248930/101‑1 (MQ=255)
   acaATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGcccc                                                                                               >  1:5249480/1‑101 (MQ=255)
     aaTACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGccccat                                                                                             >  1:5249483/1‑99 (MQ=255)
     aaTACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGccccat                                                                                             <  1:5248981/101‑3 (MQ=255)
      aTACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGccccatc                                                                                            >  2:5249487/1‑98 (MQ=255)
       tACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGccccatcc                                                                                           >  2:5249491/1‑97 (MQ=255)
        aCGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCa                                                                                         >  1:5249492/1‑101 (MQ=255)
          gCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACt                                                                                       >  1:5249497/1‑101 (MQ=255)
             cTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTc                                                                                    >  2:5249501/1‑101 (MQ=255)
              tGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCt                                                                                   >  2:5249504/1‑101 (MQ=255)
              tGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCt                                                                                   >  1:5249507/1‑101 (MQ=255)
               gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCt                                                                                   <  2:5248862/100‑1 (MQ=255)
               gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTg                                                                                  >  2:5249512/1‑101 (MQ=255)
               gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTg                                                                                  <  1:5249015/101‑1 (MQ=255)
               gAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTg                                                                                  >  1:5249510/1‑101 (MQ=255)
                  tAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGcc                                                                               <  2:5249112/101‑1 (MQ=255)
                   aaaaaCATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGcca                                                                              >  1:5249516/1‑101 (MQ=255)
                   aaaaaCATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGcca                                                                              >  1:5249519/1‑101 (MQ=255)
                    aaaaCATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGccag                                                                             >  1:5249523/1‑101 (MQ=255)
                     aaaCATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGccaga                                                                            <  1:5248748/101‑1 (MQ=255)
                     aaaCATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGccaga                                                                            >  2:5249524/1‑101 (MQ=255)
                         aTCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGCTAATATCATCACCGCCAGCCCCATCCACTTTCTGGccaga                                                                            >  1:5249533/1‑97 (MQ=255)
                         aTCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGccagaccag                                                                        >  2:5249529/1‑101 (MQ=255)
                         aTCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGccagaccag                                                                        <  1:5249038/101‑1 (MQ=255)
                         aTCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGccagaccag                                                                        >  1:5249530/1‑101 (MQ=255)
                            gggCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACt                                                                     >  2:5249536/1‑101 (MQ=255)
                                 gCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCa                                                                >  1:5249541/1‑101 (MQ=255)
                                   cATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGc                                                              <  1:5248915/101‑1 (MQ=255)
                                    aTTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCa                                                             >  2:5249544/1‑101 (MQ=255)
                                    aTTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCa                                                             <  1:5248950/101‑1 (MQ=255)
                                    aTTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCa                                                             >  1:5249546/1‑101 (MQ=255)
                                     ttGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaa                                                            <  2:5248984/101‑1 (MQ=255)
                                     ttGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaa                                                            >  1:5249551/1‑101 (MQ=255)
                                      tGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaaa                                                           <  1:5249084/101‑1 (MQ=255)
                                        gTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaaaaa                                                         <  1:5248640/101‑1 (MQ=255)
                                        gTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaaaaa                                                         >  1:5249554/1‑101 (MQ=255)
                                        gTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaaaaa                                                         <  2:5249222/101‑1 (MQ=255)
                                         tGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaaaaaa                                                        <  1:5248999/101‑1 (MQ=255)
                                           aTCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGc                                                      >  1:5249558/1‑101 (MQ=255)
                                           aTCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGc                                                      <  2:5248965/101‑1 (MQ=255)
                                            tCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCa                                                     >  2:5249561/1‑101 (MQ=255)
                                              gAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAg                                                   >  1:5249564/1‑101 (MQ=255)
                                               aGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGccc                                                                                                <  2:5248913/56‑1 (MQ=255)
                                               aGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGt                                                  <  2:5249161/101‑1 (MQ=255)
                                               aGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGt                                                  <  2:5249008/101‑1 (MQ=255)
                                                ggTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGc                                                                                                  >  2:5249568/1‑53 (MQ=255)
                                                ggTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCACAAAAGCAAGTg                                                 <  1:5249131/101‑1 (MQ=255)
                                                 gTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGc                                                <  2:5249037/101‑1 (MQ=255)
                                                  tCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCt                                               >  1:5249569/1‑101 (MQ=255)
                                                   cGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTg                                              >  1:5249574/1‑101 (MQ=255)
                                                   cGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTg                                              <  2:5249049/101‑1 (MQ=255)
                                                   cGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTg                                              <  1:5248809/101‑1 (MQ=255)
                                                    gATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGt                                             >  1:5249577/1‑101 (MQ=255)
                                                    gATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGt                                             <  1:5249234/101‑1 (MQ=255)
                                                     aTGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTc                                            >  1:5249581/1‑101 (MQ=255)
                                                      tGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCa                                           >  2:5249583/1‑101 (MQ=255)
                                                        acacTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAg                                                               <  2:5248842/79‑1 (MQ=255)
                                                        acacTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAAt                                         <  1:5249221/101‑1 (MQ=255)
                                                         cacTCGCGGCAATCAGCCAGGCGGATAATATCATCGCCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATa                                        <  1:5249103/101‑1 (MQ=255)
                                                         cacTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATa                                        >  1:5249585/1‑101 (MQ=255)
                                                          acTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATaa                                       <  2:5249146/101‑1 (MQ=255)
                                                             cgcgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTc                                    <  1:5249366/101‑1 (MQ=255)
                                                               cgGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCAc                                                                                        <  2:5249064/47‑1 (MQ=39)
                                                                ggCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTccac                                 >  1:5249591/1‑101 (MQ=255)
                                                                 gCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTccacc                                <  2:5249298/101‑1 (MQ=255)
                                                                  cAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTccacca                               >  2:5249592/1‑101 (MQ=255)
                                                                  cAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTccacca                               <  2:5248816/101‑1 (MQ=255)
                                                                   aaTCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAg                              >  1:5249597/1‑101 (MQ=255)
                                                                    aTCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGc                             >  1:5249600/1‑101 (MQ=255)
                                                                    aTCAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGc                             <  1:5249164/101‑1 (MQ=255)
                                                                      cAGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCAACCAGCgg                           >  2:5249602/1‑101 (MQ=255)
                                                                       aGCCAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGc                             >  1:5249606/1‑98 (MQ=255)
                                                                         ccAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATa                        <  1:5249054/101‑1 (MQ=255)
                                                                         ccAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATa                        >  2:5249609/1‑101 (MQ=255)
                                                                          cAGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATaa                       <  1:5249122/101‑1 (MQ=255)
                                                                           aGGCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAAc                      >  1:5249611/1‑101 (MQ=255)
                                                                            ggCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACg                     <  2:5249258/101‑1 (MQ=255)
                                                                            ggCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACg                     >  1:5249614/1‑101 (MQ=255)
                                                                             gCGGATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGc                    >  2:5249617/1‑101 (MQ=255)
                                                                               ggATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCtt                  >  2:5249621/1‑101 (MQ=255)
                                                                               ggATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACACCACCAGCGGATAACGCtt                  >  1:5249622/1‑101 (MQ=255)
                                                                                gATAATATCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTg                 <  1:5248962/101‑1 (MQ=255)
                                                                                 ataataTCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGc                >  2:5249626/1‑101 (MQ=255)
                                                                                   gatatCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTg              >  1:5249629/2‑101 (MQ=255)
                                                                                     tatCATCACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAATGCTTGCTGAt            >  2:5249633/1‑101 (MQ=255)
                                                                                      atcatcACCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtt           >  1:5249634/1‑101 (MQ=255)
                                                                                        catcaCCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtttt         >  2:5249638/1‑101 (MQ=255)
                                                                                        catcaCCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtttt         <  1:5249143/101‑1 (MQ=255)
                                                                                        catcaCCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGAtttt         <  2:5249024/101‑1 (MQ=255)
                                                                                         atcaCCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTc        <  1:5248855/101‑1 (MQ=255)
                                                                                          tctccGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGAt                         <  1:5249065/80‑1 (MQ=255)
                                                                                          tcaCCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGAt                         <  2:5249199/83‑1 (MQ=255)
                                                                                          tcaCCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCg       <  1:5249342/101‑1 (MQ=255)
                                                                                          tcaCCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCg       >  1:5249643/1‑101 (MQ=255)
                                                                                           caCCGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCaaaa                                                          <  2:5248790/49‑1 (MQ=39)
                                                                                             ccGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGc      <  2:5249078/99‑1 (MQ=255)
                                                                                             ccGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCt    >  1:5249646/1‑101 (MQ=255)
                                                                                             ccGCCAGCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCt    <  2:5249105/101‑1 (MQ=255)
                                                                                                    |                                                                                                
TCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCAGCCCCCATCCACTTTCTGGCCAGACCAGACTCGCCAGCAAAAAAGCAAGTGCTGTCAATAACTCCACCAGCGGATAACGCTTGCTGATTTTCGCCTGA  >  NC_000913/3113875‑3114071

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