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 3,786,991 T→C D18G (GAT→GGT)  yibB ← protein YibB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,786,9910TC100.0% 72.6 / NA 23D18G (GAT→GGT) yibBprotein YibB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/17);  total (6/17)

AAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAA  >  NZ_CP009273/3786906‑3787070
                                                                                     |                                                                               
aaGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                             <  1:601818/90‑1 (MQ=255)
aaGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                             <  2:251412/90‑1 (MQ=255)
  gTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctc                                                                           <  1:558373/90‑1 (MQ=255)
      ggCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAcc                                                                       <  2:516772/90‑1 (MQ=255)
        ccAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCaa                                                                     <  1:270528/90‑1 (MQ=255)
                      tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                       <  2:29609/90‑1 (MQ=255)
                      tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                       <  2:294846/90‑1 (MQ=255)
                      tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                       <  1:467838/90‑1 (MQ=255)
                      tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                       <  1:475189/90‑1 (MQ=255)
                      tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                       <  1:603041/90‑1 (MQ=255)
                           aaaaTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGAt                                                  >  1:439725/1‑90 (MQ=255)
                              aTAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAAt                                               <  2:155154/90‑1 (MQ=255)
                                       aaCTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                    <  1:156748/76‑1 (MQ=255)
                                       aaCTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                    >  2:156748/1‑76 (MQ=255)
                                       aaCTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGga                                      >  1:207354/1‑90 (MQ=255)
                                         cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                    <  2:185126/90‑1 (MQ=255)
                                         cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                    <  1:37624/90‑1 (MQ=255)
                                          tGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatga                                   >  2:411730/1‑90 (MQ=255)
                                             aCGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGAttt                                >  1:11763/1‑90 (MQ=255)
                                                  cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                           <  1:75642/90‑1 (MQ=255)
                                                                           aGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATaaaaaa  <  2:275640/90‑1 (MQ=255)
                                                                           aGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATaaaaaa  >  2:399336/1‑90 (MQ=255)
                                                                           aGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATaaaaaa  <  2:408146/90‑1 (MQ=255)
                                                                                     |                                                                               
AAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAA  >  NZ_CP009273/3786906‑3787070

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CATCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGAC  >  NZ_CP009273/3786894‑3787055
                                                                                                 |                                                                
CATCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCC                                                                <  SRR3722088.127760/100‑1 (MQ=60)
CATCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCC                                                                <  SRR3722088.323813/100‑1 (MQ=60)
    TCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTC                                                            <  SRR3722088.303757/100‑1 (MQ=60)
       TTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAC                                                         <  SRR3722088.560289/100‑1 (MQ=60)
        TTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACC                                                        <  SRR3722088.419106/100‑1 (MQ=60)
         TTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCA                                                       <  SRR3722088.629070/100‑1 (MQ=60)
         TTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCA                                                       <  SRR3722088.478461/100‑1 (MQ=60)
            AAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATG                                                    <  SRR3722088.609229/100‑1 (MQ=60)
              GTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTC                                                  <  SRR3722088.565223/100‑1 (MQ=60)
                    CCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATA                                            <  SRR3722088.273142/100‑1 (MQ=60)
                             CAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGAT                                   >  SRR3722088.444780/1‑100 (MQ=60)
                                  TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                              <  SRR3722088.473332/100‑1 (MQ=60)
                                  TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                              <  SRR3722088.480832/100‑1 (MQ=60)
                                  TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                              <  SRR3722088.610466/100‑1 (MQ=60)
                                     acagAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGG                           <  SRR3722088.158324/96‑1 (MQ=60)
                                         AATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGA                       >  SRR3722088.209424/1‑100 (MQ=60)
                                               CATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTT                 >  SRR3722088.11896/1‑100 (MQ=60)
                                                     CTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAC           <  SRR3722088.38020/100‑1 (MQ=60)
                                                              CTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGAC  <  SRR3722088.76409/100‑1 (MQ=60)
                                                                                                 |                                                                
CATCTCATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGAC  >  NZ_CP009273/3786894‑3787055

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: