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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,726,562 (T)5→4 coding (971/2574 nt) clpB ← ATP‑dependent chaperone ClpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,726,5580T.100.0% 101.3 / NA 25coding (975/2574 nt)clpBATP‑dependent chaperone ClpB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (11/14);  total (11/14)

ACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCA  >  NZ_CP009273/2726480‑2726640
                                                                              |                                                                                  
aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACGCTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                        <  1:129307/90‑1 (MQ=255)
aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                        <  1:147710/90‑1 (MQ=255)
aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                        <  2:66790/90‑1 (MQ=255)
aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                        <  1:460997/90‑1 (MQ=255)
        tCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGata                                                                <  1:390044/90‑1 (MQ=255)
          gCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATAtt                                                              >  2:158477/1‑90 (MQ=255)
               ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCa                                                         >  1:66722/1‑90 (MQ=255)
                        aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                <  2:794818/90‑1 (MQ=255)
                        aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                <  1:274237/90‑1 (MQ=255)
                         aCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg                                               >  1:266682/1‑90 (MQ=255)
                              aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                          >  1:179633/1‑90 (MQ=255)
                              aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                          >  1:24777/1‑90 (MQ=255)
                              aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                          >  2:154539/1‑90 (MQ=255)
                                gCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTAc                                        <  2:55723/90‑1 (MQ=255)
                                  tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                      <  1:56748/90‑1 (MQ=255)
                                  tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                      <  1:653062/90‑1 (MQ=255)
                                    ggCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAg                                    <  1:154539/90‑1 (MQ=255)
                                          aaacaCTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTAc                                        >  2:341769/1‑80 (MQ=255)
                                          aaacaCTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTAc                                        <  1:341769/80‑1 (MQ=255)
                                            acacTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTc                            <  1:158477/90‑1 (MQ=255)
                                                   cTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCAcgc                     >  1:351546/1‑90 (MQ=255)
                                                         acgacgTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAgc               >  1:802304/1‑90 (MQ=255)
                                                            acgTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGcc            <  2:838411/90‑1 (MQ=255)
                                                               ttCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGt         >  2:266806/1‑90 (MQ=255)
                                                                      gcagcaTCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCa  >  1:69792/1‑90 (MQ=255)
                                                                              |                                                                                  
ACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCA  >  NZ_CP009273/2726480‑2726640

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATT  >  NZ_CP009273/2726463‑2726657
                                                                                                   |                                                                                               
TAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTC                                                                                                 <  SRR3722072.663735/100‑1 (MQ=60)
                 ACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACGCTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATT                                                                               <  SRR3722072.131643/100‑1 (MQ=60)
                 ACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATT                                                                               <  SRR3722072.150330/100‑1 (MQ=60)
                 ACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATT                                                                               <  SRR3722072.469439/100‑1 (MQ=60)
                      GAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCA                                                                          >  SRR3722072.67969/1‑100 (MQ=60)
                         TCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGC                                                                       <  SRR3722072.397041/100‑1 (MQ=60)
                                GGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGG                                                                >  SRR3722072.271226/1‑100 (MQ=60)
                                     CTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCT                                                           >  SRR3722072.182743/1‑100 (MQ=60)
                                     CTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCT                                                           >  SRR3722072.25214/1‑100 (MQ=60)
                                         AACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGC                                                       <  SRR3722072.278927/100‑1 (MQ=60)
                                                 GCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAAT                                               <  SRR3722072.347793/100‑1 (MQ=60)
                                                   TCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTC                                             <  SRR3722072.57785/100‑1 (MQ=60)
                                                   TCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTC                                             <  SRR3722072.664719/100‑1 (MQ=60)
                                                     GGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCAC                                           <  SRR3722072.157260/100‑1 (MQ=60)
                                                          CAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGC                                      >  SRR3722072.357743/1‑100 (MQ=60)
                                                             ACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCC                                   <  SRR3722072.161243/100‑1 (MQ=60)
                                                                CTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGC                                >  SRR3722072.816317/1‑100 (MQ=60)
                                                                             ACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCA                   >  SRR3722072.71092/1‑100 (MQ=60)
                                                                                              tcTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCtgtctcttatacacatctccgagccc   >  SRR3722072.147375/3‑74 (MQ=60)
                                                                                               CTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATT  >  SRR3722072.490213/1‑100 (MQ=60)
                                                                                                   |                                                                                               
TAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCATGTTTCCGGCGTCCATT  >  NZ_CP009273/2726463‑2726657

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: