Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A12 F56 I0 R1 312 170.3 8276235 93.7% 7754832 99.8

Breseq alignment

BRESEQ :: Evidence
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)
       aTGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATg                                                                                             >  2:656817/1‑100 (MQ=255)
        tGAAACCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATgg                                                                                            <  1:809668/100‑1 (MQ=255)
           aaCCTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGc                                                                                          >  1:2496986/1‑100 (MQ=255)
             ccTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTg                                                                                       >  2:1716463/1‑100 (MQ=255)
              cTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGa                                                                                      <  1:77705/100‑1 (MQ=255)
              cTGGCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGa                                                                                      <  2:3827176/100‑1 (MQ=255)
                ggCAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATg                                                                                    >  2:3753427/1‑100 (MQ=255)
                  cAGAGCCTGAAAGACGCCTTCAGGCTTAGCGTACTGATGAACCCCATGTAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATgg                                                                                   <  2:1607871/100‑1 (MQ=255)
                  cAGAGCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATgg                                                                                   >  1:2064611/1‑100 (MQ=255)
                    gagCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATc                                                                                <  1:3704030/100‑1 (MQ=255)
                     agCCTGAAAGACACCTTCAGCCTTAGCGTACCGATGCACCCCATCGAGGCCATCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTTATGGAt                                                                                 <  2:3444527/98‑1 (MQ=255)
                      gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAg                                                                              <  1:1010853/100‑1 (MQ=255)
                      gCCTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAg                                                                              >  2:555039/1‑100 (MQ=255)
                        cTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAACAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGaa                                                                            >  2:3416007/1‑100 (MQ=255)
                        cTGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTCAAGCCAGCAAAAGCGACGGCGAACTTGATACC‑GGGCGATGGACGTGATGGATCAGaa                                                                            >  2:2601727/1‑100 (MQ=255)
                         tGAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGaaa                                                                           >  2:3671119/1‑100 (MQ=255)
                          gAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCACAGACGACGGCGAAATGGGTACC‑TGGCGCTGGGCGTGATGGATCaaaa                                                                            >  2:3272072/1‑98 (MQ=255)
                          gAAAGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAAt                                                                          >  1:1131421/1‑100 (MQ=255)
                            aaGACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTc                                                                        <  2:1269219/100‑1 (MQ=255)
                              gACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCg                                                                       <  2:4024685/100‑1 (MQ=255)
                              gACACCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCgg                                                                      >  1:1597607/1‑100 (MQ=255)
                               acacCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCgg                                                                      <  2:2184776/100‑1 (MQ=255)
                                cacCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTa                                                                    >  2:767932/1‑100 (MQ=255)
                                 acCTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAg                                                                   <  1:1379614/100‑1 (MQ=255)
                                  ccTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAgc                                                                  >  1:1873845/1‑100 (MQ=255)
                                  ccTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAgc                                                                  <  1:493348/100‑1 (MQ=255)
                                   cTTCAGCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGACGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAgcg                                                                 <  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)
                                        gCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagca                                                            >  1:1895315/1‑100 (MQ=255)
                                        gCCTTAGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCgcagca                                                            <  2:2695804/100‑1 (MQ=255)
                                             aGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCagcagcagc                                                       <  2:1032984/100‑1 (MQ=255)
                                             aGCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCagcagcagc                                                       <  2:1037089/100‑1 (MQ=255)
                                              gCGTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCagcagcagc                                                       >  1:396258/1‑100 (MQ=255)
                                                gTACTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTa                                                    <  1:3189369/100‑1 (MQ=255)
                                                  aCTGATGAATCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGTATCATAAATTTGGTAGCGCACCAGCAGCTGACa                                                  <  2:1017340/100‑1 (MQ=255)
                                                  aCTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTAc                                                   <  1:146504/100‑1 (MQ=255)
                                                   cTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAg                                                 <  1:3166764/100‑1 (MQ=255)
                                                   cTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAg                                                 <  2:2431797/100‑1 (MQ=255)
                                                   cTGATGAACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAg                                                 <  2:3326851/100‑1 (MQ=255)
                                                     gatgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGc                                                >  1:634987/1‑100 (MQ=255)
                                                     gatgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGc                                                <  2:1272304/100‑1 (MQ=255)
                                                     gatgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGc                                                <  2:4082351/99‑1 (MQ=255)
                                                     gatgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCt                                               >  2:1381029/1‑100 (MQ=255)
                                                     gatgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCt                                               >  1:3591250/1‑100 (MQ=255)
                                                      atgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTa                                              <  1:1437256/100‑1 (MQ=255)
                                                       tgaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTa                                              >  2:2937371/1‑100 (MQ=255)
                                                        gaACCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTAcc                                            <  2:2342436/100‑1 (MQ=255)
                                                          aCCCCATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTg                                          >  1:2009688/1‑100 (MQ=255)
                                                            cccATTGAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGAt                                        >  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)
                                                                  gAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGt                                  >  1:2328399/1‑100 (MQ=255)
                                                                  gAAGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGt                                  >  2:1433664/1‑100 (MQ=255)
                                                                   aaGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGt                                  >  2:1185984/1‑100 (MQ=255)
                                                                   aaGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTc                                 >  1:1951288/1‑100 (MQ=255)
                                                                    aGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTc                                 >  2:1928092/1‑100 (MQ=255)
                                                                    aGCCAGCAAAGGCGACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTc                                 >  1:3864603/1‑100 (MQ=255)
                                                                    aGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTcc                                >  2:2146816/1‑100 (MQ=255)
                                                                    aGCCAGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTcc                                <  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)
                                                                        aGCAAAGGCGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCctgct                            <  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)
                                                                                cGACGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTAc                    <  1:3909052/100‑1 (MQ=255)
                                                                                 gACGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTAc                    <  1:1659486/100‑1 (MQ=255)
                                                                                  aCGGCGAAATGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTAcc                   >  2:3826614/1‑100 (MQ=255)
                                                                                  aCGGCGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCa                  <  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)
                                                                                      cGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGtn              <  2:2214258/100‑2 (MQ=255)
                                                                                      cGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTa              <  1:1264701/100‑1 (MQ=255)
                                                                                      cGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTa              <  2:2180258/100‑1 (MQ=255)
                                                                                      cGAAATGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCAGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTa              <  2:2904436/100‑1 (MQ=255)
                                                                                          aTGGGTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTAtct           >  2:2473947/1‑100 (MQ=255)
                                                                                           tGGGTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCg         >  1:2259877/1‑100 (MQ=255)
                                                                                            gggTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCg         >  1:2586101/1‑100 (MQ=255)
                                                                                            gggTACCGGGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCg         >  1:3250000/1‑100 (MQ=255)
                                                                                            gggTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCgg        <  1:3479667/100‑1 (MQ=255)
                                                                                             ggTACCGGGGCGATGGGCGTGATGGACCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCgg        >  1:2523705/1‑100 (MQ=255)
                                                                                             ggTACC‑GGGCGATGGGCGTGATGGATCAGAAATTCGGTAGCGCAGCAGCAGCTTACAGCTACCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGt       <  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: ‑

GATK/CNVnator alignment

N/A