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 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,161,702 C→G A162G (GCA→GGA)  btuB → BtuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,7020CG99.2% 407.4 / NA 124A162G (GCA→GGA) btuBBtuB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (71/52);  total (71/53)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATT  >  NC_000913/4161612‑4161791
                                                                                          |                                                                                         
ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTcaggag                                                                                         <  1:3329798/93‑1 (MQ=255)
ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTcaggag                                                                                         <  1:457714/93‑1 (MQ=255)
 gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                                        <  2:3527856/93‑1 (MQ=255)
 gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCaggagg                                                                                        >  2:4459456/1‑93 (MQ=255)
 gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCCCGATGAACCCGGAACGGAAATTTCaggagg                                                                                        >  2:4042789/1‑93 (MQ=255)
  ccGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAggg                                                                                       <  1:3435815/93‑1 (MQ=255)
   cgcgCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGt                                                                                      >  1:1266860/1‑93 (MQ=255)
   cgcgCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGt                                                                                      <  1:1872263/93‑1 (MQ=255)
    gcgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGTGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                                     <  2:1417288/93‑1 (MQ=255)
     cgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgg                                                                                    >  1:1208887/1‑93 (MQ=255)
        tCCGCTGTTTATGGTTCCTATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                                                 >  1:2385063/1‑93 (MQ=255)
        tCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                                                 <  1:374696/93‑1 (MQ=255)
         ccGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                <  1:4130223/93‑1 (MQ=255)
         ccGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                >  2:2564272/1‑93 (MQ=255)
         ccGCTGTTTATGGTTCCGATGCAATACGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                >  1:4726386/1‑93 (MQ=255)
          cGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                                               >  2:3452210/1‑93 (MQ=255)
          cGCTGTTTATGGGTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                                               >  2:4026769/1‑93 (MQ=255)
            cTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCa                                                                             >  1:4401328/1‑93 (MQ=255)
            cTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCa                                                                             >  2:151689/1‑93 (MQ=255)
             tGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCaa                                                                            >  2:4291733/1‑93 (MQ=255)
               tttATGGTTCCGATGCAATAGTCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                          >  1:834945/1‑93 (MQ=255)
               tttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                          >  2:3716063/1‑93 (MQ=255)
               tttATGGTTCCGATGCAATAGGCGGGGCGGTGAATATCATCACGACCCGCGATGAACCCGGAACGGAAATTTCAGGAGGGCGGGGAAGCAAtg                                                                          >  1:271924/1‑92 (MQ=255)
                ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                         <  1:2918584/93‑1 (MQ=255)
                ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                         <  2:4382103/93‑1 (MQ=255)
                 tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                        >  2:289214/1‑93 (MQ=255)
                 tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                        <  1:3911622/93‑1 (MQ=255)
                 tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                        >  1:844277/1‑93 (MQ=255)
                  aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGtt                                                                       >  1:2083885/1‑93 (MQ=255)
                  aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGtt                                                                       >  2:1267253/1‑93 (MQ=255)
                   tGGTTGCGATTCAATAGGCGGGGTGGGTAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTa                                                                      <  1:1422160/93‑1 (MQ=255)
                   tGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTa                                                                      <  1:385570/93‑1 (MQ=255)
                     gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                                    >  1:2064161/1‑93 (MQ=255)
                     gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                                    >  2:4055095/1‑93 (MQ=255)
                      ttCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCa                                                                   >  1:903246/1‑93 (MQ=255)
                       tCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAg                                                                  >  2:2660095/1‑93 (MQ=255)
                          gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAAc                                                               >  1:541639/1‑93 (MQ=255)
                           aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACt                                                              >  2:1455964/1‑93 (MQ=255)
                            tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                                             >  2:1282922/1‑93 (MQ=255)
                            tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                                             >  2:139552/1‑93 (MQ=255)
                            tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                                             <  1:3633605/93‑1 (MQ=255)
                             gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                            <  1:2825158/93‑1 (MQ=255)
                             gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                            >  1:2632241/1‑93 (MQ=255)
                             gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                            <  2:415145/93‑1 (MQ=255)
                              cAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                                           >  1:4001115/1‑93 (MQ=255)
                               aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatga                                                          <  2:336841/93‑1 (MQ=255)
                               aaTAGGCGGGGTGGTGAAAATCATCATGACGCGCGATGAACCCGGAACGGACATTTCAGGAGGGTGGGGAAACCATAGTTATAAGAACTatga                                                          >  2:3869327/1‑93 (MQ=255)
                                aTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAATTatgag                                                         >  1:201604/1‑92 (MQ=255)
                                aTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAGGCAATAGTTATCAGAACTatgat                                                         >  2:1260502/1‑93 (MQ=255)
                                 tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                                        <  2:830484/93‑1 (MQ=255)
                                 tAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                                        <  1:436648/93‑1 (MQ=255)
                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                      <  1:1296270/93‑1 (MQ=255)
                                   ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                      <  2:3179852/93‑1 (MQ=255)
                                   ggCGGGGTGCTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                      <  2:2253809/93‑1 (MQ=255)
                                      ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                   <  1:4506155/93‑1 (MQ=255)
                                       gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                  <  2:1084315/93‑1 (MQ=255)
                                        ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                                                 >  1:4381188/1‑93 (MQ=255)
                                        ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                                                 <  1:633382/93‑1 (MQ=255)
                                          tggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                                               >  2:125369/1‑93 (MQ=255)
                                          tggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                                               >  1:927174/1‑93 (MQ=255)
                                               aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaa                                          >  2:3231410/1‑93 (MQ=255)
                                               aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaa                                          >  1:638093/1‑93 (MQ=255)
                                                 tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                        >  2:3894274/1‑93 (MQ=255)
                                                 tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                        >  1:450838/1‑93 (MQ=255)
                                                 tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGCGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                        >  2:3609612/1‑93 (MQ=255)
                                                  atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaacaa                                       <  2:2327337/93‑1 (MQ=255)
                                                  atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaacaa                                       >  2:3554359/1‑93 (MQ=255)
                                                  atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaacaa                                       >  1:2271745/1‑93 (MQ=255)
                                                   tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                      <  1:2202912/93‑1 (MQ=255)
                                                   gcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                      <  2:2499392/92‑1 (MQ=255)
                                                      tcaCGACGCGCGCTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAGGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgg                                   <  2:4213213/93‑1 (MQ=255)
                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgg                                   <  1:4168992/93‑1 (MQ=255)
                                                      tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgg                                   >  1:976275/1‑93 (MQ=255)
                                                        acgtcgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg                                 >  2:3980502/1‑93 (MQ=255)
                                                        acgacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg                                 >  1:720963/1‑93 (MQ=255)
                                                             gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa                            <  2:505273/93‑1 (MQ=255)
                                                             gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa                            <  1:1455062/93‑1 (MQ=255)
                                                               gcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGa                          >  1:2094398/1‑93 (MQ=255)
                                                                cgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                         <  2:2895205/93‑1 (MQ=255)
                                                                  aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac                       >  1:1511546/1‑93 (MQ=255)
                                                                    gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACgg                     <  2:4603688/93‑1 (MQ=255)
                                                                      aCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGt                   >  1:28945/1‑93 (MQ=255)
                                                                        ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTaa                 >  2:3836624/1‑93 (MQ=255)
                                                                        ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTaa                 >  1:2790513/1‑93 (MQ=255)
                                                                         cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTGTCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                <  2:1937871/93‑1 (MQ=255)
                                                                         cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                >  2:2943915/1‑93 (MQ=255)
                                                                           gaactgaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAGAAGACACGGGTAACGc              <  2:1898643/93‑1 (MQ=255)
                                                                           gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc              >  2:2445879/1‑93 (MQ=255)
                                                                           gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc              <  1:3623686/93‑1 (MQ=255)
                                                                            aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt             <  2:370453/93‑1 (MQ=255)
                                                                            aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt             <  2:2758909/93‑1 (MQ=255)
                                                                              cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGt           >  2:1989316/1‑93 (MQ=255)
                                                                              cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGt           >  2:3505314/1‑93 (MQ=255)
                                                                              cggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGt           >  2:3341454/1‑93 (MQ=255)
                                                                               ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtt          <  2:2308638/93‑1 (MQ=255)
                                                                                   aTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGc      <  1:2754331/93‑1 (MQ=255)
                                                                                    tttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCg     >  2:457636/1‑93 (MQ=255)
                                                                                       caggagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTACCGCTGTTGGGCGAtt  <  2:3545568/93‑1 (MQ=255)
                                                                                       caggagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt  <  1:1716638/93‑1 (MQ=255)
                                                                                       caggagGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt  <  1:900446/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATT  >  NC_000913/4161612‑4161791

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGC  >  NC_000913/4161611‑4161795
                                                                                           |                                                                                             
GGGCCGCGCTGCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAG                                                                                             <  SRR4302544.1866230/94‑1 (MQ=60)
 GGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGG                                                                                            <  SRR4302544.462220/94‑1 (MQ=60)
 GGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGG                                                                                            <  SRR4302544.3361508/94‑1 (MQ=60)
   CCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGT                                                                                          >  SRR4302544.1279315/1‑94 (MQ=60)
   CCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGT                                                                                          <  SRR4302544.3468663/94‑1 (MQ=60)
    CGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTG                                                                                         <  SRR4302544.1890359/94‑1 (MQ=60)
     GCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGG                                                                                        >  SRR4302544.1220861/1‑94 (MQ=60)
       CCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGG                                                                                      <  SRR4302544.3477736/94‑1 (MQ=60)
        CTCCGCTGTTTATGGTTCCTATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGA                                                                                     >  SRR4302544.2407651/1‑94 (MQ=60)
         TCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAA                                                                                    <  SRR4302544.378521/94‑1 (MQ=60)
         TCCGCTGTTTATGGTTCCGATGCAATACGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAA                                                                                    >  SRR4302544.4773034/1‑94 (MQ=60)
          CCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAG                                                                                   <  SRR4302544.4171165/94‑1 (MQ=60)
            GCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCA                                                                                 >  SRR4302544.4445801/1‑94 (MQ=60)
             CTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAA                                                                                >  SRR4302544.1709272/1‑94 (MQ=60)
              TGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAT                                                                               >  SRR4302544.3102132/1‑94 (MQ=60)
               GTTTATGGTTCCGATGCAATAGTCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATA                                                                              >  SRR4302544.843172/1‑94 (MQ=60)
               GTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTTGGGAAGCAATA                                                                              >  SRR4302544.4572171/1‑94 (MQ=60)
               GTTTATGGTTCCGATGCAATAGGCGGGGCGGTGAATATCATCACGACCCGCGATGAACCCGGAACGGAAATTTCAGGAGGGCGGGGAAGCAATG                                                                              >  SRR4302544.274786/1‑94 (MQ=60)
                 TTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGT                                                                            >  SRR4302544.852611/1‑94 (MQ=60)
                 TTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGT                                                                            <  SRR4302544.2946644/94‑1 (MQ=60)
                  TATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTT                                                                           >  SRR4302544.2103926/1‑94 (MQ=60)
                  TATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTT                                                                           <  SRR4302544.3950480/94‑1 (MQ=60)
                  TATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTT                                                                           <  SRR4302544.1933714/94‑1 (MQ=60)
                   ATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTA                                                                          <  SRR4302544.967432/94‑1 (MQ=60)
                    TGGTTGCGATTCAATAGGCGGGGTGGGTAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTAT                                                                         <  SRR4302544.1436314/94‑1 (MQ=60)
                    TGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTAT                                                                         <  SRR4302544.389489/94‑1 (MQ=60)
                     GGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATC                                                                        >  SRR4302544.2084002/1‑94 (MQ=60)
                      GTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCA                                                                       >  SRR4302544.912185/1‑94 (MQ=60)
                          CGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAAC                                                                   >  SRR4302544.546971/1‑94 (MQ=60)
                             TGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTAT                                                                >  SRR4302544.2657469/1‑94 (MQ=60)
                             TGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTAT                                                                <  SRR4302544.3669837/94‑1 (MQ=60)
                              GCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATG                                                               <  SRR4302544.2852454/94‑1 (MQ=60)
                              GCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATG                                                               >  SRR4302544.4040937/1‑94 (MQ=60)
                               CAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGA                                                              >  SRR4302544.1317867/1‑94 (MQ=60)
                                AATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAATTATGAG                                                             >  SRR4302544.203641/1‑94 (MQ=60)
                                  TAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGT                                                           >  SRR4302544.866696/1‑94 (MQ=60)
                                  TAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGT                                                           <  SRR4302544.441042/94‑1 (MQ=60)
                                    GGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCT                                                         <  SRR4302544.1309038/94‑1 (MQ=60)
                                       GGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTA                                                      >  SRR4302544.4518126/1‑94 (MQ=60)
                                       GGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTA                                                      <  SRR4302544.4551341/94‑1 (MQ=60)
                                        GGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAC                                                     >  SRR4302544.4425569/1‑94 (MQ=60)
                                        GGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAC                                                     <  SRR4302544.73913/94‑1 (MQ=60)
                                         GGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACG                                                    <  SRR4302544.639620/94‑1 (MQ=60)
                                          GTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGC                                                   >  SRR4302544.936311/1‑94 (MQ=60)
                                               GAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAA                                              >  SRR4302544.2104521/1‑94 (MQ=60)
                                               GAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAA                                              >  SRR4302544.644382/1‑94 (MQ=60)
                                                 ATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACA                                            >  SRR4302544.455310/1‑94 (MQ=60)
                                                  TATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAA                                           >  SRR4302544.2293126/1‑94 (MQ=60)
                                                    TCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACT                                         <  SRR4302544.2223663/94‑1 (MQ=60)
                                                     CATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTG                                        >  SRR4302544.1396875/1‑94 (MQ=60)
                                                      ATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGG                                       >  SRR4302544.985850/1‑94 (MQ=60)
                                                       TCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGG                                      <  SRR4302544.4210294/94‑1 (MQ=60)
                                                        cagcacccgcgtcGAACCCGGAACGGAAATTTCAGGAGGGTGGTGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGG                                     <  SRR4302544.3125344/81‑1 (MQ=60)
                                                        CACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGG                                     >  SRR4302544.728071/1‑94 (MQ=60)
                                                              GCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAG                               <  SRR4302544.1469516/94‑1 (MQ=60)
                                                               CGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGA                              >  SRR4302544.2114553/1‑94 (MQ=60)
                                                                  GATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACAC                           >  SRR4302544.1526526/1‑94 (MQ=60)
                                                                      AACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGT                       >  SRR4302544.29221/1‑94 (MQ=60)
                                                                       ACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTA                      <  SRR4302544.2822513/94‑1 (MQ=60)
                                                                        CCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAA                     >  SRR4302544.2817381/1‑94 (MQ=60)
                                                                            GAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCT                 <  SRR4302544.3659793/94‑1 (MQ=60)
                                                                            GAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCT                 >  SRR4302544.2349936/1‑94 (MQ=60)
                                                                              ACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGT               <  SRR4302544.3830115/94‑1 (MQ=60)
                                                                              ACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGT               >  SRR4302544.568489/1‑94 (MQ=60)
                                                                               CGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTT              <  SRR4302544.4004998/94‑1 (MQ=60)
                                                                                GGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTG             >  SRR4302544.2717405/1‑94 (MQ=60)
                                                                                    ATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTTGGGGATAAGACACGGGTAACGCTGTTGGGCG         >  SRR4302544.4358198/1‑94 (MQ=60)
                                                                                    ATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCG         <  SRR4302544.2780844/94‑1 (MQ=60)
                                                                                        CAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTA     <  SRR4302544.909356/94‑1 (MQ=60)
                                                                                        CAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTA     >  SRR4302544.865291/1‑94 (MQ=60)
                                                                                        CAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTA     <  SRR4302544.1733514/94‑1 (MQ=60)
                                                                                          GGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATG   <  SRR4302544.4504346/94‑1 (MQ=60)
                                                                                           GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGC  >  SRR4302544.4280139/1‑94 (MQ=60)
                                                                                           |                                                                                             
GGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGC  >  NC_000913/4161611‑4161795

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