Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 943,049 (A)5→4 coding (413/618 nt) dmsB → anaerobic dimethyl sulfoxide reductase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913943,0450A.85.7% 783.0 / 151.0 241coding (409/618 nt)dmsBanaerobic dimethyl sulfoxide reductase subunit B
Reads supporting (aligned to +/- strand):  ref base A (17/17);  new base . (107/97);  total (127/114)
Fisher's exact test for biased strand distribution p-value = 8.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 261 total aligned reads displayed.

ATATGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGC  >  NC_000913/942953‑943133
                                                                                            |                                                                                        
atatGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTa                                                                                          <  2:6036574/93‑1 (MQ=255)
atatGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTa                                                                                          <  1:6709575/93‑1 (MQ=255)
atatGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTa                                                                                          <  1:4721172/93‑1 (MQ=255)
atatGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTa                                                                                          >  2:7000366/1‑93 (MQ=255)
 tatGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTaa                                                                                         >  2:4637528/1‑93 (MQ=255)
  atGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTaaa                                                                                        <  2:2852257/93‑1 (MQ=255)
  atGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTaaa                                                                                        <  1:2459913/93‑1 (MQ=255)
    ggCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTaaaag                                                                                      <  1:6144259/93‑2 (MQ=255)
     gCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTaaaagc                                                                                     <  1:3092051/93‑3 (MQ=255)
      cTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGc                                                                                    >  1:1611772/1‑93 (MQ=255)
       ttGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGcc                                                                                   >  1:6485007/1‑93 (MQ=255)
        tGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGa                                                                                 >  1:1437674/1‑93 (MQ=255)
           ccGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATc                                                                               <  2:6550413/93‑1 (MQ=255)
           ccGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCt                                                                              >  1:2353297/1‑93 (MQ=255)
           ccGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCt                                                                              <  2:2171680/93‑1 (MQ=255)
             gTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCtgt                                                                            <  1:7457323/93‑1 (MQ=255)
              tATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCtgtg                                                                           <  1:4214856/93‑1 (MQ=255)
              tATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCtgtg                                                                           >  2:6237398/1‑93 (MQ=255)
                tGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGtt                                                                         >  2:3311638/1‑93 (MQ=255)
                tGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGtt                                                                         >  2:3935820/1‑93 (MQ=255)
                 ggCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTg                                                                        <  1:2114192/93‑1 (MQ=255)
                  gcgcACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTg                                                                        <  1:4590460/93‑1 (MQ=255)
                   cgcACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGaa                                                                      <  1:559642/93‑1 (MQ=255)
                      aCCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCt                                                                   >  1:440897/1‑93 (MQ=255)
                       ccGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCtt                                                                  <  2:5325717/93‑1 (MQ=255)
                       ccGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCtt                                                                  >  1:901996/1‑93 (MQ=255)
                        cGCAATATAACGAGACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCtt                                                                  <  2:2989837/93‑1 (MQ=255)
                        cGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTg                                                                 <  2:3209422/93‑1 (MQ=255)
                        cGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTg                                                                 >  1:5423272/1‑93 (MQ=255)
                         gCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGt                                                                >  2:2755754/1‑93 (MQ=255)
                         gCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGt                                                                <  2:6917489/93‑1 (MQ=255)
                           aaTATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTcc                                                              >  2:2499856/1‑93 (MQ=255)
                            atatAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCg                                                             <  2:6088229/93‑1 (MQ=255)
                               taACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTg                                                          >  1:3628308/1‑93 (MQ=255)
                                 acgaaacgaaAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGc                                                         >  2:7104710/1‑93 (MQ=255)
                                   gaaacgaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCgg                                                       <  1:7435687/93‑1 (MQ=255)
                                   gaaacgaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCggg                                                      <  2:1414449/93‑1 (MQ=255)
                                   gaaacgaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCggg                                                      >  1:3841188/1‑93 (MQ=255)
                                   gaaacgaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCggg                                                      <  1:3343605/93‑1 (MQ=255)
                                   gaaacgaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCggg                                                      >  2:7454677/1‑93 (MQ=255)
                                    aaacgaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGgc                                                     >  1:1639630/1‑93 (MQ=255)
                                      acgaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGgcgc                                                   >  1:5654697/1‑93 (MQ=255)
                                      acgaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGgcgc                                                   <  1:769261/93‑1 (MQ=255)
                                        gaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCt                                                  >  1:4499634/1‑93 (MQ=255)
                                        gaaaGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTg                                                 <  1:5898474/93‑1 (MQ=255)
                                            ggCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGAtt                                             <  2:4198570/93‑1 (MQ=255)
                                            ggCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGAtt                                             <  2:234643/93‑1 (MQ=255)
                                             gCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGAttt                                            >  1:2976366/1‑93 (MQ=255)
                                             gCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGAttt                                            <  2:3645623/93‑1 (MQ=255)
                                              ccATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTc                                           >  1:734544/1‑93 (MQ=255)
                                              ccATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTc                                           <  2:6312946/93‑1 (MQ=255)
                                              ccATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTc                                           <  2:1212048/93‑1 (MQ=255)
                                              ccATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTc                                           <  2:1113636/93‑1 (MQ=255)
                                               cATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCg                                          <  1:259379/93‑1 (MQ=255)
                                               cATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCg                                          >  2:4390939/1‑93 (MQ=255)
                                               cATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCg                                          >  1:6780134/1‑93 (MQ=255)
                                                 tatGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc                                        >  2:5111792/1‑93 (MQ=255)
                                                 tatGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc                                        <  2:5433176/93‑1 (MQ=255)
                                                   tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGccc                                      >  2:2034165/1‑93 (MQ=255)
                                                   tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGccc                                      >  1:5073641/1‑93 (MQ=255)
                                                    gACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCt                                     >  2:5450847/1‑93 (MQ=255)
                                                       cAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                                  <  1:1866757/93‑1 (MQ=255)
                                                       cAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                                  <  2:3526579/93‑1 (MQ=255)
                                                       cAAATGCGATGGTTGTTATAACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                                  >  1:7126033/1‑93 (MQ=255)
                                                         aaTGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATcga                                <  2:2707830/93‑1 (MQ=255)
                                                         aaTGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATcga                                >  2:1936440/1‑93 (MQ=255)
                                                          aTGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATcgac                               >  2:1775014/1‑93 (MQ=255)
                                                            gCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATcgacga                             >  1:4478421/1‑93 (MQ=255)
                                                              gATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAg                            >  2:3769605/1‑93 (MQ=255)
                                                              gATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGc                           <  1:5787000/93‑1 (MQ=255)
                                                               aTGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCt                          >  1:4014747/1‑93 (MQ=255)
                                                               aTGGTTGTTATGACCGTGTTGCGGAGGGGAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCt                          >  1:1148685/1‑93 (MQ=255)
                                                                tGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTg                         >  2:3708543/1‑93 (MQ=255)
                                                                tGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTg                         >  1:1775759/1‑93 (MQ=255)
                                                                tGGTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTg                         >  2:4796394/1‑93 (MQ=255)
                                                                 ggTTGTTATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGc                        <  2:2608487/93‑1 (MQ=255)
                                                                  gttgttATGACCGTGTTGCGGAGGGTAAAAGCGGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg                       >  2:1051753/1‑93 (MQ=255)
                                                                  gttgttATGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg                       >  1:4901709/1‑93 (MQ=255)
                                                                     gttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTa                     >  2:7342204/1‑93 (MQ=255)
                                                                        aTGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATCTCGGCCCTATCGACGAGCTGCGTaaaaa                 <  1:5273414/93‑1 (MQ=255)
                                                                         tGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa                <  2:5599183/93‑1 (MQ=255)
                                                                         tGACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa                >  2:1643150/1‑93 (MQ=255)
                                                                          gACCGTGTTGCGGATGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAAc               >  2:6938341/1‑93 (MQ=255)
                                                                          gACCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAAc               <  1:2650289/93‑1 (MQ=255)
                                                                           aCCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACa              >  2:229823/1‑93 (MQ=255)
                                                                           aCCGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACa              <  1:1039792/93‑1 (MQ=255)
                                                                            ccGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAt             <  1:7193358/93‑1 (MQ=255)
                                                                            ccGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAt             <  1:6314263/93‑1 (MQ=255)
                                                                            ccGTGTTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAt             >  1:5135569/1‑93 (MQ=255)
                                                                               tgtTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGc          <  2:5703927/93‑1 (MQ=255)
                                                                               tgtTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGc          <  2:4540097/93‑1 (MQ=255)
                                                                               tgtTGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGc          <  1:3592474/93‑1 (MQ=255)
                                                                                gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGc          <  2:743755/93‑1 (MQ=255)
                                                                                  tGCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGAt       <  2:5507266/93‑1 (MQ=255)
                                                                                   gCGGAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATc      <  1:2223303/93‑1 (MQ=255)
                                                                                     ggAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTg    <  2:896390/93‑1 (MQ=255)
                                                                                      gAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTgg   >  1:1726546/1‑93 (MQ=255)
                                                                                      gAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTgg   <  1:5173517/93‑1 (MQ=255)
                                                                                      gAGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTgg   <  2:5584398/93‑1 (MQ=255)
                                                                                       aGGGTAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTggc  <  2:3595807/93‑1 (MQ=255)
                                                                                            |                                                                                        
ATATGGCTTGCCCGTATGGCGCACCGCAATATAACGAAACGAAAGGCCATATGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGC  >  NC_000913/942953‑943133

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