Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,542,354 Δ1 bp 100% coding (1948/2322 nt) feoB → ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,542,3510G.54.7% 10.5 / inf 172coding (1945/2322 nt)feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (40/37);  new base . (39/55);  total (79/93)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

Only 100 of 182 total aligned reads displayed.

GCGATAGATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGG  >  NC_000913/3542252‑3542450
                                                                                                   |                                                                                                   
 cGATAGATGAAACATGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCgg                                                                                                    <  1:2677529/100‑1 (MQ=255)
  gcttgacgtaaCCTGGCCGAGCCTGAAAGACACCTGCCGCCTTAGCGGACTGATGAACCCCATTGACGCCTGCAAAGGCGCCGGCGAAATGGGTACCgg                                                                                                    <  2:1463493/90‑1 (MQ=255)
  gATAGATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCggg                                                                                                   >  1:2600796/1‑100 (MQ=255)
   aTAGATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCgggc                                                                                                  <  1:3588668/100‑2 (MQ=255)
    tAGATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGc                                                                                                 >  1:1942463/1‑100 (MQ=255)
     aGATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCgggcga                                                                                                <  1:2750391/100‑4 (MQ=255)
      gATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCgggcga                                                                                                >  2:1860464/1‑96 (MQ=255)
      gATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGa                                                                                               >  1:3337309/1‑100 (MQ=255)
       aTGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGAt                                                                                              <  2:2131094/100‑1 (MQ=255)
       aTGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGAt                                                                                              >  2:363129/1‑100 (MQ=255)
       aTGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATg                                                                                             >  2:656817/1‑100 (MQ=255)
        tGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATgg                                                                                            <  1:809668/100‑1 (MQ=255)
           aaCCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGc                                                                                          >  1:2496986/1‑100 (MQ=255)
             ccTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTg                                                                                       >  2:1716463/1‑100 (MQ=255)
              cTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGa                                                                                      <  1:77705/100‑1 (MQ=255)
              cTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGa                                                                                      <  2:3827176/100‑1 (MQ=255)
                ggCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATg                                                                                    >  2:3753427/1‑100 (MQ=255)
                  cAGAGCCTGAAAGACGCCTTCAGGCTTAGCGTACTGATGAACCCCATGTAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATgg                                                                                   <  2:1607871/100‑1 (MQ=255)
                  cAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATgg                                                                                   >  1:2064611/1‑100 (MQ=255)
                    gagCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATc                                                                                <  1:3704030/100‑1 (MQ=255)
                     agCCTGAAAGACACCTTCAGCCTTAGCGTACCGATGCACCCCATCGAGGCCATCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTTATGGAt                                                                                 <  2:3444527/98‑1 (MQ=255)
                      gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAg                                                                              <  1:1010853/100‑1 (MQ=255)
                      gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAg                                                                              >  2:555039/1‑100 (MQ=255)
                        cTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAACAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGaa                                                                            >  2:3416007/1‑100 (MQ=255)
                        cTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTCAAGCCAGCAAAAGCGACGGCGAACTTGATACCGGGCGATGGACGTGATGGATCAGaa                                                                            >  2:2601727/1‑100 (MQ=255)
                         tGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGaaa                                                                           >  2:3671119/1‑100 (MQ=255)
                          gAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCACAGACGACGGCGAAATGGGTACCTGGCGCTGGGCGTGATGGATCaaaa                                                                            >  2:3272072/1‑98 (MQ=255)
                          gAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAAt                                                                          >  1:1131421/1‑100 (MQ=255)
                            aaGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTc                                                                        <  2:1269219/100‑1 (MQ=255)
                              gACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCg                                                                       <  2:4024685/100‑1 (MQ=255)
                              gACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCgg                                                                      >  1:1597607/1‑100 (MQ=255)
                               acacCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCgg                                                                      <  2:2184776/100‑1 (MQ=255)
                                cacCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTa                                                                    >  2:767932/1‑100 (MQ=255)
                                 acCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAg                                                                   <  1:1379614/100‑1 (MQ=255)
                                  ccTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAgc                                                                  >  1:1873845/1‑100 (MQ=255)
                                  ccTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAgc                                                                  <  1:493348/100‑1 (MQ=255)
                                   cTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGACGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAgcg                                                                 <  1:2604242/100‑1 (MQ=255)
                                    ttCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAgcg                                                                 >  2:1801160/1‑100 (MQ=255)
                                       aGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcag                                                              <  1:987935/100‑1 (MQ=255)
                                        gCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagc                                                             <  1:301693/100‑1 (MQ=255)
                                        gCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagca                                                            >  1:1895315/1‑100 (MQ=255)
                                        gCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagca                                                            <  2:2695804/100‑1 (MQ=255)
                                             aGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCagcagcagc                                                       <  2:1032984/100‑1 (MQ=255)
                                             aGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCagcagcagc                                                       <  2:1037089/100‑1 (MQ=255)
                                              gCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCagcagcagc                                                       >  1:396258/1‑100 (MQ=255)
                                                gTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTa                                                    <  1:3189369/100‑1 (MQ=255)
                                                  aCTGATGAATCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGTATCATAAATTTGGTAGCGCACCAGCAGCTGACa                                                  <  2:1017340/100‑1 (MQ=255)
                                                  aCTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTAc                                                   <  1:146504/100‑1 (MQ=255)
                                                   cTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAg                                                 <  1:3166764/100‑1 (MQ=255)
                                                   cTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAg                                                 <  2:2431797/100‑1 (MQ=255)
                                                   cTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAg                                                 <  2:3326851/100‑1 (MQ=255)
                                                     gatgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGc                                                >  1:634987/1‑100 (MQ=255)
                                                     gatgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGc                                                <  2:1272304/100‑1 (MQ=255)
                                                     gatgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGc                                                <  2:4082351/99‑1 (MQ=255)
                                                     gatgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCt                                               >  2:1381029/1‑100 (MQ=255)
                                                     gatgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCt                                               >  1:3591250/1‑100 (MQ=255)
                                                      atgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTa                                              <  1:1437256/100‑1 (MQ=255)
                                                       tgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTa                                              >  2:2937371/1‑100 (MQ=255)
                                                        gaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTAcc                                            <  2:2342436/100‑1 (MQ=255)
                                                          aCCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTg                                          >  1:2009688/1‑100 (MQ=255)
                                                            cccATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGAt                                        >  2:2242759/1‑100 (MQ=255)
                                                             ccATTGAAGCCAGCAAAGGCGACGGCGAAATGGATACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGAt                                        >  2:2473923/1‑100 (MQ=255)
                                                                  gAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCg                                   <  2:943037/100‑1 (MQ=255)
                                                                  gAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCg                                   >  2:632436/1‑100 (MQ=255)
                                                                  gAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGt                                  >  1:2328399/1‑100 (MQ=255)
                                                                  gAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGt                                  >  2:1433664/1‑100 (MQ=255)
                                                                   aaGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGt                                  >  2:1185984/1‑100 (MQ=255)
                                                                   aaGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTc                                 >  1:1951288/1‑100 (MQ=255)
                                                                    aGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTc                                 >  2:1928092/1‑100 (MQ=255)
                                                                    aGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTc                                 >  1:3864603/1‑100 (MQ=255)
                                                                    aGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTcc                                >  2:2146816/1‑100 (MQ=255)
                                                                    aGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTcc                                <  1:2204544/100‑1 (MQ=255)
                                                                     gCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTcc                                <  1:2911792/100‑1 (MQ=255)
                                                                     gCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTcc                                >  2:3035323/1‑100 (MQ=255)
                                                                        aggaaaGGCGACGGGGTAATGGGTACCGGGTGGATGGGTGGGATAGATCGGAAATTCGGTAGCGCAGCAGCAGCTTACAGCGACCTGATTTTCGt                                  <  2:824152/92‑1 (MQ=255)
                                                                        aGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCctgc                             <  1:827563/100‑1 (MQ=255)
                                                                        aGCAAAGGCGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCctgct                            <  1:36327/100‑1 (MQ=255)
                                                                         gCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCctgct                            <  2:2064598/100‑1 (MQ=255)
                                                                         gCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCctgct                            >  1:4044585/1‑100 (MQ=255)
                                                                         gCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCctgct                            >  2:1214807/1‑100 (MQ=255)
                                                                             aGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCtgtat                        >  1:740951/1‑100 (MQ=255)
                                                                               gCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCtgtatgt                      <  2:2534394/100‑1 (MQ=255)
                                                                                cGACGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTAc                    <  1:3909052/100‑1 (MQ=255)
                                                                                 gACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTAc                    <  1:1659486/100‑1 (MQ=255)
                                                                                  aCGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTAcc                   >  2:3826614/1‑100 (MQ=255)
                                                                                  aCGGCGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCa                  <  2:491621/100‑1 (MQ=255)
                                                                                   cGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCa                  >  2:883560/1‑100 (MQ=255)
                                                                                   cGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCa                  >  1:1752368/1‑100 (MQ=255)
                                                                                    ggCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCAt                 <  2:2781708/100‑1 (MQ=255)
                                                                                      cGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGtn              <  2:2214258/100‑2 (MQ=255)
                                                                                      cGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTa              <  1:1264701/100‑1 (MQ=255)
                                                                                      cGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTa              <  2:2180258/100‑1 (MQ=255)
                                                                                      cGAAATGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCAGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTa              <  2:2904436/100‑1 (MQ=255)
                                                                                          aTGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTAtct           >  2:2473947/1‑100 (MQ=255)
                                                                                           tGGGTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCg         >  1:2259877/1‑100 (MQ=255)
                                                                                            gggTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCg         >  1:2586101/1‑100 (MQ=255)
                                                                                            gggTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCg         >  1:3250000/1‑100 (MQ=255)
                                                                                            gggTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCgg        <  1:3479667/100‑1 (MQ=255)
                                                                                             ggTACCGGGGCGATGGGCGTGATGGACCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCgg        >  1:2523705/1‑100 (MQ=255)
                                                                                             ggTACCGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGt       <  1:1388509/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCGATAGATGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGG  >  NC_000913/3542252‑3542450

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