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 I209 R1 216 14.4 815930 96.1% 784108 85.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,494,539 T→C F222S (TTC→TCC)  nanX → MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,494,5390TC100.0% 63.7 / NA 20F222S (TTC→TCC) nanXMFS transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/11);  total (9/11)

AGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA  >  NZ_CP009273/4494476‑4494624
                                                               |                                                                                     
aGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTg                                                             >  2:374897/1‑90 (MQ=255)
   aTAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCaa                                                          <  2:50100/90‑1 (MQ=255)
     aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                        >  1:113098/1‑90 (MQ=255)
     aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                        >  2:394562/1‑90 (MQ=255)
     aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                        >  2:149944/1‑90 (MQ=255)
                 ttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAAc                                            <  1:265751/90‑1 (MQ=255)
                     ctgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                        <  1:149944/90‑1 (MQ=255)
                         ctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGactact                                    <  2:93948/90‑1 (MQ=255)
                                 aaCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGttt                                                                           >  1:212961/1‑43 (MQ=255)
                                 aaCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGttt                                                                           <  2:212961/43‑1 (MQ=255)
                                 aaCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTtccttcct                            <  1:187033/90‑1 (MQ=255)
                                         cTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCa                    <  1:61604/90‑1 (MQ=255)
                                              aaTCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATaa               <  1:12929/90‑1 (MQ=255)
                                                tctcTTTGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  2:21488/90‑1 (MQ=255)
                                                tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  1:29701/90‑1 (MQ=255)
                                                   ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg          >  2:356155/1‑90 (MQ=255)
                                                        aTCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAAt     <  1:384112/90‑1 (MQ=255)
                                                           gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa  >  1:76729/1‑90 (MQ=255)
                                                           gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa  >  2:403185/1‑90 (MQ=255)
                                                           gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa  >  1:326096/1‑90 (MQ=255)
                                                               |                                                                                     
AGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA  >  NZ_CP009273/4494476‑4494624

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ 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
CAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCA  >  NZ_CP009273/4494443‑4494631
                                                                                                |                                                                                            
CAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCT                                                                                           <  SRR3722091.126205/100‑1 (MQ=60)
CAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCT                                                                                           <  SRR3722091.158527/100‑1 (MQ=60)
                            ATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAC                                                               >  SRR3722091.114818/1‑100 (MQ=60)
                                                  TTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTC                                         <  SRR3722091.270055/100‑1 (MQ=60)
                                                      CTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTC                                     <  SRR3722091.152153/100‑1 (MQ=60)
                                                        GTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTctgtctcttatacacatctgacgctgccgacgaccat                                   >  SRR3722091.216232/1‑63 (MQ=60)
                                                                  AACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGA                         <  SRR3722091.189770/100‑1 (MQ=60)
                                                                          taagagacagCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAG                 <  SRR3722091.396723/90‑1 (MQ=60)
                                                                          CTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAG                 <  SRR3722091.62528/100‑1 (MQ=60)
                                                                           TTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCctgtctcttatacacatctgacgctgccgacgaccatcttagtgtagatctcggtggtcgc                >  SRR3722091.324828/1‑39 (MQ=60)
                                                                               AATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAAT            <  SRR3722091.13117/100‑1 (MQ=60)
                                                                                 TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC          <  SRR3722091.30113/100‑1 (MQ=60)
                                                                                  CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA         >  SRR3722091.331742/1‑100 (MQ=60)
                                                                                  CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA         >  SRR3722091.77919/1‑100 (MQ=60)
                                                                                         ATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCA  <  SRR3722091.390889/100‑1 (MQ=60)
                                                                                                |                                                                                            
CAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCA  >  NZ_CP009273/4494443‑4494631

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: