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% 675.3 / NA 195R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (97/98);  total (97/98)
*NC_0009133,972,6380GC100.0% 673.2 / NA 194R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (96/98);  total (96/98)

Only 100 of 198 total aligned reads displayed.

CTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAG  >  NC_000913/3972550‑3972724
                                                                                       ||                                                                                      
cTCACTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTcgc                                                                                    <  1:4612079/93‑1 (MQ=255)
 tCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTcgcg                                                                                   >  1:1065475/1‑93 (MQ=255)
  cAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCgg                                                                                  >  2:4453457/1‑93 (MQ=255)
  cAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCgg                                                                                  >  1:3279992/1‑93 (MQ=255)
   aaTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCTCGGt                                                                                 >  2:888687/1‑93 (MQ=255)
   aaTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGt                                                                                 <  2:5981487/93‑1 (MQ=255)
    aTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTa                                                                                <  2:4946253/93‑1 (MQ=255)
    aTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTa                                                                                >  2:3651575/1‑93 (MQ=255)
    aTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTa                                                                                <  2:456071/93‑1 (MQ=255)
       cGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTaaaa                                                                             <  1:2329901/93‑1 (MQ=255)
       cGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTaaaa                                                                             <  2:1772491/93‑1 (MQ=255)
       cGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTaaaa                                                                             <  1:2970220/93‑1 (MQ=255)
         aCAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaaca                                                                           >  1:5918350/1‑93 (MQ=255)
          cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaa                                                                          >  2:2955908/1‑93 (MQ=255)
          cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaa                                                                          <  1:5611383/93‑1 (MQ=255)
           aGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaac                                                                         >  2:3286422/1‑93 (MQ=255)
           aGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaac                                                                         >  1:2246725/1‑93 (MQ=255)
            gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACt                                                                        >  2:5483285/1‑93 (MQ=255)
              tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACttt                                                                      <  1:915655/93‑1 (MQ=255)
                  cgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCgg                                                                  <  2:5363933/93‑1 (MQ=255)
                  cgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCgg                                                                  >  1:370944/1‑93 (MQ=255)
                    aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCgggg                                                                >  2:733872/1‑92 (MQ=255)
                    aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCgggg                                                                >  1:5787809/1‑92 (MQ=255)
                    aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGt                                                                <  2:4724679/93‑1 (MQ=255)
                    aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGGATGCTCGCGGTCAAACAACTTTGCgggg                                                                >  1:3378955/1‑92 (MQ=255)
                       tGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGc                                                             >  2:3967555/1‑93 (MQ=255)
                       tGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGc                                                             <  2:4164931/93‑1 (MQ=255)
                       tGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGc                                                             <  2:4129842/93‑1 (MQ=255)
                        ggCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCCGGTGGCg                                                            >  1:4409646/1‑93 (MQ=255)
                         gCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAAAAACTTTGCggggggggc                                                           >  2:4991932/1‑90 (MQ=255)
                           ctgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGAcc                                                         >  2:2311762/1‑93 (MQ=255)
                            tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGAccc                                                        >  1:2628513/1‑93 (MQ=255)
                            tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGAccc                                                        >  1:1085424/1‑93 (MQ=255)
                             gctggctgGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCa                                                       >  1:1943568/1‑93 (MQ=255)
                             gctggctgGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCa                                                       >  1:327599/1‑93 (MQ=255)
                                ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGTGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAAt                                                    <  1:3006797/93‑1 (MQ=255)
                                 gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                   <  2:2420455/93‑1 (MQ=255)
                                  ctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATgg                                                  <  2:2332511/93‑1 (MQ=255)
                                    ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                                >  1:409471/1‑93 (MQ=255)
                                    ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGCTGGCGACCCAAATGGGc                                                <  1:2767650/93‑1 (MQ=255)
                                       cgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAac                                             <  2:3560777/93‑1 (MQ=255)
                                        gcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAaca                                            <  1:1424794/93‑1 (MQ=255)
                                         cgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAacac                                           <  2:1724729/93‑1 (MQ=255)
                                          ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                          >  2:3321515/1‑93 (MQ=255)
                                            tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        <  1:3573646/93‑1 (MQ=255)
                                            tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        <  1:3980058/93‑1 (MQ=255)
                                              cAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgccg                                      <  2:4776031/93‑1 (MQ=255)
                                              cAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgccg                                      <  2:1305372/93‑1 (MQ=255)
                                               aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                     <  1:5833188/93‑1 (MQ=255)
                                               aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                     >  2:2420646/1‑93 (MQ=255)
                                                 gTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCgcgc                                   <  2:1619524/93‑1 (MQ=255)
                                                 gTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCgcgc                                   <  2:4064151/93‑1 (MQ=255)
                                                  tGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                                  <  2:2035486/93‑1 (MQ=255)
                                                   gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt                                 <  1:5551633/93‑1 (MQ=255)
                                                    ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                <  1:4151688/93‑1 (MQ=255)
                                                    ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                >  1:1555883/1‑93 (MQ=255)
                                                    ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                <  2:5927406/93‑1 (MQ=255)
                                                     tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa                               <  1:2285948/93‑1 (MQ=255)
                                                      gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                              >  2:3851133/1‑93 (MQ=255)
                                                      gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                              >  2:6256537/1‑93 (MQ=255)
                                                       agctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc                             >  2:3458728/1‑93 (MQ=255)
                                                       agctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc                             <  1:3857638/93‑1 (MQ=255)
                                                       agctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc                             <  2:3199437/93‑1 (MQ=255)
                                                       agctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATCGGCACCACCGCCGGGCTTAAAc                             >  2:3349944/1‑93 (MQ=255)
                                                        gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                            <  2:5805894/93‑1 (MQ=255)
                                                         ctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGa                           <  2:5407271/93‑1 (MQ=255)
                                                          tgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACATCTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGat                          <  1:1712794/93‑1 (MQ=255)
                                                          tgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGat                          >  2:3804020/1‑93 (MQ=255)
                                                           gaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata                         >  2:2184243/1‑93 (MQ=255)
                                                           gaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata                         >  1:5996774/1‑93 (MQ=255)
                                                           gaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata                         >  2:6439116/1‑93 (MQ=255)
                                                            aTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                        >  2:1401395/1‑93 (MQ=255)
                                                            aTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                        <  1:255484/93‑1 (MQ=255)
                                                             tGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGataca                       >  1:5985038/1‑93 (MQ=255)
                                                               cAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                     <  1:2463734/93‑1 (MQ=255)
                                                               cAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                     <  1:2898164/93‑1 (MQ=255)
                                                                 aaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                   <  2:2485998/93‑1 (MQ=255)
                                                                 aaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                   >  2:2731838/1‑93 (MQ=255)
                                                                   cGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaaa                 >  2:4312851/1‑93 (MQ=255)
                                                                   cGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaaa                 >  1:3241700/1‑93 (MQ=255)
                                                                    ggCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg                <  1:2677355/93‑1 (MQ=255)
                                                                     gcgcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgt               <  2:5623618/93‑1 (MQ=255)
                                                                      cgcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg              <  2:489807/93‑1 (MQ=255)
                                                                       gcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg             <  2:4600179/93‑1 (MQ=255)
                                                                        cTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt            <  1:5214625/93‑1 (MQ=255)
                                                                        cTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt            >  1:5743168/1‑93 (MQ=255)
                                                                         tAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggtg           >  1:1840737/1‑93 (MQ=255)
                                                                         tAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggtg           >  1:436299/1‑93 (MQ=255)
                                                                           aaCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCg         >  2:2928745/1‑93 (MQ=255)
                                                                            aCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa        >  1:6295841/1‑93 (MQ=255)
                                                                             cTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa       >  1:4837865/1‑93 (MQ=255)
                                                                             cTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa       >  1:1023043/1‑93 (MQ=255)
                                                                              tGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt      <  2:5058156/93‑1 (MQ=255)
                                                                              tGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt      <  1:3458908/93‑1 (MQ=255)
                                                                               gggCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg     >  2:5326530/1‑93 (MQ=255)
                                                                               gggCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg     >  1:5446055/1‑93 (MQ=255)
                                                                               gggCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg     <  1:4717613/93‑1 (MQ=255)
                                                                                 gCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa   <  2:2523724/93‑1 (MQ=255)
                                                                                 gCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa   <  1:3072720/93‑1 (MQ=255)
                                                                                  cGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAg  >  2:3700695/1‑93 (MQ=255)
                                                                                       ||                                                                                      
CTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAG  >  NC_000913/3972550‑3972724

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