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 F3 I96 R1 172 34.7 2005794 96.2% 1929573 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,632,437 (T)7→6 intergenic (‑71/+8) BW25113_RS08135 ← / ← BW25113_RS25620 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,632,4310T.100.0% 139.2 / NA 33intergenic (‑65/+14)BW25113_RS08135/BW25113_RS25620hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (13/20);  total (13/20)

ATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAGGAGCCTTCATTTTCA  >  NZ_CP009273/1632349‑1632515
                                                                                  |                                                                                    
aTTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttttttca                                                                               <  2:504823/90‑3 (MQ=255)
    ggATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCatgat                                                                          <  2:422684/90‑1 (MQ=255)
    ggATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCatgat                                                                          <  1:699890/90‑1 (MQ=255)
        atTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTa                                                                      >  2:901406/1‑90 (MQ=255)
          tCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTAcc                                                                    <  2:964365/90‑1 (MQ=255)
          tCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTAcc                                                                    <  2:355755/90‑1 (MQ=255)
            aCGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                  <  1:671878/90‑1 (MQ=255)
            aCGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                  <  2:435514/90‑1 (MQ=255)
            aCGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                  <  1:253414/90‑1 (MQ=255)
                atgatgGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgtagta                                                              <  2:525796/90‑1 (MQ=255)
                          tAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATc                                                    <  2:540237/90‑1 (MQ=255)
                          tAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATc                                                    <  1:882627/90‑1 (MQ=255)
                           aaCTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCatgat                                                                          >  1:770399/1‑67 (MQ=255)
                           aaCTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCatgat                                                                          <  2:770399/67‑1 (MQ=255)
                            aCTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCa                                                  <  2:676056/90‑1 (MQ=255)
                              tCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATc                                                <  1:649223/90‑1 (MQ=255)
                              tCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATc                                                <  1:777058/90‑1 (MQ=255)
                                cAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTAcc                                                                    >  1:192332/1‑68 (MQ=255)
                                cAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTAcc                                                                    <  2:192332/68‑1 (MQ=255)
                                cAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGa                                              >  2:264473/1‑90 (MQ=255)
                                cAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGa                                              >  1:926816/1‑90 (MQ=255)
                                  gTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAAt                                            >  2:414127/1‑90 (MQ=255)
                                                tACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                  >  2:23003/1‑54 (MQ=255)
                                                tACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                  <  1:23003/54‑1 (MQ=255)
                                                tACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCaaa                              >  1:121867/1‑90 (MQ=255)
                                                     ctgctgGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                  >  2:404868/1‑49 (MQ=39)
                                                     ctgctgGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                  <  1:404868/49‑1 (MQ=39)
                                                     ctgctgGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGc                         <  1:982058/90‑1 (MQ=255)
                                                       gctgGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTc                       >  2:720625/1‑90 (MQ=255)
                                                        ctgGTTCAGAGTTAAAAAAACTAAGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTcc                      <  2:655303/90‑1 (MQ=255)
                                                         tgGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTccc                     >  1:817217/1‑90 (MQ=255)
                                                                          aaCTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATc                                                    >  2:671985/1‑42 (MQ=38)
                                                                          aaCTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATc                                                    <  1:671985/42‑1 (MQ=38)
                                                                            cTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAGGAGCCTTCATTTTCa  >  1:927964/1‑90 (MQ=255)
                                                                                  |                                                                                    
ATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAGGAGCCTTCATTTTCA  >  NZ_CP009273/1632349‑1632515

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ 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
TTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAGGAGCCTTCATTTTCA  >  NZ_CP009273/1632351‑1632515
                                                                                      |                                                                              
agatgtgtataagagacagGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGT                                                                  <  SRR3722236.610049/81‑1 (MQ=60)
  GGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAG                                                                <  SRR3722236.711704/100‑1 (MQ=60)
          gtcggcagcgtcagatgtgtataagagacagGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAG                                                        <  SRR3722236.410535/69‑1 (MQ=60)
          cagcgtcagatgtgtataagagacagTCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAG                                                        <  SRR3722236.23345/74‑1 (MQ=60)
          ACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAG                                                        <  SRR3722236.257156/100‑1 (MQ=60)
          ACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAG                                                        <  SRR3722236.683066/100‑1 (MQ=60)
               TGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGctgtctcttatac                                                   >  SRR3722236.783650/1‑87 (MQ=60)
                    GTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGctgtctcttata                                              >  SRR3722236.195183/1‑88 (MQ=60)
                    GTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGA                                              >  SRR3722236.943184/1‑100 (MQ=60)
                        agccttcgtcggcagcgtcagatgtgtataagagacagGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAG                                          <  SRR3722236.683174/62‑1 (MQ=60)
                        TAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAG                                          <  SRR3722236.898148/100‑1 (MQ=60)
                            TCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACAC                                      <  SRR3722236.790450/100‑1 (MQ=60)
                            TCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACAC                                      <  SRR3722236.659938/100‑1 (MQ=60)
                                    TCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAA                              >  SRR3722236.123580/1‑100 (MQ=60)
                                             CTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCC                     >  SRR3722236.831458/1‑100 (MQ=60)
                                                   CTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAGGA               <  SRR3722236.999364/100‑1 (MQ=60)
                                                                GTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAGGAGCCTTCATTTTCA  >  SRR3722236.944353/1‑100 (MQ=60)
                                                                                      |                                                                              
TTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAGGAGCCTTCATTTTCA  >  NZ_CP009273/1632351‑1632515

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