Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,593,254 2 bp→GC coding (59‑60/807 nt) yjjM ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,593,2540CG100.0% 936.2 / NA 286A20G (GCG→GGCyjjMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (157/129);  total (157/129)
*NC_0009134,593,2550GC100.0% 933.0 / NA 286A20G (GCG→GGC) yjjMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (157/129);  total (157/129)

Only 100 of 290 total aligned reads displayed.

ACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAG  >  NC_000913/4593167‑4593339
                                                                                       ||                                                                                    
aCTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCa                                                                                  >  1:4882552/1‑93 (MQ=255)
aCTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCa                                                                                  >  1:5895827/1‑93 (MQ=255)
  tCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCAtt                                                                                >  1:3232856/1‑93 (MQ=255)
   cATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTc                                                                               <  1:1285595/93‑1 (MQ=255)
    aTTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcg                                                                              >  2:1314336/1‑93 (MQ=255)
    aTTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcg                                                                              >  1:6460448/1‑93 (MQ=255)
    aTTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcg                                                                              <  1:6418459/93‑1 (MQ=255)
    aTTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcg                                                                              >  2:5338588/1‑93 (MQ=255)
     ttGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcgc                                                                             >  1:610523/1‑93 (MQ=255)
          gCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAg                                                                        <  2:3097588/93‑1 (MQ=255)
          gCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAg                                                                        <  1:6728039/93‑1 (MQ=255)
          gCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAg                                                                        >  2:5584787/1‑93 (MQ=255)
           cGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGt                                                                       >  2:4535710/1‑93 (MQ=255)
           cGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGt                                                                       <  1:1259155/93‑1 (MQ=255)
            ggTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTg                                                                      >  2:3187091/1‑93 (MQ=255)
             gtagtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTgg                                                                     >  1:5239253/1‑93 (MQ=255)
              tagtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGc                                                                    >  1:3089459/1‑93 (MQ=255)
              tagtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGc                                                                    >  2:4348617/1‑93 (MQ=255)
               agtaCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCt                                                                   <  1:4233255/93‑1 (MQ=255)
                 taCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTga                                                                 <  1:4625080/93‑1 (MQ=255)
                 taCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTga                                                                 <  1:4842891/93‑1 (MQ=255)
                  aCAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTgag                                                                >  1:4930858/1‑93 (MQ=255)
                   cAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTgaga                                                               >  2:4755278/1‑93 (MQ=255)
                    aGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGaa                                                              <  1:4727074/93‑1 (MQ=255)
                    aGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGaa                                                              >  2:3013333/1‑93 (MQ=255)
                     ggCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAat                                                             >  2:2633715/1‑93 (MQ=255)
                     ggCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAat                                                             >  2:1728454/1‑93 (MQ=255)
                      gCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAata                                                            <  2:3926840/93‑1 (MQ=255)
                        gAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATg                                                          >  1:514377/1‑93 (MQ=255)
                         aaaCCGTCATCATGGTCAGGTTTACGCGCAATAGCATAGTCGTTCCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGa                                                         <  1:880067/93‑1 (MQ=255)
                         aaaCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGa                                                         <  2:25054/93‑1 (MQ=255)
                          aaCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGAc                                                        <  1:2645739/93‑1 (MQ=255)
                           aCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACg                                                       >  1:4771904/1‑93 (MQ=255)
                            ccGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGc                                                      <  2:3378789/93‑1 (MQ=255)
                            ccGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGc                                                      >  2:6597046/1‑93 (MQ=255)
                            ccGTCATCATGGTCAGGGTTACGCGCAATAACATAGTCGTTTCCAACCTGCGGCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGc                                                      >  2:1058588/1‑93 (MQ=255)
                             cGTCATCATGGTCAGGTGTACGTGCAATAACATAGTCGTTTCCAAAATCAGTAAAGACGCCGCATTCGCGTCAGAGTGTGAGAATAGGACGCa                                                     >  2:6285428/1‑93 (MQ=255)
                              gTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCAc                                                    >  2:2102431/1‑93 (MQ=255)
                                 atcatGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtgt                                                 <  1:2507325/93‑1 (MQ=255)
                                 atcatGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtgt                                                 <  2:6530169/93‑1 (MQ=255)
                                   catGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtgg                                               <  1:3684429/93‑1 (MQ=255)
                                   catGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtgg                                               >  1:2147795/1‑93 (MQ=255)
                                    atGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCATCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggt                                              <  1:6625288/93‑1 (MQ=255)
                                     tGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggtg                                             <  1:335811/93‑1 (MQ=255)
                                      ggTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGc                                            >  2:4602000/1‑93 (MQ=255)
                                      ggTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGc                                            <  2:220141/93‑1 (MQ=255)
                                       gTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCg                                           <  1:1385656/93‑1 (MQ=255)
                                       gTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCg                                           <  2:786730/93‑1 (MQ=255)
                                        tCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCgg                                          >  1:2409549/1‑93 (MQ=255)
                                          aGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCt                                        <  2:111154/93‑1 (MQ=255)
                                          aGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCt                                        >  2:1045479/1‑93 (MQ=255)
                                          aGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCt                                        >  1:1157371/1‑93 (MQ=255)
                                           ggTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTg                                       <  1:1040916/93‑1 (MQ=255)
                                           ggTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTg                                       <  1:5722284/93‑1 (MQ=255)
                                           ggTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGAAGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTg                                       >  1:2960537/1‑93 (MQ=255)
                                            gTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGa                                      >  2:4872621/1‑93 (MQ=255)
                                               tACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgt                                   >  2:2392702/1‑93 (MQ=255)
                                                aCGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgtt                                  >  2:6135688/1‑93 (MQ=255)
                                                aCGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgtt                                  <  2:2184088/93‑1 (MQ=255)
                                                 cgcgCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttg                                 <  1:224247/93‑1 (MQ=255)
                                                  gcgcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttgt                                <  2:5804060/93‑1 (MQ=255)
                                                     cAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACa                             >  2:3589526/1‑93 (MQ=255)
                                                      aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTg                                       <  1:3440575/82‑1 (MQ=255)
                                                      aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAt                            <  1:4350940/93‑1 (MQ=255)
                                                      aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAt                            >  2:5507689/1‑93 (MQ=255)
                                                      aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAt                            <  2:32925/93‑1 (MQ=255)
                                                      aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATACGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAt                            <  1:1595892/93‑1 (MQ=255)
                                                       aTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAtt                           <  2:1650218/93‑1 (MQ=255)
                                                       aTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAtt                           <  2:509219/93‑1 (MQ=255)
                                                        tAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAttt                          >  2:6468763/1‑93 (MQ=255)
                                                          aCATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTcc                        <  2:5527692/93‑1 (MQ=255)
                                                          aCATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTcc                        >  2:4085147/1‑93 (MQ=255)
                                                           cATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCg                       >  2:1434272/1‑93 (MQ=255)
                                                            aTAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGc                      >  2:5975400/1‑93 (MQ=255)
                                                            aTAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGc                      >  2:3472911/1‑93 (MQ=255)
                                                              aGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCa                    >  1:732571/1‑93 (MQ=255)
                                                               gTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAg                   <  1:5484372/93‑1 (MQ=255)
                                                               gTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAg                   <  1:1684312/93‑1 (MQ=255)
                                                                tCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgc                  <  2:3092169/93‑1 (MQ=255)
                                                                tCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAgc                  >  1:377557/1‑93 (MQ=255)
                                                                   tttCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCa               <  1:901373/93‑1 (MQ=255)
                                                                   tttCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCa               >  1:5086218/1‑93 (MQ=255)
                                                                      ccAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACtt            <  1:3843233/93‑1 (MQ=255)
                                                                      ccAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACtt            <  1:5434760/93‑1 (MQ=255)
                                                                      ccAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACtt            >  1:4504039/1‑93 (MQ=255)
                                                                       cAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACttt           >  2:1186060/1‑93 (MQ=255)
                                                                       cAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACttt           >  2:4427094/1‑93 (MQ=255)
                                                                        aaCCGGCGTCAAGACGCCGCATTCGCGGAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTg          >  2:6047055/1‑93 (MQ=255)
                                                                         aCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGc         >  2:3849975/1‑93 (MQ=255)
                                                                         aCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGc         >  2:2741184/1‑93 (MQ=255)
                                                                         aCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGc         >  1:5927834/1‑93 (MQ=255)
                                                                         aCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGc         >  2:1381349/1‑93 (MQ=255)
                                                                          ccTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCg        >  1:3987291/1‑93 (MQ=255)
                                                                           cTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGa       <  2:5745134/93‑1 (MQ=255)
                                                                            tGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGAt      <  1:6140318/93‑1 (MQ=255)
                                                                               gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga   >  1:2558678/1‑91 (MQ=255)
                                                                               gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga   >  2:5672418/1‑91 (MQ=255)
                                                                               gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATgga   >  1:5608598/1‑91 (MQ=255)
                                                                                tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa  <  1:1808741/93‑4 (MQ=255)
                                                                                tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATggaa  >  2:4977524/1‑90 (MQ=255)
                                                                                       ||                                                                                    
ACTCATTGAAGCGGTAGTACAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTTGCGATGAAG  >  NC_000913/4593167‑4593339

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