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 I200 R1 216 30.8 1683622 97.2% 1636480 87.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,061,342 T→A intergenic (‑576/‑774) BW25113_RS10415 ← / → BW25113_RS24005 IS5‑like element IS5 family transposase/TonB‑dependent receptor plug domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,061,3420TA83.9% 82.3 / 10.2 31intergenic (‑576/‑774)BW25113_RS10415/BW25113_RS24005IS5‑like element IS5 family transposase/TonB‑dependent receptor plug domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base A (10/16);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.17e-02

AAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGT  >  NZ_CP009273/2061258‑2061426
                                                                                    |                                                                                    
aaGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTa                                                                                 >  1:567580/1‑90 (MQ=255)
  gCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAgt                                                                               <  1:739646/90‑1 (MQ=255)
             tctcAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGaa                                                                    <  2:398146/90‑1 (MQ=255)
             tctcAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGaa                                                                    <  1:153145/90‑1 (MQ=255)
                            tctcAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGc                                                     <  1:542165/90‑1 (MQ=255)
                                 gACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGt                                                <  1:153974/90‑1 (MQ=255)
                                 gACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGt                                                <  2:35052/90‑1 (MQ=255)
                                   cacaGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTAc                                              <  2:81311/90‑1 (MQ=255)
                                       gCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACt                                          >  1:557267/1‑90 (MQ=255)
                                         tttCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCg                                        >  2:332863/1‑90 (MQ=255)
                                            cTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGtca                                     >  1:605619/1‑90 (MQ=255)
                                               cTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCAc                                           >  2:530524/1‑81 (MQ=255)
                                               cTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCAc                                           <  1:530524/81‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                               <  1:677326/90‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                               <  2:173292/90‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                               <  2:35089/90‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                               <  2:539150/90‑1 (MQ=255)
                                                  cATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGt                               <  2:557267/90‑1 (MQ=255)
                                                   aTGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGtcatc                                   >  1:604665/1‑85 (MQ=255)
                                                   aTGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGtcatc                                   <  2:604665/85‑1 (MQ=255)
                                                   aTGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTa                              >  2:242832/1‑90 (MQ=255)
                                                         aTGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtata                        >  2:607842/1‑90 (MQ=255)
                                                         aTGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtata                        >  2:655424/1‑90 (MQ=255)
                                                           ggCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatata                      >  1:2095/1‑90 (MQ=255)
                                                           ggCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTtatata                      >  2:759523/1‑90 (MQ=255)
                                                                    aCGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGc             <  2:459384/90‑1 (MQ=255)
                                                                       aaGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAt          <  2:2653/90‑1 (MQ=255)
                                                                            gCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAttgttg     >  2:767286/1‑90 (MQ=255)
                                                                            gCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAttgttg     <  1:463919/90‑1 (MQ=255)
                                                                            gCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTAttgttg     <  2:836742/90‑1 (MQ=255)
                                                                               cAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTgtgt  <  2:673315/90‑1 (MQ=255)
                                                                                    |                                                                                    
AAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGT  >  NZ_CP009273/2061258‑2061426

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTTTGTACATCAGTAG  >  NZ_CP009273/2061248‑2061441
                                                                                              |                                                                                                   
CTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTA                                                                                                >  SRR3722079.573137/1‑100 (MQ=60)
 TCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAG                                                                                               <  SRR3722079.837293/100‑1 (MQ=60)
    CCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGG                                                                                            <  SRR3722079.443690/100‑1 (MQ=60)
            GCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGA                                                                                    <  SRR3722079.747553/100‑1 (MQ=60)
                       TCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAAC                                                                         <  SRR3722079.154474/100‑1 (MQ=60)
                                      TCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCAC                                                          <  SRR3722079.547305/100‑1 (MQ=60)
                                       CTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACT                                                         >  SRR3722079.562656/1‑100 (MQ=60)
                                           GACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTC                                                     <  SRR3722079.155310/100‑1 (MQ=60)
                                            ACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCA                                                    >  SRR3722079.611727/1‑100 (MQ=60)
                                                AGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCAC                                                <  SRR3722079.535472/100‑1 (MQ=60)
                                                   TTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTA                                             >  SRR3722079.610763/1‑100 (MQ=60)
                                                           CCATGGTAATGGCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATA                                     >  SRR3722079.2112/1‑100 (MQ=60)
                                                            CATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATAT                                    <  SRR3722079.684434/100‑1 (MQ=60)
                                                                      GCCCAAATACGAAGTTGCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCct                          >  SRR3722079.463821/1‑98 (MQ=60)
                                                                                      GCTCAGGAACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTTTGTAC          <  SRR3722079.467979/100‑1 (MQ=60)
                                                                                           GGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTTTGTACATCAG     >  SRR3722079.295378/1‑100 (MQ=60)
                                                                                              ACATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTTTGTACATCAGTAG  <  SRR3722079.775496/100‑1 (MQ=60)
                                                                                              |                                                                                                   
CTCTCCAAGAAAGCGAATTCTTTTCTCAAATTCTGCATTCTCAGACACAGCTTTCTACTCCATGGTAATGGCCCAAATACGAAGTTGCTCAGGATCATTAGTGGTGCTGAGAAGTTAAATAACCATGCGGAGTACCACTCGTCATCACGTATTTATATATCATGAAGCTATTGTTGTGTTTTGTACATCAGTAG  >  NZ_CP009273/2061248‑2061441

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