Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,458,277 T→C V634A (GTG→GCG)  gspD → Type II secretion system protein GspD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,458,2770TC100.0% 923.3 / NA 226V634A (GTG→GCG) gspDType II secretion system protein GspD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (136/90);  total (136/90)

Only 100 of 226 total aligned reads displayed.

TCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCG  >  NC_000913/3458180‑3458374
                                                                                                 |                                                                                                 
  aCTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTgg                                                                                              >  1:5071678/1‑101 (MQ=255)
   cTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGa                                                                                             <  2:5071210/101‑1 (MQ=255)
    tGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGAt                                                                                            <  2:5071390/101‑1 (MQ=255)
     gTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATg                                                                                           <  2:5071113/101‑1 (MQ=255)
     gTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATg                                                                                           >  1:5071681/1‑101 (MQ=255)
       cAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGaa                                                                                         >  1:5071684/1‑101 (MQ=255)
       cAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGaa                                                                                         >  2:5071683/1‑101 (MQ=255)
       cAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGaa                                                                                         >  1:5071687/1‑101 (MQ=255)
       cAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGaa                                                                                         <  1:5071344/101‑1 (MQ=255)
        aaaaGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGaaa                                                                                        >  1:5071689/1‑101 (MQ=255)
        aaaaGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGaaa                                                                                        <  1:5071239/101‑1 (MQ=255)
        aaaaGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGaaa                                                                                        >  2:5071690/1‑101 (MQ=255)
         aaaGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGaaaa                                                                                       <  2:5071287/101‑1 (MQ=255)
          aaGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAac                                                                                      >  1:5071694/1‑101 (MQ=255)
           agagAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAaca                                                                                     >  1:5071696/1‑101 (MQ=255)
           agagAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAaca                                                                                     <  1:5071289/101‑1 (MQ=255)
           agagAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAaca                                                                                     <  2:5071365/101‑1 (MQ=255)
            gagaAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAacac                                                                                    <  2:5071159/101‑1 (MQ=255)
            gagaAATACACCCGTTACCGTCAGGAACAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAacac                                                                                    >  2:5071698/1‑101 (MQ=255)
               aaaTACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGtt                                                                                 >  1:5071701/1‑101 (MQ=255)
                aaTACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTc                                                                                >  2:5071705/1‑101 (MQ=255)
                 aTACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCa                                                                               >  2:5071708/1‑101 (MQ=255)
                    cacCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACa                                                                            >  2:5071711/1‑101 (MQ=255)
                     acCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAg                                                                           >  2:5071714/1‑101 (MQ=255)
                      cccGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGt                                                                          >  1:5071715/1‑101 (MQ=255)
                       ccGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTc                                                                                <  1:5071218/94‑1 (MQ=255)
                       ccGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTc                                                                         >  1:5071719/1‑101 (MQ=255)
                       ccGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTc                                                                         <  1:5071275/101‑1 (MQ=255)
                       ccGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTc                                                                         >  2:5071718/1‑101 (MQ=255)
                       ccGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTc                                                                         <  1:5071250/101‑1 (MQ=255)
                        cGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCa                                                                        <  1:5071265/101‑1 (MQ=255)
                         gTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCAc                                                                       <  2:5071135/101‑1 (MQ=255)
                          ttACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAg                                                                           >  2:5071725/1‑96 (MQ=255)
                          ttACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACg                                                                      >  1:5071724/1‑101 (MQ=255)
                          ttACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACg                                                                      >  2:5071722/1‑101 (MQ=255)
                           tACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGc                                                                     <  2:5071434/101‑1 (MQ=255)
                              cGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAg                                                                           <  2:5071138/92‑1 (MQ=255)
                              cGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGcccc                                                                  >  1:5071729/1‑101 (MQ=255)
                               gTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCt                                                                 >  2:5071732/1‑101 (MQ=255)
                                 cAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACg                                                                      >  2:5071737/1‑94 (MQ=255)
                                 cAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTgc                                                               >  1:5071733/1‑101 (MQ=255)
                                 cAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTgc                                                               >  1:5071735/1‑101 (MQ=255)
                                  aGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCt                                                                 <  1:5071101/98‑1 (MQ=255)
                                  aGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTgcg                                                              >  2:5071739/1‑101 (MQ=255)
                                   ggAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTgcgc                                                             >  1:5071742/1‑101 (MQ=255)
                                      gCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCAt                                                          <  2:5071598/101‑1 (MQ=255)
                                      gCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCAt                                                          <  1:5071311/101‑1 (MQ=255)
                                       caacaaCAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATc                                                         >  2:5071743/1‑101 (MQ=255)
                                         acaacaGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGt                                                       >  2:5071748/1‑101 (MQ=255)
                                         acaacaGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGt                                                       <  2:5071184/101‑1 (MQ=255)
                                         acaacaGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGt                                                       >  2:5071746/1‑101 (MQ=255)
                                           aacaGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTc                                                      >  2:5071752/1‑100 (MQ=255)
                                           aacaGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCa                                                     >  1:5071750/1‑101 (MQ=255)
                                            acaGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGaaa                                                                                        <  2:5071226/65‑1 (MQ=255)
                                            acaGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCAc                                                    <  1:5070975/101‑1 (MQ=255)
                                              aGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACgg                                                  <  2:5071353/101‑1 (MQ=255)
                                                 ggATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGa                                               >  2:5071756/1‑101 (MQ=255)
                                                   aTCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAgg                                             >  1:5071757/1‑101 (MQ=255)
                                                        cGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACgg                                                  >  2:5071761/1‑91 (MQ=255)
                                                         gggAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCa                                       >  2:5071763/1‑101 (MQ=255)
                                                          ggAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCat                                      >  1:5071764/1‑101 (MQ=255)
                                                            aaaTCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCt                                                                 <  2:5071228/72‑1 (MQ=255)
                                                            aaaTCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCatat                                    >  2:5071766/1‑101 (MQ=255)
                                                             aaTCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATg                                   <  1:5071428/101‑1 (MQ=255)
                                                              aTCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATga                                  >  1:5071768/1‑101 (MQ=255)
                                                               tCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATgag                                 >  2:5071771/1‑101 (MQ=255)
                                                               tCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATgag                                 >  1:5071772/1‑101 (MQ=255)
                                                                cAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATgaga                                >  2:5071775/1‑101 (MQ=255)
                                                                  aaaGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAAt                              >  2:5071780/1‑101 (MQ=255)
                                                                  aaaGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAAt                              >  1:5071779/1‑101 (MQ=255)
                                                                   aaGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAAtt                             >  1:5071782/1‑101 (MQ=255)
                                                                    aGCGCTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTc                                                      >  1:5071784/1‑75 (MQ=255)
                                                                     gcgcTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATtca                           <  2:5071357/101‑1 (MQ=255)
                                                                     gcgcTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATtca                           <  1:5071141/101‑1 (MQ=255)
                                                                      cgcTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATtcac                          >  2:5071787/1‑101 (MQ=255)
                                                                       gcTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCat                                      <  2:5071343/88‑1 (MQ=255)
                                                                       gcTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATtcact                         >  1:5071789/1‑101 (MQ=255)
                                                                        cTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATtcactc                        <  1:5071482/101‑1 (MQ=255)
                                                                        cTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATtcactc                        <  1:5071368/101‑1 (MQ=255)
                                                                        cTGGTTGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATtcactc                        >  2:5071790/1‑101 (MQ=255)
                                                                             tGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGt                   <  2:5071304/101‑1 (MQ=255)
                                                                             tGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGt                   >  2:5071796/1‑101 (MQ=255)
                                                                             tGGCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGt                   >  1:5071795/1‑101 (MQ=255)
                                                                              ggCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTa                  <  2:5071266/101‑1 (MQ=255)
                                                                               gCTCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTAc                 >  1:5071798/1‑101 (MQ=255)
                                                                                 tCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATtc                            <  2:5071270/88‑1 (MQ=255)
                                                                                 tCGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTAccc               >  1:5071800/1‑101 (MQ=255)
                                                                                  cGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTAcccc              <  1:5071117/101‑1 (MQ=255)
                                                                                  cGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTAcccc              <  1:5071318/101‑1 (MQ=255)
                                                                                  cGGAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTAcccc              <  2:5071542/101‑1 (MQ=255)
                                                                                   ggAAGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCg             >  2:5071803/1‑101 (MQ=255)
                                                                                     aaGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGcc           <  1:5071556/101‑1 (MQ=255)
                                                                                      aGATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCa          >  1:5071805/1‑101 (MQ=255)
                                                                                       gATTTGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAg         >  1:5071808/1‑101 (MQ=255)
                                                                                          ttGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTg      <  1:5071161/101‑1 (MQ=255)
                                                                                           tGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTgg     <  2:5071487/101‑1 (MQ=255)
                                                                                           tGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTgg     >  1:5071810/1‑101 (MQ=255)
                                                                                           tGCCGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTgg     >  2:5071811/1‑101 (MQ=255)
                                                                                             ccGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCAtt                                         <  2:5071240/63‑1 (MQ=255)
                                                                                              cGGCGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGt                                                 <  1:5071182/54‑1 (MQ=255)
                                                                                                 |                                                                                                 
TCACTGTCAAAAGAGAAATACACCCGTTACCGTCAGGAGCAACAACAGCGGATCGACGGGAAATCAAAAGCGCTGGTTGGCTCGGAAGATTTGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCG  >  NC_000913/3458180‑3458374

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: ‑