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 I206 R1 206 23.7 1314362 96.8% 1272302 86.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,174,381 T→G C270G (TGC→GGC)  yegU → ADP‑ribosylglycohydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,174,3810TG100.0% 88.4 / NA 28C270G (TGC→GGC) yegUADP‑ribosylglycohydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/13);  total (15/13)

GTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGC  >  NZ_CP009273/2174297‑2174464
                                                                                    |                                                                                   
gTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGcc                                                                                >  1:186383/1‑90 (MQ=255)
   gTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTc                                                                             >  1:19195/1‑90 (MQ=255)
    ttGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTcc                                                                            >  1:559474/1‑90 (MQ=255)
      ggCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCtg                                                                          <  1:68211/90‑1 (MQ=255)
                      cTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcg                                                          >  1:203505/1‑90 (MQ=255)
                      cTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcg                                                          >  1:90933/1‑90 (MQ=255)
                        aTTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTggcggc                                                        >  2:200152/1‑90 (MQ=255)
                             gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac                                                   >  2:298538/1‑90 (MQ=255)
                             gTCCGTTCCGTGCGCCATTGCGCTGGTTGAACCGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGacac                                                   >  1:656204/1‑90 (MQ=255)
                                cgttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGa                                                >  2:469239/1‑90 (MQ=255)
                                  ttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGaca                                              >  2:472371/1‑90 (MQ=255)
                                  ttccgtGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGaca                                              <  2:254578/90‑1 (MQ=255)
                                       tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                         <  2:186383/90‑1 (MQ=255)
                                       tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                         <  1:469239/90‑1 (MQ=255)
                                       tGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                         <  1:213117/90‑1 (MQ=255)
                                       tGCGCCATCGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATc                                         <  2:205502/90‑1 (MQ=255)
                                                   cTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCg                             <  2:19195/90‑1 (MQ=255)
                                                     ggTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGAc                           >  1:113184/1‑90 (MQ=255)
                                                           aCTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTTTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAAt                     >  2:301475/1‑90 (MQ=255)
                                                                        gACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTg                                                             >  2:358658/1‑37 (MQ=255)
                                                                        gACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTg                                                             <  1:358658/37‑1 (MQ=255)
                                                                        gACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTg        >  2:529991/1‑90 (MQ=255)
                                                                           ccGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCAt     <  1:533908/90‑1 (MQ=255)
                                                                            ccaaTCGCGGCGCCGTCCGGTGCGCTAACCTTGGCGGCGACACCGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg    <  2:575409/88‑1 (MQ=255)
                                                                            cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg    <  1:627228/90‑1 (MQ=255)
                                                                            cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg    <  2:250450/90‑1 (MQ=255)
                                                                            cGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATg    <  1:523940/90‑1 (MQ=255)
                                                                              aaTCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGc  >  2:427730/1‑90 (MQ=255)
                                                                                    |                                                                                   
GTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGC  >  NZ_CP009273/2174297‑2174464

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ 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
GCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGC  >  NZ_CP009273/2174287‑2174472
                                                                                              |                                                                                           
GCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCC                                                                                        >  SRR3722088.188254/1‑100 (MQ=60)
   TTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTC                                                                                     >  SRR3722088.19399/1‑100 (MQ=60)
    TACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCC                                                                                    <  SRR3722088.384249/100‑1 (MQ=60)
    TACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCC                                                                                    >  SRR3722088.566352/1‑100 (MQ=60)
                GGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACC                                                                        <  SRR3722088.68921/100‑1 (MQ=60)
                      GGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCG                                                                  >  SRR3722088.205543/1‑100 (MQ=60)
                      GGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCG                                                                  >  SRR3722088.91854/1‑100 (MQ=60)
                             tctagtacggtctcgtgggctcggagatgtgtataagagacagGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACAC                                                           <  SRR3722088.362459/57‑1 (MQ=60)
                             GCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACCGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACAC                                                           >  SRR3722088.664331/1‑100 (MQ=60)
                                                 TGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG                                       <  SRR3722088.215256/100‑1 (MQ=60)
                                                 TGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGG                                       <  SRR3722088.474760/100‑1 (MQ=60)
                                                     CCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGAC                                   >  SRR3722088.114362/1‑100 (MQ=60)
                                                                                     CCGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACG   <  SRR3722088.540387/100‑1 (MQ=60)
                                                                                      CGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGC  <  SRR3722088.530278/100‑1 (MQ=60)
                                                                                      CGAATCGCGGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGC  <  SRR3722088.634982/100‑1 (MQ=60)
                                                                                              |                                                                                           
GCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGCGACACAGACACCATCGGTGCTATGGCGACGGCAATTTGCGGCGCGTTGCATGGCGTTAACGC  >  NZ_CP009273/2174287‑2174472

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