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 2,298,561 T→C H738R (CAC→CGC)  napA ← putative nitrate reductase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,298,5610TC99.0% 287.1 / ‑4.7 97H738R (CAC→CGC) napAputative nitrate reductase 3
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (55/41);  total (55/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

Only 100 of 113 total aligned reads displayed.

GAAACCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGTGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTCATCC  >  NC_000913/2298471‑2298650
                                                                                          |                                                                                         
gaAACCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAggcgcggcgc                                                                                         >  2:1322709/1‑93 (MQ=255)
 aaaCCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCa                                                                                        >  1:3359211/1‑93 (MQ=255)
  aaCCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAg                                                                                       >  2:1119265/1‑93 (MQ=255)
  aaCCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAg                                                                                       >  1:275991/1‑93 (MQ=255)
  aaCCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAg                                                                                       >  2:4276172/1‑93 (MQ=255)
   aCCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGt                                                                                      >  1:3046536/1‑93 (MQ=255)
    ccACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGtt                                                                                     >  2:2988619/1‑93 (MQ=255)
     cactttcactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTc                                                                                    >  1:129378/1‑93 (MQ=255)
      cctttcgctttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCt                                                                                   <  2:737446/92‑1 (MQ=255)
      actttcactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCt                                                                                   >  1:3264518/1‑93 (MQ=255)
      actttcactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCt                                                                                   >  1:3395035/1‑93 (MQ=255)
       ctttcactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTg                                                                                  >  1:2734294/1‑93 (MQ=255)
       ctttcactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTg                                                                                  <  1:1721346/93‑1 (MQ=255)
        tttcactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGAGCGGCGCAGTTCTgg                                                                                 >  1:4014644/1‑93 (MQ=255)
        tttcactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGCAAGGCGCGGCGCAGTTCTgg                                                                                 >  1:404394/1‑93 (MQ=255)
         ttcactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGt                                                                                <  1:1754650/93‑1 (MQ=255)
         ttcactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGt                                                                                >  1:1716774/1‑93 (MQ=255)
          tcactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTa                                                                               >  2:1801615/1‑93 (MQ=255)
           cactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTac                                                                              <  1:1124772/93‑1 (MQ=255)
           cactttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGAACAGCTTACGCGAAGGCGCGGCGCAGTTCTGGTac                                                                              >  2:2370802/1‑93 (MQ=255)
            tctttGTCGCCAGGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAAAAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTaca                                                                             <  2:2124384/92‑1 (MQ=255)
            actttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTaca                                                                             >  1:1741520/1‑93 (MQ=255)
             ctttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTacac                                                                            <  2:249015/93‑1 (MQ=255)
              tttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACg                                                                           <  2:3396138/93‑1 (MQ=255)
              tttGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACg                                                                           >  2:1969665/1‑93 (MQ=255)
                tGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGc                                                                         >  1:1864630/1‑93 (MQ=255)
                tGTCGCCACCGCGCAGATCGCGCGCTTTCGCATCCAGCCGGGGAAAAAACAGGACCGCTTCAGGGAAAGCGCGGCGCAGTTCTGGTACACggt                                                                         >  2:412175/1‑92 (MQ=255)
                 gTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCg                                                                        >  2:1680704/1‑93 (MQ=255)
                     ccACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTc                                                                    >  2:1463687/1‑93 (MQ=255)
                      cACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCa                                                                   >  2:1450809/1‑93 (MQ=255)
                      cACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCa                                                                   <  1:2685417/93‑1 (MQ=255)
                       aCGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGCCCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCAt                                                                  <  2:3316458/93‑1 (MQ=255)
                       aCGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCAt                                                                  >  1:2342577/1‑93 (MQ=255)
                         ggCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATAc                                                                >  1:3370084/1‑93 (MQ=255)
                          gcgcAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACt                                                               <  1:343392/93‑1 (MQ=255)
                          gcgcAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACt                                                               <  2:1811256/93‑1 (MQ=255)
                            gcAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGc                                                             >  1:1454501/1‑93 (MQ=255)
                            gcAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGc                                                             >  2:2272207/1‑93 (MQ=255)
                            gcAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGc                                                             >  1:2162527/1‑93 (MQ=255)
                              aGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCg                                                           >  2:803441/1‑93 (MQ=255)
                               gATCGCGCGCTTTCGCATCCAGCGGGTGAATACACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTCTTACACGGCGAGTCATACTCCCgc                                                          >  2:3141717/1‑92 (MQ=255)
                               gATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGAGAGTCATACTGCCgg                                                          >  1:4508972/1‑93 (MQ=255)
                               gATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCGTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCgg                                                          >  1:253997/1‑93 (MQ=255)
                                      cgcTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCa                                                   >  2:1534545/1‑93 (MQ=255)
                                      cgcTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCa                                                   <  2:3402532/93‑1 (MQ=255)
                                       gcTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAg                                                  >  2:906187/1‑93 (MQ=255)
                                       gcTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAg                                                  >  1:3134325/1‑93 (MQ=255)
                                       gcTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAg                                                  >  1:62041/1‑93 (MQ=255)
                                       gcTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAg                                                  >  1:212077/1‑93 (MQ=255)
                                       gcTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAg                                                  <  1:2350032/93‑1 (MQ=255)
                                       gcTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAg                                                  >  1:3641100/1‑93 (MQ=255)
                                       gcTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCCAGTCATACTGCCGGTGTGCCAg                                                  >  2:2517508/1‑93 (MQ=255)
                                        cTTTCGGATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGAGCAGTTCTGGTACACGGCGAGTCATACTCACGGTGTGCCAgg                                                 >  2:4449979/1‑92 (MQ=255)
                                         tttCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTg                                                >  2:476685/1‑93 (MQ=255)
                                            cGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTc                                             <  2:1905932/93‑1 (MQ=255)
                                             gCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTcc                                            <  2:3880670/93‑1 (MQ=255)
                                               aTCCAGCGGGGGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGGCATACTGCCGGTGTGCCAGTGCTCCAg                                          >  2:1356330/1‑93 (MQ=255)
                                                tCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGa                                         >  1:3242260/1‑93 (MQ=255)
                                                tCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGa                                         >  1:29249/1‑93 (MQ=255)
                                                 ccAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGaa                                        <  2:3578191/93‑1 (MQ=255)
                                                 ccAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGaa                                        <  1:3357379/93‑1 (MQ=255)
                                                  cAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAAc                                       <  2:1990391/93‑1 (MQ=255)
                                                  cAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAAc                                       >  1:4183771/1‑93 (MQ=255)
                                                  cAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAAc                                       <  2:455098/93‑1 (MQ=255)
                                                    gCGGGTGAATAAACAGGACCGCTGCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAAcgc                                     <  2:1169031/93‑1 (MQ=255)
                                                      gggTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGt                                   >  2:2471487/1‑93 (MQ=255)
                                                       ggTGAATAAACAGGACCGCTTCCGGGAAGGCGGGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTc                                  <  2:222202/93‑1 (MQ=255)
                                                       ggTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTc                                  <  2:3351368/93‑1 (MQ=255)
                                                       ggTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTc                                  <  1:1330556/93‑1 (MQ=255)
                                                         tGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCg                                >  2:3550419/1‑93 (MQ=255)
                                                           aaTAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGt                              >  2:3022990/1‑93 (MQ=255)
                                                           aaTAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGt                              >  1:2944523/1‑93 (MQ=255)
                                                           aaTAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACCCGGCCgtg                              >  2:3313507/1‑91 (MQ=255)
                                                            aTAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTa                             <  2:2877609/93‑1 (MQ=255)
                                                            aTAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTa                             <  1:4115572/93‑1 (MQ=255)
                                                            aTAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTa                             <  2:3255781/93‑1 (MQ=255)
                                                             tAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTag                            >  1:735222/1‑93 (MQ=255)
                                                             tAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTag                            >  2:4013297/1‑93 (MQ=255)
                                                             tAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTag                            >  2:1908126/1‑93 (MQ=255)
                                                                aCAGGACCGCTTCCGGGAAGGCGCGGTGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTagaga                         <  1:4666563/93‑1 (MQ=255)
                                                                aCAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTagaga                         <  1:2568388/93‑1 (MQ=255)
                                                                 cAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTagagag                        <  2:897354/93‑1 (MQ=255)
                                                                 cAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTagagag                        <  2:558579/93‑1 (MQ=255)
                                                                 cAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTagagag                        >  2:4343945/1‑93 (MQ=255)
                                                                   ggACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGcc                      <  2:3380574/93‑1 (MQ=255)
                                                                    gACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCca                     >  1:3522886/1‑93 (MQ=255)
                                                                      ccGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCcaca                   <  1:2890489/93‑1 (MQ=255)
                                                                       cGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAg                  <  1:614399/93‑1 (MQ=255)
                                                                        gCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAgg                 >  2:2532808/1‑93 (MQ=255)
                                                                        gCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAgg                 <  1:4602159/93‑1 (MQ=255)
                                                                        gCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAgg                 >  2:2332322/1‑93 (MQ=255)
                                                                        gCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAgg                 >  2:704943/1‑93 (MQ=255)
                                                                         cTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGt                <  2:2751271/93‑1 (MQ=255)
                                                                          ttCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTc               >  2:3263834/1‑93 (MQ=255)
                                                                          tgccGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTc               <  2:4449668/91‑1 (MQ=255)
                                                                             cGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTa            <  1:3567394/93‑1 (MQ=255)
                                                                             cGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTa            <  1:2310945/93‑1 (MQ=255)
                                                                                   ggcgcggcgcAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTc      <  1:4614596/93‑1 (MQ=255)
                                                                                   ggcgcggcgcAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTc      <  1:638786/93‑1 (MQ=255)
                                                                                       cgGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTCATcc  <  1:4063524/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAAACCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGTGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTCATCC  >  NC_000913/2298471‑2298650

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 26 ≤ 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
CGAGAAACCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGTGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTCATCCGGT  >  NC_000913/2298468‑2298653
                                                                                             |                                                                                            
CGAGAAACCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGC                                                                                              <  SRR4302535.2549662/94‑1 (MQ=60)
 GAGAAACCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCG                                                                                             <  SRR4302535.244423/94‑1 (MQ=60)
   GAAACCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCA                                                                                           >  SRR4302535.3382593/1‑94 (MQ=60)
    AAACCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAG                                                                                          >  SRR4302535.277914/1‑94 (MQ=60)
     AACCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGT                                                                                         >  SRR4302535.3067879/1‑94 (MQ=60)
       CCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTC                                                                                       >  SRR4302535.130289/1‑94 (MQ=60)
        CACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCT                                                                                      >  SRR4302535.3418655/1‑94 (MQ=60)
        CACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCT                                                                                      >  SRR4302535.3287190/1‑94 (MQ=60)
         ACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTG                                                                                     >  SRR4302535.2753805/1‑94 (MQ=60)
          CTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGG                                                                                    <  SRR4302535.1733795/94‑1 (MQ=60)
          CTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGAGCGGCGCAGTTCTGG                                                                                    >  SRR4302535.4042101/1‑94 (MQ=60)
          CTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGCAAGGCGCGGCGCAGTTCTGG                                                                                    >  SRR4302535.407226/1‑94 (MQ=60)
           TTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGT                                                                                   >  SRR4302535.1729195/1‑94 (MQ=60)
            TTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTA                                                                                  <  SRR4302535.1767331/94‑1 (MQ=60)
              CACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACA                                                                                <  SRR4302535.1132797/94‑1 (MQ=60)
              CACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACA                                                                                >  SRR4302535.4524678/1‑94 (MQ=60)
              CACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACA                                                                                >  SRR4302535.1754107/1‑94 (MQ=60)
                  TTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGC                                                                            >  SRR4302535.1878074/1‑94 (MQ=60)
                      CGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGT                                                                        >  SRR4302535.3225637/1‑94 (MQ=60)
                         CACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCAT                                                                     >  SRR4302535.2359340/1‑94 (MQ=60)
                         CACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCAT                                                                     <  SRR4302535.2704619/94‑1 (MQ=60)
                           CGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATAC                                                                   >  SRR4302535.3393531/1‑94 (MQ=60)
                             GCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTG                                                                 <  SRR4302535.345764/94‑1 (MQ=60)
                              CGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGC                                                                >  SRR4302535.1464911/1‑94 (MQ=60)
                              CGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGC                                                                >  SRR4302535.2178044/1‑94 (MQ=60)
                                 AGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGAGAGTCATACTGCCGG                                                             >  SRR4302535.4539912/1‑94 (MQ=60)
                                 AGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCGTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGG                                                             >  SRR4302535.255750/1‑94 (MQ=60)
                                        GCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCA                                                      >  SRR4302535.3532016/1‑94 (MQ=60)
                                         CGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAG                                                     >  SRR4302535.62505/1‑94 (MQ=60)
                                         CGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAG                                                     >  SRR4302535.3666321/1‑94 (MQ=60)
                                         CGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAG                                                     >  SRR4302535.213535/1‑94 (MQ=60)
                                         CGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAG                                                     >  SRR4302535.3156146/1‑94 (MQ=60)
                                          GCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGT                                                    <  SRR4302535.2366879/94‑1 (MQ=60)
                                             TTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTAATACTGCCGGTGTGCCAGTGCT                                                 <  SRR4302535.1853272/94‑1 (MQ=60)
                                                  ATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGA                                            >  SRR4302535.29476/1‑94 (MQ=60)
                                                  ATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGA                                            >  SRR4302535.3264780/1‑94 (MQ=60)
                                                    CCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAAC                                          >  SRR4302535.4212398/1‑94 (MQ=60)
                                                    CCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAAC                                          <  SRR4302535.3380753/94‑1 (MQ=60)
                                                       GCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCG                                       >  SRR4302535.2590399/1‑94 (MQ=60)
                                                          GGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCC                                    <  SRR4302535.1340046/94‑1 (MQ=60)
                                                             GAATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGT                                 >  SRR4302535.2965221/1‑94 (MQ=60)
                                                               ATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAG                               >  SRR4302535.740391/1‑94 (MQ=60)
                                                               ATAAACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAG                               <  SRR4302535.4143743/94‑1 (MQ=60)
                                                                   ACAGGACCGCTTCCGGGAAGGCGCGGTGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAG                           <  SRR4302535.4698106/94‑1 (MQ=60)
                                                                   ACAGGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAG                           <  SRR4302535.2586863/94‑1 (MQ=60)
                                                                      GGACCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCA                        >  SRR4302535.3547352/1‑94 (MQ=60)
                                                                         CCGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAG                     <  SRR4302535.2910939/94‑1 (MQ=60)
                                                                          CGCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGG                    <  SRR4302535.618671/94‑1 (MQ=60)
                                                                           GCTTCCGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGT                   <  SRR4302535.4633456/94‑1 (MQ=60)
                                                                              TGCGGGAGGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATGCTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGT                <  SRR4302535.1004078/94‑1 (MQ=60)
                                                                                CGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTAC              <  SRR4302535.3592130/94‑1 (MQ=60)
                                                                                CGGGAAGGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTAC              <  SRR4302535.2327482/94‑1 (MQ=60)
                                                                                      GGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTCA        <  SRR4302535.4645950/94‑1 (MQ=60)
                                                                                      GGCGCGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTCA        <  SRR4302535.643251/94‑1 (MQ=60)
                                                                                          CGGCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTCATCCG    <  SRR4302535.4091333/94‑1 (MQ=60)
                                                                                            GCGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTCATCCGGT  <  SRR4302535.2209504/94‑1 (MQ=60)
                                                                                             |                                                                                            
CGAGAAACCACTTTCACTTTGTCGCCACGGCGCAGATCGCGCGCTTTCGCATCCAGCGGGTGAATAAACAGGACCGCTTCCGGGAAGGCGCGGTGCAGTTCTGGTACACGGCGAGTCATACTGCCGGTGTGCCAGTGCTCCAGAACGCGTCCGGTAGAGAGCCACAGGTCGTACTCTTCATCCGGT  >  NC_000913/2298468‑2298653

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