Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,972,637 2 bp→GC coding (430‑431/546 nt) wecD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6370CG100.0% 439.8 / NA 135R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (63/72);  total (63/72)
*NC_0009133,972,6380GC100.0% 440.2 / NA 135R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (63/72);  total (63/72)

Only 100 of 136 total aligned reads displayed.

AATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTA  >  NC_000913/3972553‑3972723
                                                                                    ||                                                                                     
aaTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGt                                                                                >  2:837228/1‑93 (MQ=255)
  tGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTaa                                                                              >  2:77711/1‑93 (MQ=255)
    cGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTaaaa                                                                            <  2:4377821/93‑1 (MQ=255)
     gACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaac                                                                           <  2:713934/93‑1 (MQ=255)
     gACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaac                                                                           >  2:714933/1‑93 (MQ=255)
      aCAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaaca                                                                          <  2:1893752/93‑1 (MQ=255)
       cAGATGCGCGAATTGGCCTGCTGGATGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaa                                                                         <  2:2068746/93‑1 (MQ=255)
         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                       <  1:2419382/93‑1 (MQ=255)
          aTGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACtt                                                                      <  2:2279567/93‑1 (MQ=255)
           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                     <  2:3087503/93‑1 (MQ=255)
           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                     <  1:732681/93‑1 (MQ=255)
             cgcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGc                                                                   >  1:1727031/1‑93 (MQ=255)
              gcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCg                                                                  <  1:3545532/93‑1 (MQ=255)
              gcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCg                                                                  <  1:4412953/93‑1 (MQ=255)
               cgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCgg                                                                 <  1:3817775/93‑1 (MQ=255)
               cgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCgg                                                                 <  2:2943874/93‑1 (MQ=255)
               cgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCgg                                                                 <  2:1768405/93‑1 (MQ=255)
               cgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCgg                                                                 <  1:2682541/93‑1 (MQ=255)
                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                               <  1:1007802/93‑1 (MQ=255)
                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAAAAACTTTGCGGGt                                                               >  1:2487752/1‑93 (MQ=255)
                  aTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTg                                                              >  1:2399090/1‑93 (MQ=255)
                   ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTgg                                                             >  2:1837520/1‑93 (MQ=255)
                    tGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGc                                                            <  2:3977005/93‑1 (MQ=255)
                      gCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGa                                                          <  1:1923853/93‑1 (MQ=255)
                       ccTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGAc                                                         >  1:1562711/1‑93 (MQ=255)
                        ctgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAcc                                                        >  2:2142458/1‑93 (MQ=255)
                         tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAcca                                                       >  2:4082498/1‑92 (MQ=255)
                          gctggctgGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCa                                                      >  1:3696565/1‑93 (MQ=255)
                          gctggctgGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCa                                                      <  2:1962240/93‑1 (MQ=255)
                            tggctggACGCGGTGCAGGTGCTGAGCTGATGCGAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCaaa                                                    <  2:1889880/93‑1 (MQ=255)
                            tggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCaaa                                                    >  2:4694655/1‑93 (MQ=255)
                             ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAAt                                                   <  2:4102920/93‑1 (MQ=255)
                             ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAAt                                                   >  2:634276/1‑93 (MQ=255)
                               ctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATgg                                                 <  2:2492885/93‑1 (MQ=255)
                                tggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATggg                                                <  1:329435/93‑1 (MQ=255)
                                tggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATggg                                                <  2:1900877/93‑1 (MQ=255)
                                tggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATggg                                                <  2:3656465/93‑1 (MQ=255)
                                 ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                               >  1:3524649/1‑93 (MQ=255)
                                 ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                               >  1:366312/1‑93 (MQ=255)
                                  ggcgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCa                                              <  1:1978165/91‑1 (MQ=255)
                                  gACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCa                                              >  2:4166534/1‑93 (MQ=255)
                                   gcgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCaa                                             <  1:1776425/92‑1 (MQ=255)
                                   aCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCaa                                             <  2:1065641/93‑1 (MQ=255)
                                    cgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAac                                            <  1:1964789/93‑1 (MQ=255)
                                    agcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAac                                            <  1:1246544/92‑1 (MQ=255)
                                     gcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAaca                                           >  1:2164794/1‑93 (MQ=255)
                                      cgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAacac                                          <  1:1026794/93‑1 (MQ=255)
                                        gTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                                        >  1:2544178/1‑93 (MQ=255)
                                        gTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAAAAACTTTGCGGGGGGCGACCCAAATGGGCAAAAccc                                        >  2:900069/1‑92 (MQ=255)
                                         tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                       <  1:496393/93‑1 (MQ=255)
                                           cAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgccg                                     <  1:1782699/93‑1 (MQ=255)
                                           cAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgccg                                     <  2:4678687/93‑1 (MQ=255)
                                            aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                    >  2:3655744/1‑93 (MQ=255)
                                            aGGTGCTGAGCTGATGCAAACGGCACTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGAGAGGCGCGACCCAAATGGGAAACACggacgc                                    >  1:1685225/1‑93 (MQ=255)
                                             ggTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCcgcg                                   >  2:2817587/1‑93 (MQ=255)
                                              gTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCgcgc                                  <  1:4487126/93‑1 (MQ=255)
                                                gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt                                <  1:1291900/93‑1 (MQ=255)
                                                gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt                                <  1:1194088/93‑1 (MQ=255)
                                                gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt                                >  1:241585/1‑93 (MQ=255)
                                                  tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa                              >  2:123774/1‑93 (MQ=255)
                                                  tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa                              >  1:2954099/1‑93 (MQ=255)
                                                   gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                             <  1:3511545/93‑1 (MQ=255)
                                                     gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                           >  1:3745711/1‑93 (MQ=255)
                                                     gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                           <  2:941538/93‑1 (MQ=255)
                                                       tgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGat                         >  2:661486/1‑93 (MQ=255)
                                                       tgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGat                         >  1:3446630/1‑93 (MQ=255)
                                                       tgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGat                         >  1:2401629/1‑93 (MQ=255)
                                                          tGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGataca                      <  1:1034647/93‑1 (MQ=255)
                                                            cAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                    >  1:2230781/1‑93 (MQ=255)
                                                             aaaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                   >  2:70231/1‑93 (MQ=255)
                                                               aCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaa                 >  2:4727736/1‑93 (MQ=255)
                                                               aCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaa                 >  1:4714746/1‑93 (MQ=255)
                                                                cGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaaa                <  2:2766158/93‑1 (MQ=255)
                                                                 ggCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg               <  2:51875/93‑1 (MQ=255)
                                                                 ggCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg               >  1:2957766/1‑93 (MQ=255)
                                                                 ggCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg               >  1:3920485/1‑93 (MQ=255)
                                                                   cgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg             <  2:3036209/93‑1 (MQ=255)
                                                                   cgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg             <  1:1935961/93‑1 (MQ=255)
                                                                   cgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg             <  2:579285/93‑1 (MQ=255)
                                                                    gcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg            >  2:2525389/1‑93 (MQ=255)
                                                                    gcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg            >  1:428457/1‑93 (MQ=255)
                                                                     cTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt           >  1:1719463/1‑93 (MQ=255)
                                                                      tAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgaggtg          >  2:3919287/1‑93 (MQ=255)
                                                                      tAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggtg          <  1:1548714/93‑1 (MQ=255)
                                                                      tAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggtg          <  1:1548713/93‑1 (MQ=255)
                                                                       aaaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc         <  1:3328796/93‑1 (MQ=255)
                                                                       aaaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc         <  1:3328913/93‑1 (MQ=255)
                                                                       aaaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc         <  2:1253443/93‑1 (MQ=255)
                                                                        aaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCg        >  1:4499983/1‑93 (MQ=255)
                                                                         aCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa       >  2:3370768/1‑93 (MQ=255)
                                                                         aCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa       >  1:2398290/1‑93 (MQ=255)
                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACATCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa      >  2:4593967/1‑93 (MQ=255)
                                                                           tGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt     >  1:3511649/1‑93 (MQ=255)
                                                                           tGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt     >  1:3283797/1‑93 (MQ=255)
                                                                            gggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg    <  1:930812/93‑1 (MQ=255)
                                                                            gggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg    >  2:943851/1‑93 (MQ=255)
                                                                             ggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAACGTGGTGCGAATGt   <  1:3136308/93‑1 (MQ=255)
                                                                              gCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa  <  1:737618/93‑1 (MQ=255)
                                                                              gCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa  <  1:3732674/93‑1 (MQ=255)
                                                                              gCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa  <  1:2508901/93‑1 (MQ=255)
                                                                                    ||                                                                                     
AATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTA  >  NC_000913/3972553‑3972723

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