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,425,778 T→C E215G (GAA→GGA)  folC ← bifunctional tetrahydrofolate synthase/dihydrofolate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,425,7780TC93.5% 85.8 / ‑0.5 31E215G (GAA→GGA) folCbifunctional tetrahydrofolate synthase/dihydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (20/9);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCT  >  NZ_CP009273/2425693‑2425858
                                                                                     |                                                                                
gCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTg                                                                              >  2:617604/1‑90 (MQ=255)
    tGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGcc                                                                          <  2:611310/90‑1 (MQ=255)
        cATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACAt                                                                      <  1:36818/90‑1 (MQ=255)
          tCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCa                                                                    <  1:112339/90‑1 (MQ=255)
          tCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCa                                                                    >  1:329940/1‑90 (MQ=255)
                     ccAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGa                                                         <  1:556507/90‑1 (MQ=255)
                     ccAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGa                                                         <  2:620753/90‑1 (MQ=255)
                       aGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAg                                                       <  1:689264/90‑1 (MQ=255)
                         tCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGc                                                     >  1:491151/1‑90 (MQ=255)
                            tGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCAtt                                                  >  2:492272/1‑90 (MQ=255)
                                  gTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTgc                                                                             >  2:836250/1‑57 (MQ=255)
                                  gTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTgc                                                                             <  1:836250/57‑1 (MQ=255)
                                               tCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTg                               >  2:797052/1‑90 (MQ=255)
                                               tCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTg                               >  2:122911/1‑90 (MQ=255)
                                                 cACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGcc                             <  1:162392/90‑1 (MQ=255)
                                                  aCTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCg                            >  2:239903/1‑90 (MQ=255)
                                                  aCTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCg                            >  2:718192/1‑90 (MQ=255)
                                                  aCTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCg                            >  1:249503/1‑90 (MQ=255)
                                                       aCGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGtttt                       >  2:443270/1‑90 (MQ=255)
                                                            ccGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCt                  >  1:644611/1‑90 (MQ=255)
                                                               aCGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCg               >  2:197332/1‑90 (MQ=255)
                                                               aCGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCg               >  2:132191/1‑90 (MQ=255)
                                                                          gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTg    >  2:175710/1‑90 (MQ=255)
                                                                          gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTg    >  2:65662/1‑90 (MQ=255)
                                                                          gTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTg    >  1:35104/1‑90 (MQ=255)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGc                                                     >  2:595408/1‑39 (MQ=38)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGc                                                     <  1:595408/39‑1 (MQ=255)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGcc      >  2:559010/1‑86 (MQ=255)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGcc      <  1:559010/86‑1 (MQ=255)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCt  >  2:55770/1‑90 (MQ=255)
                                                                            gCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCt  >  2:562131/1‑90 (MQ=255)
                                                                                     |                                                                                
GCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCT  >  NZ_CP009273/2425693‑2425858

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATTTTCCAGCGTGCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCTTTCTCG  >  NZ_CP009273/2425681‑2425864
                                                                                                 |                                                                                      
gacggcatacgagatcctctctggtctcgtgggctcggagatgtgtataagagacagGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCC                                                                                      <  SRR3722079.474535/43‑1 (MQ=60)
ATTTTCCAGCGTGCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCC                                                                                      <  SRR3722079.764080/100‑1 (MQ=60)
            GCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCA                                                                          >  SRR3722079.332575/1‑100 (MQ=60)
             tcggagatgtgtataagagacagGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAG                                                                         <  SRR3722079.845200/77‑1 (MQ=60)
                    CATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGT                                                                  <  SRR3722079.37086/100‑1 (MQ=60)
                      TCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAG                                                                <  SRR3722079.113253/100‑1 (MQ=60)
                           AAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGC                                                           >  SRR3722079.495480/1‑100 (MQ=60)
                                 CCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCC                                                     <  SRR3722079.561883/100‑1 (MQ=60)
                                   AGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGG                                                   <  SRR3722079.696511/100‑1 (MQ=60)
                                     ctctctggtctcgtgggctcggagatgtgtataagagacaGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCT                                                 <  SRR3722079.601373/60‑1 (MQ=60)
                                                    GAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCG                                  >  SRR3722079.251592/1‑100 (MQ=60)
                                                             CACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGC                         <  SRR3722079.163817/100‑1 (MQ=60)
                                                              ACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCT                        >  SRR3722079.651273/1‑100 (MQ=60)
                                                                            CGTTGTAACAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTG          >  SRR3722079.35361/1‑100 (MQ=60)
                                                                                    CAGTGCACCTTTTCCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCTTTCTCG  <  SRR3722079.564423/100‑1 (MQ=60)
                                                                                                 |                                                                                      
ATTTTCCAGCGTGCCGTGAGCATCGCTAAACGCCCAGTCATGATCGGTGACGGAATAGTTCCACTCAACGCCCCGACGTTGTAACAGTGCACCTTTTTCCTGCGCCACATCAGCAATGGTAGAAGGCATTTCCGGCTCACCGACAATTGCCGGTTTTTCGCTGCGGAAGATGCCTGCTTTCTCG  >  NZ_CP009273/2425681‑2425864

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