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 I225 R1 227 21.4 1179702 97.1% 1145490 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 674,442 T→G T141P (ACT→CCT)  ybeT ← tetratricopeptide repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273674,4420TG100.0% 104.2 / NA 32T141P (ACT→CCT) ybeTtetratricopeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (14/18);  total (14/18)

TGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGA  >  NZ_CP009273/674357‑674516
                                                                                     |                                                                          
tgttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGatat                                                                        <  2:223382/90‑1 (MQ=255)
 gttTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATAtt                                                                       >  1:34415/1‑90 (MQ=255)
      gCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagca                                                                  <  1:150945/90‑1 (MQ=255)
       cGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcat                                                                 <  1:453904/90‑1 (MQ=255)
             gTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTgagcatgagcat                                                           <  1:360323/90‑1 (MQ=255)
                       tttttGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcat                                                                 <  2:584560/74‑1 (MQ=255)
                       tttttGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATtgagcat                                                                 >  1:584560/1‑74 (MQ=255)
                       tttttGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTc                                                 >  2:270551/1‑90 (MQ=255)
                       tttttGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTc                                                 >  1:227989/1‑90 (MQ=255)
                             agaTTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgc                                           <  1:474099/90‑1 (MQ=255)
                              gaTTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgct                                          <  2:480013/90‑1 (MQ=255)
                              gaTTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgct                                          <  2:381693/90‑1 (MQ=255)
                              gaTTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgct                                          <  1:126330/90‑1 (MQ=255)
                                 tCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTcc                                       >  2:291983/1‑90 (MQ=255)
                                     aaTGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTtat                                   <  2:334823/90‑1 (MQ=255)
                                        ggTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGatat                                                                        >  1:304284/1‑50 (MQ=255)
                                        ggTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGatat                                                                        <  2:304284/50‑1 (MQ=255)
                                         gTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTgagcatgagcat                                                           <  1:580724/62‑1 (MQ=255)
                                         gTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTgagcatgagcat                                                           >  2:580724/1‑62 (MQ=255)
                                                ttAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAAtt                        <  2:64534/90‑1 (MQ=255)
                                                ttAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAAtt                        <  2:127070/90‑1 (MQ=255)
                                                ttAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAAtt                        <  1:270551/90‑1 (MQ=255)
                                                        tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgc                                           <  2:297702/63‑1 (MQ=255)
                                                        tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAgctgc                                           >  1:297702/1‑63 (MQ=255)
                                                        tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                >  2:287844/1‑90 (MQ=255)
                                                        tGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACg                >  2:111572/1‑90 (MQ=255)
                                                         gAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACgg               >  2:423983/1‑90 (MQ=255)
                                                         gAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACgg               <  1:111572/90‑1 (MQ=255)
                                                              aTTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTc          >  1:529227/1‑90 (MQ=255)
                                                                    gTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGttt    >  2:240393/1‑90 (MQ=255)
                                                                    gTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGttt    >  1:496642/1‑90 (MQ=255)
                                                                      aGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGa  <  2:389036/90‑1 (MQ=255)
                                                                                     |                                                                          
TGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTTGA  >  NZ_CP009273/674357‑674516

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTT  >  NZ_CP009273/674344‑674514
                                                                                                  |                                                                        
ATTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGA                                                                         <  SRR3722109.546934/100‑1 (MQ=60)
    TTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATT                                                                     >  SRR3722109.34818/1‑100 (MQ=60)
                   GCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACC                                                      <  SRR3722109.152810/100‑1 (MQ=60)
                    CGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCC                                                     <  SRR3722109.460716/100‑1 (MQ=60)
                          gatgtgtataagagacAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTC                                               <  SRR3722109.589517/84‑1 (MQ=60)
                          GTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCctgtct                                               >  SRR3722109.593402/1‑94 (MQ=60)
                          GTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTC                                               >  SRR3722109.230696/1‑100 (MQ=60)
                          GTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTC                                               <  SRR3722109.365298/100‑1 (MQ=60)
                                          AGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATAC                               <  SRR3722109.481238/100‑1 (MQ=60)
                                           GATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCtgtctcttatacacatctgacgctgccga                              >  SRR3722109.308238/1‑71 (MQ=60)
                                           GATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACC                              <  SRR3722109.127887/100‑1 (MQ=60)
                                                           CATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCtgtctcttatacacat              >  SRR3722109.301546/1‑84 (MQ=60)
                                                             TTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGG            <  SRR3722109.273915/100‑1 (MQ=60)
                                                                 TACCTGAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTC        >  SRR3722109.537281/1‑100 (MQ=60)
                                                                      GAGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTT   <  SRR3722109.112958/100‑1 (MQ=60)
                                                                       AGCCATTTCTGTAGATCAAACCTAAAGGATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTT  >  SRR3722109.504138/1‑100 (MQ=60)
                                                                                                  |                                                                        
ATTTTTGTAATGTTGTTTGGCGGTCAGTTTTAACCATTTTTGAGATTCATAATGGTTAACATTAATACCTGAGCCATTTCTGTAGATCAAACCTAAAGTATATTGAGCATGAGCATACCCTTGTTCAGCTGCTTCCTTATACCAGAGAATTGCTTTACGGGTGTCGGGTTT  >  NZ_CP009273/674344‑674514

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