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 I229 R1 214 17.4 943020 97.5% 919444 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,859,998 T→C I305V (ATT→GTT)  rpoS ← RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,859,9980TC100.0% 83.5 / NA 27I305V (ATT→GTT) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/13);  total (14/13)

TACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAAC  >  NZ_CP009273/2859915‑2860075
                                                                                   |                                                                             
tacttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCg                                                                         >  2:91273/1‑90 (MQ=255)
tacttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCg                                                                         >  2:78232/1‑90 (MQ=255)
tacttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCg                                                                         >  2:459685/1‑90 (MQ=255)
tacttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCg                                                                         >  2:243759/1‑90 (MQ=255)
    tactTGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                     <  2:356956/90‑1 (MQ=255)
    tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                     <  1:137292/90‑1 (MQ=255)
    tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                     <  1:30960/90‑1 (MQ=255)
    tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                     <  2:383740/90‑1 (MQ=255)
    tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                     <  1:62290/90‑1 (MQ=255)
    tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                     >  1:470864/1‑90 (MQ=255)
         gcgGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTc                                                                >  1:176933/1‑90 (MQ=255)
           ggAGCAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCAc                                                              <  1:379281/90‑1 (MQ=255)
                      tCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCa                                                   >  1:62559/1‑90 (MQ=255)
                             tCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc                                            <  1:389656/90‑1 (MQ=255)
                             tCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc                                            <  1:190378/90‑1 (MQ=255)
                              cAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACg                                           >  2:237124/1‑90 (MQ=255)
                                gCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc                                            >  2:315318/1‑87 (MQ=255)
                                gCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc                                            <  1:315318/87‑1 (MQ=255)
                                gCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc                                         >  2:333899/1‑90 (MQ=255)
                                   ccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                     >  1:449144/1‑59 (MQ=255)
                                   ccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                     <  2:449144/59‑1 (MQ=255)
                                     tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa                                    <  1:79895/90‑1 (MQ=255)
                                     tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa                                    <  2:167193/90‑1 (MQ=255)
                                          ttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTc                               >  2:229234/1‑90 (MQ=255)
                                           tGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTcc                              <  2:128174/90‑1 (MQ=255)
                                                      gcgcAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGc                   >  1:429317/1‑90 (MQ=255)
                                                                       cTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAAc  >  1:364647/1‑90 (MQ=255)
                                                                                   |                                                                             
TACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAAC  >  NZ_CP009273/2859915‑2860075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCC  >  NZ_CP009273/2859909‑2860089
                                                                                         |                                                                                           
GATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACA                                                                                   >  SRR3722113.475862/1‑100 (MQ=60)
     TTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTC                                                                              >  SRR3722113.178343/1‑100 (MQ=60)
          TACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGG                                                                         <  SRR3722113.138439/100‑1 (MQ=60)
          TACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGG                                                                         <  SRR3722113.31218/100‑1 (MQ=60)
          TACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGG                                                                         <  SRR3722113.62840/100‑1 (MQ=60)
                 GGAGCAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCC                                                                  <  SRR3722113.383229/100‑1 (MQ=60)
                  GAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCA                                                                 >  SRR3722113.63110/1‑100 (MQ=60)
                               ATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGctgtctcttatacacatctcc                                                    >  SRR3722113.453896/1‑79 (MQ=60)
                                   TCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATC                                                <  SRR3722113.191866/100‑1 (MQ=60)
                                   TCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATC                                                <  SRR3722113.318403/100‑1 (MQ=60)
                                   TCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATC                                                <  SRR3722113.393737/100‑1 (MQ=60)
                                           TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTG                                        <  SRR3722113.80600/100‑1 (MQ=60)
                                                  GCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGC                                 >  SRR3722113.433860/1‑100 (MQ=60)
                                                                   CGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAAC                >  SRR3722113.368405/1‑100 (MQ=60)
                                                                                 AACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCC  >  SRR3722113.261030/1‑100 (MQ=60)
                                                                                         |                                                                                           
GATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGACGTGCC  >  NZ_CP009273/2859909‑2860089

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: