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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

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% 141.3 / NA 40D18G (GAT→GGT) yibBprotein YibB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/28);  total (12/28)

TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA  >  NZ_CP009273/3786904‑3787073
                                                                                       |                                                                                  
tCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAAccc                                                                                  <  1:685758/90‑1 (MQ=255)
  aaGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACcccc                                                                                <  2:523959/90‑1 (MQ=255)
    gTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctc                                                                              <  1:156644/90‑1 (MQ=255)
     tGCGGCAAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCctct                                                                             >  2:413122/1‑90 (MQ=255)
      gCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTa                                                                            >  1:253295/1‑90 (MQ=255)
       cGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAc                                                                           >  1:678118/1‑90 (MQ=255)
        ggCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAcc                                                                          >  2:118363/1‑90 (MQ=255)
        ggCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAcc                                                                          <  1:549253/90‑1 (MQ=255)
             gACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGt                                                                     <  1:222674/90‑1 (MQ=255)
                       gTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGc                                                           <  1:288440/90‑1 (MQ=255)
                       gTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGc                                                           <  2:642508/90‑1 (MQ=255)
                       gTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACAATGTCAAAATATGCa                                                          <  2:462809/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  2:532209/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  2:253295/90‑1 (MQ=255)
                        tAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCa                                                          <  1:36199/90‑1 (MQ=255)
                         aaCTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAg                                                         >  2:526250/1‑90 (MQ=255)
                           cTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       <  2:678118/90‑1 (MQ=255)
                           cTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       <  2:430903/90‑1 (MQ=255)
                                     cATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTgtgg                                             <  1:301537/90‑1 (MQ=255)
                                         aaCTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       >  2:325020/1‑76 (MQ=255)
                                         aaCTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTg                                                       <  1:325020/76‑1 (MQ=255)
                                           cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  2:609121/90‑1 (MQ=255)
                                           cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  2:748553/90‑1 (MQ=255)
                                           cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  1:61624/90‑1 (MQ=255)
                                           cTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGgatg                                       <  1:451795/90‑1 (MQ=255)
                                                cGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTc                                  <  1:118363/90‑1 (MQ=255)
                                                  aGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCAt                                <  1:206517/90‑1 (MQ=255)
                                                    cTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTa                              <  1:68167/90‑1 (MQ=255)
                                                       gTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACta                           >  1:592351/1‑90 (MQ=255)
                                                          tttCACGGACCCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                        <  2:174479/90‑1 (MQ=255)
                                                          tttCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                        <  1:348074/90‑1 (MQ=255)
                                                          tttCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACtaata                        <  1:322529/90‑1 (MQ=255)
                                                            tCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAg                      <  1:242552/90‑1 (MQ=255)
                                                             cACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTCATAAGa                     <  2:413395/90‑1 (MQ=255)
                                                             cACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGa                     <  1:675495/90‑1 (MQ=255)
                                                                ggAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCa                  >  1:115140/1‑90 (MQ=255)
                                                                               cAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAg   >  1:602394/1‑90 (MQ=255)
                                                                               cAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAg   >  2:656039/1‑90 (MQ=255)
                                                                                aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAatat           <  1:260908/81‑1 (MQ=255)
                                                                                aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAatat           >  2:260908/1‑81 (MQ=255)
                                                                                aGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGa  >  2:527262/1‑90 (MQ=255)
                                                                                       |                                                                                  
TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGA  >  NZ_CP009273/3786904‑3787073

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ 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
ATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATTATATCATA  >  NZ_CP009273/3786900‑3787084
                                                                                           |                                                                                             
ATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTA                                                                                       >  SRR3722073.258617/1‑100 (MQ=60)
 TTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTAC                                                                                      >  SRR3722073.691976/1‑100 (MQ=60)
    TCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAA                                                                                   <  SRR3722073.699730/100‑1 (MQ=60)
        GTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTC                                                                               <  SRR3722073.159815/100‑1 (MQ=60)
            GGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAA                                                                           <  SRR3722073.560721/100‑1 (MQ=60)
                 GACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGC                                                                      <  SRR3722073.227293/100‑1 (MQ=60)
                           GTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATT                                                            <  SRR3722073.294543/100‑1 (MQ=60)
                            TAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTG                                                           <  SRR3722073.37092/100‑1 (MQ=60)
                               acagAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGG                                                        <  SRR3722073.331924/96‑1 (MQ=60)
                                         CATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTT                                              <  SRR3722073.307912/100‑1 (MQ=60)
                                              ACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATnnnnnnnnnnnnnnnnnnnn                                         >  SRR3722073.244851/1‑80 (MQ=60)
                                               CTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAC                                        <  SRR3722073.63069/100‑1 (MQ=60)
                                               CTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTAC                                        <  SRR3722073.461452/100‑1 (MQ=60)
                                                 GAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTA                                      >  SRR3722073.604718/1‑100 (MQ=60)
                                                    CGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATA                                   <  SRR3722073.120879/100‑1 (MQ=60)
                                                      AGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAN                                 <  SRR3722073.210733/100‑1 (MQ=60)
                                                        CTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGAC                               <  SRR3722073.69763/100‑1 (MQ=60)
                                                          AGTTTTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCA                             >  SRR3722073.117610/1‑100 (MQ=60)
                                                              TTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAA                         <  SRR3722073.355445/100‑1 (MQ=60)
                                                              TTTCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAA                         <  SRR3722073.329376/100‑1 (MQ=60)
                                                                TCACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATA                       <  SRR3722073.247640/100‑1 (MQ=60)
                                                                 CACGGAACCCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATAT                      <  SRR3722073.689308/100‑1 (MQ=60)
                                                                         CCTTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAG              >  SRR3722073.614951/1‑100 (MQ=60)
                                                                           TTTATTAGCAGTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAA            <  SRR3722073.266409/100‑1 (MQ=60)
                                                                                     GTCCAACCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATActgtctcttatacacatctgacgctgccgacga  >  SRR3722073.705740/1‑67 (MQ=60)
                                                                                           |                                                                                             
ATTTTCAAGTGCGGCCAGACGCTCAAAGTAACTAAAATAAACATCAACTGAACGAGCTAGTTTTTCACGGAACCCTTTATTAGCAGTCCAATCCCCTCTACCAATGTCAAAATATGCAGTGATAATTGTGGTGGATGATTTCATTACTAATAAGACCATTAATATAAAAAAAGAATTATATCATA  >  NZ_CP009273/3786900‑3787084

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