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,102,646 T→C G140G (GGA→GGG wzx ← O16 family O‑antigen flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,102,6460TC91.7% 60.8 / ‑0.7 24G140G (GGA→GGGwzxO16 family O‑antigen flippase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (11/11);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTC  >  NZ_CP009273/2102557‑2102732
                                                                                         |                                                                                      
gACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTttttctttt                                                                                        >  1:780984/1‑90 (MQ=255)
      gtgtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGa                                                                                  >  1:13688/1‑90 (MQ=255)
      gtgtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGa                                                                                  >  2:290329/1‑90 (MQ=255)
        gtTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGaaa                                                                                <  2:382436/90‑1 (MQ=255)
          tCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACa                                                                              >  1:601141/1‑90 (MQ=255)
            acaAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc                                                                            <  2:452012/90‑1 (MQ=255)
            acaAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc                                                                            <  2:601141/90‑1 (MQ=255)
                 aataaatGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGAt                                                                       <  2:410360/90‑1 (MQ=255)
                                         cAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTaa                                               >  1:729033/1‑90 (MQ=255)
                                          aaGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc                                                                            >  1:391251/1‑60 (MQ=255)
                                          aaGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGc                                                                            <  2:391251/60‑1 (MQ=255)
                                                gATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGa                                        >  1:411099/1‑90 (MQ=255)
                                                       tactacTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACg                                 <  2:729033/90‑1 (MQ=255)
                                                         ctactTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGAc                               >  1:582412/1‑90 (MQ=255)
                                                         ctactTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGAc                               <  1:11456/90‑1 (MQ=255)
                                                            ctTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCat                            <  1:172202/90‑1 (MQ=255)
                                                              tAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCatat                          <  1:547548/90‑1 (MQ=255)
                                                              tAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCatat                          >  2:114662/1‑90 (MQ=255)
                                                                     aGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACg                   >  1:80532/1‑90 (MQ=255)
                                                                     aGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACg                   >  2:70390/1‑90 (MQ=255)
                                                                      gCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACgg                  >  1:120646/1‑90 (MQ=255)
                                                                          tttcatttTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATg              <  1:153260/90‑1 (MQ=255)
                                                                            tcatttTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCt            <  2:745253/90‑1 (MQ=255)
                                                                            tcatttTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCt            <  2:512925/90‑1 (MQ=255)
                                                                               ttttttCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGt         <  2:15678/90‑1 (MQ=255)
                                                                                   ttctttccCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTtata     >  2:814075/1‑90 (MQ=255)
                                                                                      tttccCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTc  <  1:140425/90‑1 (MQ=255)
                                                                                         |                                                                                      
GACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTC  >  NZ_CP009273/2102557‑2102732

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ 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
TGCTATATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCAGGAACAAGCAAG  >  NZ_CP009273/2102547‑2102745
                                                                                                   |                                                                                                   
TGCTATATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTT                                                                                                     >  SRR3722079.789354/1‑100 (MQ=60)
TGCTATATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTC                                                                                                     >  SRR3722079.79675/1‑100 (MQ=60)
      ATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGA                                                                                               >  SRR3722079.13787/1‑100 (MQ=60)
          GACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACA                                                                                           >  SRR3722079.607193/1‑100 (MQ=60)
                                         GGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAA                                                            >  SRR3722079.736824/1‑100 (MQ=60)
                                          GATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTctgtctcttatacacatctc                                                           >  SRR3722079.394505/1‑80 (MQ=60)
                                                AGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGA                                                     >  SRR3722079.414614/1‑100 (MQ=60)
                                                         GGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGAC                                            >  SRR3722079.588178/1‑100 (MQ=60)
                                                                   CTACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTA                                  <  SRR3722079.11542/100‑1 (MQ=60)
                                                                     ACTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACG                                >  SRR3722079.81186/1‑100 (MQ=60)
                                                                      CTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGG                               >  SRR3722079.121639/1‑100 (MQ=60)
                                                                      CTTAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGG                               <  SRR3722079.173699/100‑1 (MQ=60)
                                                                        TAAAGTCAGCCATTTCATTTTTTCTTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCA                             <  SRR3722079.552781/100‑1 (MQ=60)
                                                                                    TTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATAT                 <  SRR3722079.154590/100‑1 (MQ=60)
                                                                                                TTTCCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCAGGAACAAGC     <  SRR3722079.141619/100‑1 (MQ=60)
                                                                                                   CCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCAGGAACAAGCAAG  >  SRR3722079.339626/1‑100 (MQ=60)
                                                                                                   |                                                                                                   
TGCTATATCTGACTTTGTGTTCACAAAAATAAATGTTAGAGGGATAATAGCCAAGCGGGATAAAATACTACTTAAAGTCAGCCATTTCATTTTTTCTTTTCCCTGAAACAGCCAGATAGGGTAGATTAAATTCCCGACTAATGCAGGAACAAACGACCATATAATTACGGCATGCTTGTTATATTCAGGAACAAGCAAG  >  NZ_CP009273/2102547‑2102745

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: