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,131,004 A→G intergenic (‑280/‑379) wza ← / → yegH polysaccharide export protein Wza/TerC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,131,0040AG89.7% 71.4 / 3.0 29intergenic (‑280/‑379)wza/yegHpolysaccharide export protein Wza/TerC family protein
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base G (11/15);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

TGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACG  >  NZ_CP009273/2130917‑2131087
                                                                                       |                                                                                   
tGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGgaa                                                                                   >  1:809085/1‑90 (MQ=255)
    tttttAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTaa                                                                               <  1:314914/90‑1 (MQ=255)
               cAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAatat                                                                    <  1:517772/90‑1 (MQ=255)
                    cacTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGt                                                               <  2:830810/90‑1 (MQ=255)
                    cacTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGt                                                               <  2:307822/90‑1 (MQ=255)
                    cacTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGt                                                               <  2:289304/90‑1 (MQ=255)
                    cacTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGt                                                               <  1:709925/90‑1 (MQ=255)
                    cacTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGt                                                               <  1:120098/90‑1 (MQ=255)
                     acTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTg                                                              <  1:397340/90‑1 (MQ=255)
                        gATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACa                                                           >  2:281034/1‑90 (MQ=255)
                               tttAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGt                                                    <  1:709466/90‑1 (MQ=255)
                                ttAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTa                                                   >  1:623933/1‑90 (MQ=255)
                                ttAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTa                                                   <  2:553120/90‑1 (MQ=255)
                                  aaTTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAAt                                                 >  1:387566/1‑90 (MQ=255)
                                      tttCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGaa                                             >  1:402330/1‑90 (MQ=255)
                                           cacCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTg                                        <  1:435775/90‑1 (MQ=255)
                                                tctcAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCttt                                   <  1:657660/90‑1 (MQ=255)
                                                  tcAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTcc                                 <  2:272394/90‑1 (MQ=255)
                                                     gCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCaat                              <  1:310767/90‑1 (MQ=255)
                                                     gCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCgat                              <  1:512455/90‑1 (MQ=255)
                                                         gCAGTCGTTGATGAGGAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATg                          >  1:196221/1‑90 (MQ=255)
                                                         gCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATg                          >  2:380766/1‑90 (MQ=255)
                                                          cAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGc                         >  2:534048/1‑90 (MQ=255)
                                                             tCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCagg                      >  2:758806/1‑90 (MQ=255)
                                                               gTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCaggag                    <  2:551398/90‑1 (MQ=255)
                                                                 tgatgaGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCaggagca                  >  2:495975/1‑88 (MQ=255)
                                                                 tgatgaGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAggagga                  >  2:791212/1‑90 (MQ=255)
                                                                   atgaGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAg                <  1:534048/90‑1 (MQ=255)
                                                                                attaCGGGAATTAACTTCCGAATATAAGGTGACATTGTGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAAc   <  1:319267/90‑1 (MQ=255)
                                                                                 ttaCGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACg  >  1:178193/1‑90 (MQ=255)
                                                                                       |                                                                                   
TGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACG  >  NZ_CP009273/2130917‑2131087

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ 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
CTACCGCCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGA  >  NZ_CP009273/2130907‑2131102
                                                                                                 |                                                                                                  
CTACCGCCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAA                                                                                                  >  SRR3722079.817767/1‑100 (MQ=60)
              TTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATA                                                                                    <  SRR3722079.317459/100‑1 (MQ=60)
                         CAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACAT                                                                         <  SRR3722079.522536/100‑1 (MQ=60)
                              CACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGG                                                                    <  SRR3722079.121084/100‑1 (MQ=60)
                              CACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGG                                                                    <  SRR3722079.717464/100‑1 (MQ=60)
                               ACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGT                                                                   <  SRR3722079.400670/100‑1 (MQ=60)
                                CTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTA                                                                  >  SRR3722079.630293/1‑100 (MQ=60)
                                  GATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAAT                                                                >  SRR3722079.390762/1‑100 (MQ=60)
                                      TAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAA                                                            >  SRR3722079.405731/1‑100 (MQ=60)
                                         TTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATAT                                                         <  SRR3722079.717000/100‑1 (MQ=60)
                                                     CACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAAT                                             <  SRR3722079.439557/100‑1 (MQ=60)
                                                         CTCTCAGCATGCAGTCGTTGATGAGGAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATG                                         >  SRR3722079.197912/1‑100 (MQ=60)
                                                          TCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGC                                        <  SRR3722079.664483/100‑1 (MQ=60)
                                                               GCATGCAGTCGTTGATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAG                                   <  SRR3722079.313294/100‑1 (MQ=60)
                                                               GCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCGATAATGCAGGAG                                   <  SRR3722079.517131/100‑1 (MQ=60)
                                                                             ATGAGAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCT                     <  SRR3722079.539048/100‑1 (MQ=60)
                                                                                 GAAAGGGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACG                 >  SRR3722079.179748/1‑100 (MQ=60)
                                                                                      GGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATA            >  SRR3722079.503316/1‑100 (MQ=60)
                                                                                      GGTTATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATA            >  SRR3722079.437473/1‑100 (MQ=60)
                                                                                          ATTACGGGAATTAACTTCCGAATATAAGGTGACATTGTGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAG        <  SRR3722079.321835/100‑1 (MQ=60)
                                                                                          ATTACGGGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAG        >  SRR3722079.608287/1‑100 (MQ=60)
                                                                                                GGAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGA  >  SRR3722079.482653/1‑100 (MQ=60)
                                                                                                 |                                                                                                  
CTACCGCCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGA  >  NZ_CP009273/2130907‑2131102

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