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% 1678.7 / NA 489R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (229/260);  total (229/260)
*NC_0009133,972,6380GC99.8% 1679.9 / ‑5.1 489R144A (CGT→GCT) wecDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (228/260);  minor base T (1/0);  total (229/260)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 500 total aligned reads displayed.

CTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAG  >  NC_000913/3972550‑3972724
                                                                                       ||                                                                                      
  cAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCgg                                                                                  <  1:4171417/93‑1 (MQ=255)
   aaTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGt                                                                                 >  2:3276969/1‑93 (MQ=255)
    aTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTa                                                                                >  2:6820935/1‑93 (MQ=255)
    aTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTa                                                                                >  1:6149841/1‑93 (MQ=255)
    aTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTa                                                                                <  2:8140378/93‑1 (MQ=255)
        gACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTaaaaa                                                                            >  2:3060562/1‑92 (MQ=255)
        gACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaac                                                                            <  1:6775648/93‑1 (MQ=255)
         aCAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaaca                                                                           >  1:8539054/1‑93 (MQ=255)
         aCAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaaca                                                                           >  1:4961875/1‑93 (MQ=255)
          cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaa                                                                          <  1:151821/93‑1 (MQ=255)
          cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGGGCTGGGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacca                                                                          >  2:4441627/1‑91 (MQ=255)
           aGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAaacaac                                                                         <  2:5829793/93‑1 (MQ=255)
            gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACt                                                                        >  2:1352614/1‑93 (MQ=255)
            gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACt                                                                        <  1:1332954/93‑1 (MQ=255)
              tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACttt                                                                      <  1:4629500/93‑1 (MQ=255)
              tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAAAttt                                                                      >  1:5891574/1‑93 (MQ=255)
              agcgcGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGGAAAACAACttt                                                                      <  2:2609268/92‑1 (MQ=255)
                cgcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGc                                                                    <  1:100120/93‑1 (MQ=255)
                 gcgAATTGGCCTGCTGGCTGGACGCGGGGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGAGAGAATGCTCCCGGTAAAACAACTTTGCg                                                                   >  2:7948809/1‑93 (MQ=255)
                  cgAATTGGCCTGCTGGCTGGACCCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAACCAACTTTGCgg                                                                  <  1:4063005/93‑1 (MQ=255)
                   gAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCggg                                                                 <  1:5666471/93‑1 (MQ=255)
                   gAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCggg                                                                 >  2:911405/1‑93 (MQ=255)
                   gAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCggg                                                                 <  1:1671438/93‑1 (MQ=255)
                    aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTTCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGt                                                                <  1:6644412/93‑1 (MQ=255)
                     aTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                                               <  2:3487856/93‑1 (MQ=255)
                     aTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                                               >  1:8086233/1‑93 (MQ=255)
                     aTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                                               <  2:115984/93‑1 (MQ=255)
                     aTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                                               >  1:7654069/1‑93 (MQ=255)
                     aTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTg                                                               >  2:2254551/1‑93 (MQ=255)
                      ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAAAAACTTTGCGGGTgg                                                              >  1:7761226/1‑93 (MQ=255)
                        ggCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCg                                                            >  2:4477974/1‑93 (MQ=255)
                        ggCCTGCTGGCTGGACGCGGTGCAGGCGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCg                                                            >  1:2892778/1‑93 (MQ=255)
                         gCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGa                                                           <  1:3535087/93‑1 (MQ=255)
                          ccTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGAc                                                          <  2:39432/93‑1 (MQ=255)
                           ctgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGAcc                                                         <  1:640556/93‑1 (MQ=255)
                           ctgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGAcc                                                         <  2:4288249/93‑1 (MQ=255)
                               tggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCaaa                                                     >  2:2910411/1‑93 (MQ=255)
                               tggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCaaa                                                     >  1:3470009/1‑93 (MQ=255)
                                 gctggACGCGGTGCAGGTGCTGATCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                   <  2:6276015/93‑1 (MQ=255)
                                 gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                   <  1:368761/93‑1 (MQ=255)
                                 gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                   <  1:7556284/93‑1 (MQ=255)
                                  ctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATgg                                                  >  2:726659/1‑93 (MQ=255)
                                   tggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATggg                                                 >  1:4412301/1‑93 (MQ=255)
                                    ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                                >  2:8696677/1‑93 (MQ=255)
                                     gACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCa                                               >  2:7591594/1‑93 (MQ=255)
                                       cgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAac                                             <  1:3365422/93‑1 (MQ=255)
                                       cgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAac                                             <  2:7704161/93‑1 (MQ=255)
                                       cgcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAac                                             <  2:8546029/93‑1 (MQ=255)
                                        gcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAaca                                            >  2:3702333/1‑93 (MQ=255)
                                         cgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAacac                                           <  1:8602633/93‑1 (MQ=255)
                                          ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                          >  1:5489745/1‑93 (MQ=255)
                                          ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                          <  2:1501575/93‑1 (MQ=255)
                                           gTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                                         >  2:2404933/1‑93 (MQ=255)
                                           gTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                                         >  1:505874/1‑93 (MQ=255)
                                            tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        <  2:5474884/93‑1 (MQ=255)
                                            tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        <  2:7183941/93‑1 (MQ=255)
                                            tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        >  1:2283320/1‑93 (MQ=255)
                                             gCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgcc                                       <  2:6183879/93‑1 (MQ=255)
                                                 gTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCgcgc                                   >  1:2132479/1‑93 (MQ=255)
                                                  tGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                                  >  1:7935225/1‑93 (MQ=255)
                                                    ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                <  2:6670233/93‑1 (MQ=255)
                                                    ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                >  1:8758037/1‑93 (MQ=255)
                                                     tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCACATGGGCAACACCGCCGCGCTTaa                               <  1:1957278/93‑1 (MQ=255)
                                                     tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa                               >  2:5658903/1‑93 (MQ=255)
                                                     tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa                               >  2:5382131/1‑93 (MQ=255)
                                                      gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                              >  2:4155294/1‑93 (MQ=255)
                                                      gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                              <  1:2741374/93‑1 (MQ=255)
                                                      gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                              <  2:5143720/93‑1 (MQ=255)
                                                       agctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc                             >  2:4885447/1‑93 (MQ=255)
                                                        gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                            <  1:2567985/93‑1 (MQ=255)
                                                          tgaTGCAAACGGCGCTAAACTGGGCGAATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGat                          <  2:6009793/93‑1 (MQ=255)
                                                           gaTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata                         <  1:4762367/93‑1 (MQ=255)
                                                            aTGCAAACGTCGCTACACTGGGCGTATGCTCGCGGTAAAACAACGTTGCGGGCGGCGACCAAAATGGGCCACACCGCCGCGCTTAAACGatac                        >  2:6978360/1‑93 (MQ=255)
                                                            aTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                        <  1:650309/93‑1 (MQ=255)
                                                            aTGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                        >  1:320392/1‑93 (MQ=255)
                                                             tGCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGataca                       <  2:1686430/93‑1 (MQ=255)
                                                              gCAAACGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCAACCCAAATGGGAAACACTCCCGCGCTTAAACGatacaa                      >  2:1763081/1‑92 (MQ=255)
                                                                aaaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                    >  2:4719987/1‑93 (MQ=255)
                                                                 aaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                   >  1:2413553/1‑93 (MQ=255)
                                                                 aaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                   >  1:1811618/1‑93 (MQ=255)
                                                                 aaCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                   >  2:1929149/1‑93 (MQ=255)
                                                                  aCGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaa                  >  2:210699/1‑93 (MQ=255)
                                                                   cGGCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaaa                 <  1:7037508/93‑1 (MQ=255)
                                                                    ggCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg                >  2:5334763/1‑93 (MQ=255)
                                                                    ggCGCTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg                <  1:8328878/93‑1 (MQ=255)
                                                                     gcgcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAga               >  1:6935242/1‑92 (MQ=255)
                                                                      ggcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg              <  1:3895462/92‑1 (MQ=255)
                                                                      cgcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg              <  2:7878866/93‑1 (MQ=255)
                                                                       gcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg             <  1:4002592/93‑1 (MQ=255)
                                                                       gcTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg             <  2:1077909/93‑1 (MQ=255)
                                                                       gcTAAACGGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg             <  1:8632865/93‑1 (MQ=255)
                                                                        cTAAACTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt            <  1:3058779/93‑1 (MQ=255)
                                                                          aaaCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc          >  1:7367721/1‑93 (MQ=255)
                                                                           aaCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCg         >  1:7201671/1‑93 (MQ=255)
                                                                            aCTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa        >  2:7413521/1‑93 (MQ=255)
                                                                             cTGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa       <  2:3541129/93‑1 (MQ=255)
                                                                              tGGGCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt      >  2:8374206/1‑93 (MQ=255)
                                                                               gggCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg     <  1:5261219/93‑1 (MQ=255)
                                                                                ggCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGt    <  2:379799/93‑1 (MQ=255)
                                                                                ggCGTATGCTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGt    <  1:1235650/93‑1 (MQ=255)
                                                                                       ||                                                                                      
CTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAG  >  NC_000913/3972550‑3972724

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