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 F2 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 539,975 A→G D258D (GAT→GAC allC ← allantoate deiminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273539,9750AG100.0% 166.8 / NA 51D258D (GAT→GACallCallantoate deiminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/14);  total (37/14)

AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGC  >  NZ_CP009273/539890‑540059
                                                                                     |                                                                                    
aatcaatGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGc                                                                                  >  2:288272/1‑90 (MQ=255)
aatcaatGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGc                                                                                  >  2:665736/1‑90 (MQ=255)
     atGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGTCAAAGGTCAGAACCAGCGGGTCGCCCATc                                                                             >  2:297797/1‑90 (MQ=255)
          gAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCtttt                                                                        >  2:410795/1‑90 (MQ=255)
          gAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCtttt                                                                        >  1:238810/1‑90 (MQ=255)
          gAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCACCCtttc                                                                        >  1:1387293/1‑89 (MQ=255)
           aaCGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTc                                                                       >  2:411105/1‑90 (MQ=255)
               tGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCtt                                                                   >  1:56626/1‑90 (MQ=255)
                  tcttACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCttttt                                                                <  1:1341205/88‑1 (MQ=255)
                    ttACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcg                                                              >  1:1119717/1‑90 (MQ=255)
                      aCCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcgac                                                            >  1:179752/1‑90 (MQ=255)
                      aCCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcgac                                                            >  2:555816/1‑90 (MQ=255)
                       cccGGCCCCCCATTTTCCCTTTTCGGGCGGGGGTCCTCCTTGCCAAAGGTCAGAACCAGGGGGTCGGCCATCCTTTTCGCTTTTTcgacc                                                           <  2:1387293/90‑1 (MQ=255)
                         cGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcgaccga                                                         >  1:1466272/1‑90 (MQ=255)
                           gCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACt                                                       >  2:186418/1‑90 (MQ=255)
                           gCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACt                                                       >  2:1261748/1‑90 (MQ=255)
                                acaTTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATgg                                                  >  1:568250/1‑90 (MQ=255)
                                     tACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCg                                                                      >  2:49255/1‑65 (MQ=255)
                                     tACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCg                                                                      <  1:49255/65‑1 (MQ=255)
                                     tACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAAt                                             >  2:900478/1‑90 (MQ=255)
                                        cGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCg                                          <  1:761770/90‑1 (MQ=255)
                                        cGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCg                                          >  2:76638/1‑90 (MQ=255)
                                        cGTATTCGGGCGCAGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCg                                          >  1:540222/1‑90 (MQ=255)
                                           aTTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACt                                       >  2:582046/1‑90 (MQ=255)
                                              cGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGaa                                    <  1:158547/90‑1 (MQ=255)
                                               gggCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGaaa                                   <  1:751389/90‑1 (MQ=255)
                                                     ggCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTaa                             <  2:1119717/90‑1 (MQ=255)
                                                     ggCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTaa                             >  1:709279/1‑90 (MQ=255)
                                                       ctctACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                           <  2:568250/90‑1 (MQ=255)
                                                       ctctACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                           <  1:1206924/90‑1 (MQ=255)
                                                              ttGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCa                    >  1:398485/1‑90 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                  >  2:1530439/1‑90 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                  >  2:1419838/1‑90 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                  >  2:1358670/1‑90 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                  >  2:1320332/1‑90 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                  <  1:900478/90‑1 (MQ=255)
                                                                  cAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacgac                >  2:511001/1‑90 (MQ=255)
                                                                   aaaGGTCAGAACCAGCGGGACGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacgacg               >  2:798674/1‑90 (MQ=255)
                                                                     aGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGAt             >  2:1155027/1‑90 (MQ=255)
                                                                         cAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAAcc         >  1:1374158/1‑90 (MQ=255)
                                                                             aCCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATc     <  1:763141/90‑1 (MQ=255)
                                                                              ccAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCg    <  2:56626/90‑1 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                     >  1:175029/1‑55 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                     <  2:175029/55‑1 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                           <  2:613908/65‑1 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                           >  1:613908/1‑65 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc  >  2:506326/1‑90 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc  >  2:416867/1‑90 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc  >  2:662696/1‑90 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc  >  1:699138/1‑90 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACGGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc  >  1:438952/1‑90 (MQ=255)
                                                                                     |                                                                                    
AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGC  >  NZ_CP009273/539890‑540059

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGTGC  >  NZ_CP009273/539890‑540066
                                                                                     |                                                                                           
AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTT                                                                               >  SRR3722078.241723/1‑100 (MQ=60)
AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCACCCTTTC                                                                               >  SRR3722078.1409702/1‑100 (MQ=60)
     ATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTT                                                                          >  SRR3722078.57296/1‑100 (MQ=60)
          GAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCG                                                                     >  SRR3722078.1136811/1‑100 (MQ=60)
            gtgtataagagacaGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGAC                                                                   <  SRR3722078.49830/86‑1 (MQ=60)
            ACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGAC                                                                   >  SRR3722078.181873/1‑100 (MQ=60)
               TGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGA                                                                >  SRR3722078.1490194/1‑100 (MQ=60)
                  TCTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTG                                                             <  SRR3722078.1362719/100‑1 (MQ=60)
                      ACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGG                                                         >  SRR3722078.574909/1‑100 (MQ=60)
                              CCACATTTACCGTATTCGGGCGCAGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCG                                                 >  SRR3722078.546597/1‑100 (MQ=60)
                                        CGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCG                                       <  SRR3722078.771561/100‑1 (MQ=60)
                                           ATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAA                                    >  SRR3722078.717911/1‑100 (MQ=60)
                                              CGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACT                                 <  SRR3722078.160397/100‑1 (MQ=60)
                                               GGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTG                                <  SRR3722078.760988/100‑1 (MQ=60)
                                                    CGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCA                           >  SRR3722078.403319/1‑100 (MQ=60)
                                                       CTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGA                        <  SRR3722078.1225876/100‑1 (MQ=60)
                                                               TGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACC                >  SRR3722078.1396303/1‑100 (MQ=60)
                                                                GCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCC               <  SRR3722078.912568/100‑1 (MQ=60)
                                                                      GGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACctgtctcttatacacatctgacgct         >  SRR3722078.177078/1‑75 (MQ=60)
                                                                      GGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGC         >  SRR3722078.707551/1‑100 (MQ=60)
                                                                      GGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACGGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGC         >  SRR3722078.444221/1‑100 (MQ=60)
                                                                      GGACAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACtgtctcttatacac         >  SRR3722078.620989/1‑86 (MQ=60)
                                                                           GAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGT    >  SRR3722078.105575/1‑100 (MQ=60)
                                                                             ACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGTGC  <  SRR3722078.772959/100‑1 (MQ=60)
                                                                                     |                                                                                           
AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGCGTGGTGC  >  NZ_CP009273/539890‑540066

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