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 F4 I96 R1 227 23.3 882352 98.0% 864704 126.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 737,899 A→G N38N (AAT→AAC dtpD ← dipeptide permease DtpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273737,8990AG93.8% 44.9 / ‑3.2 16N38N (AAT→AACdtpDdipeptide permease DtpD
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (7/8);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

CCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGC  >  NZ_CP009273/737781‑738033
                                                                                                                      |                                                                                                                                      
ccATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATGGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATcag                                                                                                                    <  1:188933/139‑1 (MQ=255)
   tGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATcagcag                                                                                                                 >  2:110923/1‑139 (MQ=255)
                                                    gACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                                                <  2:83819/139‑1 (MQ=255)
                                                     aCATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                                                <  2:425883/138‑1 (MQ=255)
                                                      cATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGa                                                               >  2:286831/1‑138 (MQ=255)
                                                                  aaCAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgc                                                  >  1:415526/1‑139 (MQ=255)
                                                                     aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCAt                                                                                                        <  2:104673/82‑1 (MQ=255)
                                                                     aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCAt                                                                                                        >  1:104673/1‑82 (MQ=255)
                                                                     aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt                                                            >  2:182252/1‑126 (MQ=255)
                                                                     aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt                                                            <  1:182252/126‑1 (MQ=255)
                                                                     aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                               <  1:61994/139‑1 (MQ=255)
                                                                         ggCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGc                                           <  1:60977/139‑1 (MQ=255)
                                                                               aaaTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGAt                                     >  2:150443/1‑139 (MQ=255)
                                                                                     cTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtgt                               <  1:110923/139‑1 (MQ=255)
                                                                                          aGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAtt                          >  2:215631/1‑139 (MQ=255)
                                                                                                 aGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCataat                    >  2:224860/1‑138 (MQ=255)
                                                                                                                   ttggttggTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc  >  1:162203/1‑138 (MQ=255)
                                                                                                                      |                                                                                                                                      
CCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGC  >  NZ_CP009273/737781‑738033

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGT  >  NZ_CP009273/737780‑738049
                                                                                                                       |                                                                                                                                                      
GCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATGGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCA                                                                                                                         <  SRR3722119.190916/151‑1 (MQ=60)
                                                        cagCACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCG                                                                 <  SRR3722119.184141/148‑1 (MQ=60)
                                                        ATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCG                                                                 >  SRR3722119.419712/1‑151 (MQ=60)
                                                           CACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCctgtctcttatacacatctccgagcccacgagactaaggcgaatctc                                                              >  SRR3722119.105697/1‑104 (MQ=60)
                                                                     CAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATG                                                    <  SRR3722119.62569/151‑1 (MQ=60)
                                                                         AGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTG                                                <  SRR3722119.61545/151‑1 (MQ=60)
                                                                                     ACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAAT                                    <  SRR3722119.112011/151‑1 (MQ=60)
                                                                                                         TTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCA                 >  SRR3722119.163854/1‑150 (MQ=60)
                                                                                                                       GTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGActgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatg  >  SRR3722119.64840/1‑99 (MQ=60)
                                                                                                                       |                                                                                                                                                      
GCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAGT  >  NZ_CP009273/737780‑738049

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