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 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% 88.7 / NA 28F222S (TTC→TCC) nanXMFS transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/15);  total (13/15)

TGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA  >  NZ_CP009273/4494457‑4494624
                                                                                  |                                                                                     
tGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTc                                                                                <  1:401618/90‑1 (MQ=255)
      gAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTc                                                                                <  1:136306/84‑1 (MQ=255)
      gAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTc                                                                                >  2:136306/1‑84 (MQ=255)
      gAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGtttt                                                                          >  1:433681/1‑90 (MQ=255)
              atatAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTAtt                                                                  <  2:163878/90‑1 (MQ=255)
                       tAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCaaa                                                         <  2:560689/90‑1 (MQ=255)
                        aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCACTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                        >  1:557084/1‑90 (MQ=255)
                            cAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGc                                                    <  1:501678/90‑1 (MQ=255)
                                    ttGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAAc                                            >  1:509566/1‑90 (MQ=255)
                                     tgtctgtcTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACg                                           <  2:379837/90‑1 (MQ=255)
                                        ctgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                        >  1:299792/1‑90 (MQ=255)
                                         tgtctTCAGAAACCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGact                                       <  1:8525/90‑1 (MQ=255)
                                         tgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGact                                       <  1:161424/90‑1 (MQ=255)
                                             ttCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActt                                   <  2:449973/90‑1 (MQ=255)
                                                    aaCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTtccttcct                            >  1:126888/1‑90 (MQ=255)
                                                    aaCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTtccttcct                            >  2:343308/1‑90 (MQ=255)
                                                        acaTCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCt                        <  2:237146/90‑1 (MQ=255)
                                                           tCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGc                     <  2:443118/90‑1 (MQ=255)
                                                               tCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGAt                 >  1:295154/1‑90 (MQ=255)
                                                                   tctcTATGATCGTTTCCCTCGTCTATTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  2:574321/90‑1 (MQ=255)
                                                                      ctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAg          >  1:14915/1‑90 (MQ=255)
                                                                       tATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGt         >  2:510418/1‑90 (MQ=255)
                                                                        aTGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGtt        <  2:488937/90‑1 (MQ=255)
                                                                        aTGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGtt        <  1:104079/90‑1 (MQ=255)
                                                                          gATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTaa      <  2:201996/90‑1 (MQ=255)
                                                                              gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa  >  2:135236/1‑90 (MQ=255)
                                                                              gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCATACCTGGCAGATAATGGAGTTAATACa  >  1:308903/1‑90 (MQ=255)
                                                                              gTTTCCCTCGTCTGTTTCTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa  >  1:160350/1‑90 (MQ=255)
                                                                                  |                                                                                     
TGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA  >  NZ_CP009273/4494457‑4494624

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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
GGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTC  >  NZ_CP009273/4494453‑4494635
                                                                                      |                                                                                                
GGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTT                                                                                     >  SRR3722109.440093/1‑100 (MQ=60)
    TGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTC                                                                                 <  SRR3722109.137994/100‑1 (MQ=60)
    TGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTC                                                                                 <  SRR3722109.407427/100‑1 (MQ=60)
                  ATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCACTCGTCTGTTTTTGTCTATTTGGTGCAAAC                                                                   >  SRR3722109.565536/1‑100 (MQ=60)
                              ATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAAC                                                       >  SRR3722109.517279/1‑100 (MQ=60)
                                CAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGG                                                     <  SRR3722109.509251/100‑1 (MQ=60)
                                  ACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGAC                                                   >  SRR3722109.303673/1‑100 (MQ=60)
                                             TGTCTTCAGAAACCCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCC                                        <  SRR3722109.8642/100‑1 (MQ=60)
                                             TGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCC                                        <  SRR3722109.163425/100‑1 (MQ=60)
                                              GTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCT                                       >  SRR3722109.128451/1‑100 (MQ=60)
                                                cgtcggcagcgtcagatgtgtataagagacagTCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTAC                                     <  SRR3722109.192963/68‑1 (MQ=60)
                                                         ACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGAT                            >  SRR3722109.298955/1‑100 (MQ=60)
                                                                CTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAG                     >  SRR3722109.15114/1‑100 (MQ=60)
                                                                        CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCATACCTGGCAGATAATGGAGTTAATACA             >  SRR3722109.312946/1‑100 (MQ=60)
                                                                        CTCTATGATCGTTTCCCTCGTCTGTTTCTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA             >  SRR3722109.162335/1‑100 (MQ=60)
                                                                          CTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGC           >  SRR3722109.468591/1‑100 (MQ=60)
                                                                            ATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGG         <  SRR3722109.105380/100‑1 (MQ=60)
                                                                                   TTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTC  >  SRR3722109.360008/1‑100 (MQ=60)
                                                                                      |                                                                                                
GGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTC  >  NZ_CP009273/4494453‑4494635

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: