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
A1 F1 I2 R1 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,421 +T coding (504/516 nt) b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4211.T100.0% 359.4 / NA 100L168F (TTG→TTTb4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (50/50);  total (50/50)

Only 100 of 101 total aligned reads displayed.

TCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCT  >  NC_000913/4008346‑4008503
                                                                                     |                                                                                         
tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                                                                   >  2:3254248/1‑93 (MQ=255)
tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                                                                   <  1:2592961/93‑1 (MQ=255)
tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                                                                   <  1:4432993/93‑1 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                                                                                 <  1:3632859/93‑3 (MQ=255)
  tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                                                                                 <  1:3793677/93‑3 (MQ=255)
   ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                                                                                >  1:1240529/1‑90 (MQ=255)
    gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                                                                               <  1:1051163/93‑5 (MQ=255)
    gAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                                                                               >  1:3702178/1‑89 (MQ=255)
          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                                                                         <  2:760103/93‑1 (MQ=255)
          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                                                                         <  2:229684/93‑1 (MQ=255)
          aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                                                                         >  2:1077404/1‑93 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                                                                        <  1:1894017/93‑1 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                                                                        <  1:1566229/93‑1 (MQ=255)
           tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                                                                        <  1:3374702/93‑1 (MQ=255)
            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                                                                       <  1:401856/93‑1 (MQ=255)
            ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCCGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                                                                       <  2:3408020/93‑1 (MQ=255)
             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                                                                      <  2:1467147/93‑1 (MQ=255)
             gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                                                                      <  1:3117785/93‑1 (MQ=255)
               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                                                                    <  1:4425676/93‑1 (MQ=255)
               gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                                                                    <  1:1119358/93‑1 (MQ=255)
                 gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                                                                  >  1:4457103/1‑93 (MQ=255)
                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga                                                                 >  1:4576177/1‑93 (MQ=255)
                  cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga                                                                 <  2:3526130/93‑1 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                                <  1:2958163/93‑1 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                                <  2:4322943/93‑1 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                                <  1:2669860/93‑1 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                                >  2:516542/1‑93 (MQ=255)
                   ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac                                                                >  2:4523682/1‑93 (MQ=255)
                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgggg                                                               >  1:2074904/1‑90 (MQ=255)
                    gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg                                                               >  2:2058391/1‑93 (MQ=255)
                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc                                                              >  2:4150932/1‑93 (MQ=255)
                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc                                                              <  2:1382101/93‑1 (MQ=255)
                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc                                                              <  1:1265185/93‑1 (MQ=255)
                     cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc                                                              >  1:2192640/1‑93 (MQ=255)
                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc                                                             <  2:1470272/93‑1 (MQ=255)
                      gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc                                                             >  2:1671891/1‑93 (MQ=255)
                      gTTTCAAAAGGCCATGCGCAGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc                                                             >  2:4571105/1‑93 (MQ=255)
                          canaaGNNCATGCGCTGGAAGNNGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                                         <  1:4559497/90‑1 (MQ=255)
                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                                         <  1:2234580/91‑1 (MQ=255)
                          caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                                         <  2:3847674/91‑1 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                                        >  2:1013108/2‑93 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                                        >  1:4466518/2‑93 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                                        >  1:2030356/2‑93 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                                        >  2:2669759/2‑93 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                                        >  1:2475525/2‑93 (MQ=255)
                           aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                                        <  2:4362651/92‑1 (MQ=255)
                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                                       >  2:470735/1‑93 (MQ=255)
                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                                       >  1:1050044/1‑93 (MQ=255)
                            aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                                       >  1:4590471/1‑93 (MQ=255)
                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                      >  1:3864513/1‑93 (MQ=255)
                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                      <  2:1078745/93‑1 (MQ=255)
                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                      <  1:4450069/93‑1 (MQ=255)
                             aaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                                      <  1:2907437/93‑1 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                     <  2:2559408/93‑1 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                     >  1:4654945/1‑93 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                     <  1:4631955/93‑1 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                     >  2:1507239/1‑93 (MQ=255)
                              aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                                     <  1:3115899/93‑1 (MQ=255)
                               ggCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg                                                    <  2:3964472/93‑1 (MQ=255)
                                                              ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGc                   <  2:2551572/90‑1 (MQ=255)
                                                              ctgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGc                   >  2:2827746/4‑93 (MQ=255)
                                                               tgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCa                  >  1:1030894/3‑93 (MQ=255)
                                                               tgggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCa                  <  1:3975408/91‑1 (MQ=255)
                                                                gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAg                 <  2:4494890/92‑1 (MQ=255)
                                                                gggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAg                 <  2:483481/92‑1 (MQ=255)
                                                                gggCTACAGCCTGAAGGATTGTATTGCGTTTGGCGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAg                 >  2:4215796/2‑93 (MQ=255)
                                                                 ggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt                <  1:2506546/93‑1 (MQ=255)
                                                                 ggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt                >  2:3217413/1‑93 (MQ=255)
                                                                 ggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGt                >  2:2396445/1‑93 (MQ=255)
                                                                   cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc              <  1:3340585/93‑2 (MQ=255)
                                                                   cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc              >  1:2035757/1‑92 (MQ=255)
                                                                   cTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgc              <  2:2953294/93‑2 (MQ=255)
                                                                    tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg             >  2:2407127/1‑91 (MQ=255)
                                                                    tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg             >  1:689379/1‑91 (MQ=255)
                                                                    tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg             >  2:2696585/1‑91 (MQ=255)
                                                                    tACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTgcg             >  1:1686572/1‑91 (MQ=255)
                                                                        gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcacc         >  2:4246988/3‑93 (MQ=255)
                                                                        gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcacc         >  2:4074519/3‑93 (MQ=255)
                                                                        gccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcacc         >  1:4332377/3‑93 (MQ=255)
                                                                         ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcaccg        >  1:750969/2‑92 (MQ=255)
                                                                          cTGAAGGATTGTATTGCGTTTGGTGACGGGGTGAACGGCGCCGTAATGCTGTCGAAGGCGGGG‑‑AAAGCCTGCATTATGGGCAGGGGGTACCAg       >  2:1076668/1‑93 (MQ=255)
                                                                          cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAg       <  1:3892728/93‑1 (MQ=255)
                                                                          cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAg       >  1:1192872/1‑93 (MQ=255)
                                                                          cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAg       <  1:960598/93‑1 (MQ=255)
                                                                           tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGc      <  1:2756240/93‑1 (MQ=255)
                                                                           tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGc      >  1:1246741/1‑93 (MQ=255)
                                                                           tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGc      >  1:336362/1‑93 (MQ=255)
                                                                            gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCg     >  2:4258583/1‑93 (MQ=255)
                                                                            gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCg     >  2:1980376/1‑93 (MQ=255)
                                                                            gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGTAAGGCTGCATTATGGGCAGTGCGCACCAGCg     >  2:3586240/1‑93 (MQ=255)
                                                                             aaGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGt    <  2:1341619/93‑1 (MQ=255)
                                                                              aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   >  2:1623503/1‑93 (MQ=255)
                                                                              aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   >  2:2104761/1‑93 (MQ=255)
                                                                              aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  2:4469741/93‑1 (MQ=255)
                                                                              aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  1:2805984/93‑1 (MQ=255)
                                                                              aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  2:3541119/93‑1 (MQ=255)
                                                                              aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc   <  2:1727104/93‑1 (MQ=255)
                                                                               ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  <  1:4531157/93‑1 (MQ=255)
                                                                               ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  >  2:4164986/1‑93 (MQ=255)
                                                                               ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt  >  2:399184/1‑93 (MQ=255)
                                                                                     |                                                                                         
TCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCT  >  NC_000913/4008346‑4008503

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGAC  >  NC_000913/4008371‑4008510
                                                         |                                                                                               
canaagnncatgcgctggaagnngtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT                                                             <  SRR4302535.4590623/56‑1 (MQ=25)
caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT                                                             >  SRR4302535.2493380/3‑94 (MQ=18)
caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT                                                             >  SRR4302535.2044988/3‑94 (MQ=18)
caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT                                                             >  SRR4302535.4497158/3‑94 (MQ=18)
caAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAT                                                             <  SRR4302535.2250476/92‑1 (MQ=18)
 AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG                                                            >  SRR4302535.1057512/1‑94 (MQ=19)
 AAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATG                                                            >  SRR4302535.4621718/1‑94 (MQ=19)
  AAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGC                                                           >  SRR4302535.3891097/1‑94 (MQ=21)
   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT                                                          <  SRR4302535.2927986/94‑1 (MQ=21)
   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT                                                          >  SRR4302535.4686452/1‑94 (MQ=21)
   AAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCT                                                          <  SRR4302535.4480591/94‑1 (MQ=21)
    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG                                                         <  SRR4302535.3137592/94‑1 (MQ=23)
    AGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTG                                                         <  SRR4302535.4663384/94‑1 (MQ=23)
     GGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGT                                                        >  SRR4302535.4137075/1‑94 (MQ=25)
       CCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCG                                                      <  SRR4302535.3291791/94‑1 (MQ=28)
        catgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGA                                                     >  SRR4302535.68451/31‑94 (MQ=42)
         atgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGAT                                                    <  SRR4302535.670392/65‑1 (MQ=45)
         atgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGAG                                                    >  SRR4302535.708396/30‑94 (MQ=42)
          tgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATG                                                   <  SRR4302535.3824190/66‑1 (MQ=47)
          tgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATG                                                   <  SRR4302535.3836667/66‑1 (MQ=47)
          tgcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATG                                                   <  SRR4302535.3607013/66‑1 (MQ=47)
           gcgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGG                                                  <  SRR4302535.1107420/67‑1 (MQ=49)
            cgctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC                                                 <  SRR4302535.3844523/68‑1 (MQ=51)
              ctggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGG                                              >  SRR4302535.4438932/25‑94 (MQ=55)
               tggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG                                             <  SRR4302535.3256569/71‑1 (MQ=57)
               tggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG                                             >  SRR4302535.1508168/24‑94 (MQ=57)
               tggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG                                             >  SRR4302535.923139/24‑94 (MQ=57)
                ggaagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑G                                          >  SRR4302535.4113452/23‑94 (MQ=56)
                  aagcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGTCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑GAA                                        >  SRR4302535.3812459/21‑94 (MQ=42)
                    gcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAG                                       <  SRR4302535.4233280/76‑1 (MQ=53)
                    gcggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAG                                       <  SRR4302535.3360272/76‑1 (MQ=53)
                     cggtgtcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaagg                                       <  SRR4302535.1449641/77‑7 (MQ=57)
                      ggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggc                                      >  SRR4302535.4371501/17‑87 (MQ=57)
                      ggtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggc                                      >  SRR4302535.843260/17‑87 (MQ=57)
                       gtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggct                                     <  SRR4302535.2857901/79‑9 (MQ=57)
                       gtggcgaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGgaaaggct                                     <  SRR4302535.3930632/79‑9 (MQ=57)
                           ccaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCAT                                 <  SRR4302535.2780208/83‑1 (MQ=49)
                            gaagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATT                                >  SRR4302535.2628570/11‑94 (MQ=49)
                             aagaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTA                               >  SRR4302535.4576429/10‑94 (MQ=48)
                              agaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTAT                              <  SRR4302535.2354039/86‑1 (MQ=48)
                               gaaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATG                             <  SRR4302535.1021142/87‑1 (MQ=49)
                                aaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGG                            <  SRR4302535.3344291/88‑1 (MQ=51)
                                aaactgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGG                            >  SRR4302535.1760338/7‑94 (MQ=51)
                                 AACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGG                           <  SRR4302535.3071604/94‑1 (MQ=50)
                                   ctgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCA                         >  SRR4302535.1038215/4‑94 (MQ=57)
                                    tgGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAG                        <  SRR4302535.4002565/92‑1 (MQ=58)
                                      GGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTG                      <  SRR4302535.2524617/94‑1 (MQ=60)
                                       GCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGC                     >  SRR4302535.2050421/1‑94 (MQ=60)
                                        CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcg                    <  SRR4302535.3363861/94‑3 (MQ=60)
                                        CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcg                    >  SRR4302535.1698782/1‑92 (MQ=60)
                                        CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcg                    >  SRR4302535.694210/1‑92 (MQ=60)
                                         TACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcgc                   >  SRR4302535.94673/1‑91 (MQ=60)
                                           CAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcgcac                 >  SRR4302535.4034112/1‑89 (MQ=60)
                                            agcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACC                >  SRR4302535.4362005/4‑94 (MQ=60)
                                             gcCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGcgcaccg               >  SRR4302535.756261/3‑87 (MQ=60)
                                             CCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAG              >  SRR4302535.1201421/1‑94 (MQ=60)
                                               CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC             <  SRR4302535.967423/94‑1 (MQ=60)
                                               CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC             >  SRR4302535.1255636/1‑94 (MQ=60)
                                               CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC             >  SRR4302535.338688/1‑94 (MQ=60)
                                               CTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGC             <  SRR4302535.3919415/94‑1 (MQ=60)
                                                TGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCG            <  SRR4302535.2775887/94‑1 (MQ=57)
                                                   AGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCT         <  SRR4302535.2825932/94‑1 (MQ=57)
                                                    GGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTG        <  SRR4302535.4562173/94‑1 (MQ=60)
                                                     GATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGA       >  SRR4302535.57736/1‑94 (MQ=60)
                                                      ATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAA      <  SRR4302535.3356104/94‑1 (MQ=60)
                                                      ATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAA      <  SRR4302535.912953/94‑1 (MQ=60)
                                                         TATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAC  <  SRR4302535.582540/94‑1 (MQ=60)
                                                         |                                                                                               
TCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGAC  >  NC_000913/4008371‑4008510

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