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,635,419 T→C D262G (GAT→GGT)  rodZ ← cytoskeleton protein RodZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,635,4190TC100.0% 72.5 / NA 24D262G (GAT→GGT) rodZcytoskeleton protein RodZ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/11);  total (13/11)

TGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCAC  >  NZ_CP009273/2635332‑2635503
                                                                                       |                                                                                    
tGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCg                                                                                    >  1:197687/1‑90 (MQ=255)
       ttAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTaa                                                                             <  1:168665/90‑1 (MQ=255)
        tAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTaaa                                                                            <  2:123968/90‑1 (MQ=255)
                aaTTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTtcatca                                                                    >  2:120933/1‑90 (MQ=255)
                          tttACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCAtt                                                          <  1:326387/90‑1 (MQ=255)
                           ttACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTc                                                         <  1:191243/90‑1 (MQ=255)
                             aCGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCgg                                                       <  2:109255/90‑1 (MQ=255)
                                  gCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAg                                                  >  2:282612/1‑90 (MQ=255)
                                  gCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAg                                                  >  1:45834/1‑90 (MQ=255)
                                         gCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTtc                                                                        <  2:453247/61‑1 (MQ=255)
                                         gCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTtc                                                                        >  1:453247/1‑61 (MQ=255)
                                            aaaCAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGt                                        >  1:387272/1‑90 (MQ=255)
                                            aaaCAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGt                                        >  1:116534/1‑90 (MQ=255)
                                            aaaCAATTTTTTACCGGTAGCATCAGCGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGt                                        >  1:237087/1‑90 (MQ=255)
                                                 aTTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCa                                   >  1:449349/1‑90 (MQ=255)
                                                 aTTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCa                                   >  2:169434/1‑90 (MQ=255)
                                                                cATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCgg                    <  2:279495/90‑1 (MQ=255)
                                                                        aCCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAAc            <  1:433728/90‑1 (MQ=255)
                                                                        aCCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAAc            <  1:282612/90‑1 (MQ=255)
                                                                        aCCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAAc            <  1:209247/90‑1 (MQ=255)
                                                                            ccagccagCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGgcg        >  1:194282/1‑90 (MQ=255)
                                                                               gccagcAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCgcag     >  2:79597/1‑90 (MQ=255)
                                                                                 cagcAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCgcagca   >  1:382981/1‑90 (MQ=255)
                                                                                  agcAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCAc  <  2:116534/90‑1 (MQ=255)
                                                                                       |                                                                                    
TGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCAC  >  NZ_CP009273/2635332‑2635503

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ 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
GTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATC  >  NZ_CP009273/2635322‑2635507
                                                                                                 |                                                                                        
GTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCG                                                                                        >  SRR3722113.199220/1‑100 (MQ=60)
         CTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAG                                                                               <  SRR3722113.323976/100‑1 (MQ=60)
                 TTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCAC                                                                       <  SRR3722113.170029/100‑1 (MQ=60)
                                  TCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAG                                                      >  SRR3722113.46218/1‑100 (MQ=60)
                                    TTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCC                                                    <  SRR3722113.329641/100‑1 (MQ=60)
                                     TTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCA                                                   <  SRR3722113.192736/100‑1 (MQ=60)
                                         GCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCtgtctcttatacacatct                                               >  SRR3722113.458036/1‑82 (MQ=60)
                                            GCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGT                                            >  SRR3722113.117539/1‑100 (MQ=60)
                                            GCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGT                                            >  SRR3722113.391325/1‑100 (MQ=60)
                                            GCATACCGCTAAACAATTTTTTACCGGTAGCATCAGCGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGT                                            >  SRR3722113.239107/1‑100 (MQ=60)
                                                 CCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCA                                       >  SRR3722113.454102/1‑100 (MQ=60)
                                                                            TCAGTGACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCG            >  SRR3722113.195794/1‑100 (MQ=60)
                                                                                 GACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCA       >  SRR3722113.386983/1‑100 (MQ=60)
                                                                                  ACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCAC      <  SRR3722113.210848/100‑1 (MQ=60)
                                                                                  ACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCAC      <  SRR3722113.285156/100‑1 (MQ=60)
                                                                                  ACCTCCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCAC      <  SRR3722113.438327/100‑1 (MQ=60)
                                                                                      CCAGCCAGCAACCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATC  >  SRR3722113.44742/1‑100 (MQ=60)
                                                                                                 |                                                                                        
GTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATC  >  NZ_CP009273/2635322‑2635507

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