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 F4 I88 R1 217 30.5 1742422 95.9% 1670982 87.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 675,373 A→G E137G (GAA→GGA)  ybeU → DUF1266 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273675,3730AG100.0% 140.9 / NA 37E137G (GAA→GGA) ybeUDUF1266 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/21);  total (16/21)

GCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATT  >  NZ_CP009273/675286‑675459
                                                                                       |                                                                                      
gCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGggaga                                                                                     >  1:576927/1‑91 (MQ=255)
 cTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGggagaa                                                                                    >  2:79592/1‑91 (MQ=255)
   gtgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGagaaga                                                                                  >  2:841494/1‑91 (MQ=255)
    tgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAagag                                                                                 <  2:202979/91‑1 (MQ=255)
      gggAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAagagag                                                                               >  1:272123/1‑91 (MQ=255)
        gAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGtt                                                                             >  1:464207/1‑91 (MQ=255)
        gAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGtt                                                                             >  2:588938/1‑91 (MQ=255)
        gAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGtt                                                                             >  1:503787/1‑91 (MQ=255)
        gAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTAAACAAGTGGCTGACGGGAGAAGAGAGtt                                                                             >  1:412625/1‑91 (MQ=255)
          aGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTa                                                                           >  2:765214/1‑91 (MQ=255)
             gggAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATgg                                                                        >  2:468658/1‑91 (MQ=255)
                        ccTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATcgcg                                                             <  2:210589/91‑1 (MQ=255)
                         cTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATcgcgc                                                            <  1:355511/91‑1 (MQ=255)
                         cTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATcgcgc                                                            >  2:748001/1‑91 (MQ=255)
                           gggATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGt                                                          >  2:867821/1‑91 (MQ=255)
                            ggATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTc                                                         <  2:525751/91‑1 (MQ=255)
                            ggATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTc                                                         <  1:150359/91‑1 (MQ=255)
                            ggATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTc                                                         <  1:765214/91‑1 (MQ=255)
                             gATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCt                                                        <  2:576927/91‑1 (MQ=255)
                                    ttCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACg                                                 <  1:748001/91‑1 (MQ=255)
                                         aTGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTc                                            <  2:464207/91‑1 (MQ=255)
                                          tGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTATATGTACGAGCTCa                                           <  1:372068/91‑1 (MQ=255)
                                           ggggTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAc                                          >  1:501302/1‑91 (MQ=255)
                                            gggTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAcc                                         <  1:559981/91‑1 (MQ=255)
                                             ggTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCa                                        <  1:468658/91‑1 (MQ=255)
                                             ggTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCa                                        <  2:680181/91‑1 (MQ=255)
                                                ttCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCACCCTTCTATGGACGGGCTCACCATTa                                     >  2:551704/1‑91 (MQ=255)
                                                 tCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTAc                                    <  1:79592/91‑1 (MQ=255)
                                                 tCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTAc                                    <  1:667461/91‑1 (MQ=255)
                                                     tCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATc                                <  2:409422/91‑1 (MQ=255)
                                                      cACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGGTTATGGCTGCAATCGCGGGGCTATGTAAGAGCTCAACATTAACATca                               >  2:717612/1‑91 (MQ=255)
                                                         gggTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAg                            >  1:553582/1‑91 (MQ=255)
                                                          ggTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGc                           >  1:728811/1‑91 (MQ=255)
                                                                        acaaGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGAt                      <  2:714011/82‑1 (MQ=255)
                                                                        acaaGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGAt                      >  1:714011/1‑82 (MQ=255)
                                                                                cTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCt     >  1:390227/1‑91 (MQ=255)
                                                                                  gACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTAt   <  2:553582/91‑1 (MQ=255)
                                                                                  gACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTAt   <  2:490865/91‑1 (MQ=255)
                                                                                  gACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTAt   >  1:248356/1‑91 (MQ=255)
                                                                                  gACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTAt   <  1:140422/91‑1 (MQ=255)
                                                                                   aCGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTAtt  <  2:228805/91‑1 (MQ=255)
                                                                                       |                                                                                      
GCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATT  >  NZ_CP009273/675286‑675459

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ 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
ACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGAC  >  NZ_CP009273/675276‑675468
                                                                                                 |                                                                                               
ACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGA                                                                                              >  SRR3722118.589148/1‑101 (MQ=60)
      ATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAG                                                                                        >  SRR3722118.277691/1‑101 (MQ=60)
        GTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTT                                                                                      >  SRR3722118.474007/1‑101 (MQ=60)
        GTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTT                                                                                      >  SRR3722118.514464/1‑101 (MQ=60)
        GTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTAAACAAGTGGCTGACGGGAGAAGAGAGTT                                                                                      >  SRR3722118.421306/1‑101 (MQ=60)
                                   CTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTAC                                                           <  SRR3722118.362903/101‑1 (MQ=60)
                                      GGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAG                                                        <  SRR3722118.782014/101‑1 (MQ=60)
                                      GGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAG                                                        <  SRR3722118.153248/101‑1 (MQ=60)
                                           ATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAC                                                   >  SRR3722118.511939/1‑101 (MQ=60)
                                              TTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCAT                                                <  SRR3722118.764426/101‑1 (MQ=60)
                                                    TGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTATATGTACGAGCTCACCATTACTAT                                          <  SRR3722118.379843/101‑1 (MQ=60)
                                                      GGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCA                                        <  SRR3722118.571848/101‑1 (MQ=60)
                                                       GGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCAC                                       <  SRR3722118.478545/101‑1 (MQ=60)
                                                         TTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAG                                     >  SRR3722118.565322/1‑101 (MQ=60)
                                                          TTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGC                                    >  SRR3722118.744820/1‑101 (MQ=60)
                                                           TCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCT                                   <  SRR3722118.81065/101‑1 (MQ=60)
                                                           TCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCT                                   <  SRR3722118.681921/101‑1 (MQ=60)
                                                                        GGGGTGCTCAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTT                      >  SRR3722118.729680/1‑101 (MQ=60)
                                                                                CAACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCT              >  SRR3722118.398420/1‑101 (MQ=60)
                                                                                  ACAAGTGGCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTAT            >  SRR3722118.253390/1‑101 (MQ=60)
                                                                                         GCTGACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAG     >  SRR3722118.132838/1‑101 (MQ=60)
                                                                                            GACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATctgtctcttatacacatctccgagcccacgagacggactcctatctcgtatgccgtcttct  >  SRR3722118.696197/1‑40 (MQ=60)
                                                                                            GACGGGAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGAC  <  SRR3722118.143104/101‑1 (MQ=60)
                                                                                                 |                                                                                               
ACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCACTATTTTTCCGCCTATTCGTTAGGAC  >  NZ_CP009273/675276‑675468

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