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 F1 I1 R1 66 46.7 5280542 98.1% 5180211 47.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_017635 2,753,670 A→G W20R (TGG→CGG)  ypfG ← DUF1176 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176352,753,6700AG100.0% 114.3 / NA 30W20R (TGG→CGG) ypfGDUF1176 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/19);  total (11/19)

GTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAA  >  NC_017635/2753627‑2753714
                                           |                                            
gTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCgcac                                           >  1:707577/1‑47 (MQ=255)
 tGACCTGCCAGTCGGAAAACGCCCGTTGCGATGGTGCCGCCCGCAcc                                          >  1:1969069/1‑47 (MQ=38)
    ccTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAg                                      <  2:2491505/48‑1 (MQ=255)
    ccTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAg                                      <  1:2226892/48‑1 (MQ=255)
     cTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGa                                     >  2:1664297/1‑48 (MQ=255)
        ccAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGAc                                    <  1:1415129/46‑1 (MQ=255)
          aGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCg                                <  2:15728/48‑1 (MQ=255)
             cGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGc                              >  2:1232883/1‑47 (MQ=255)
              ggAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGca                             <  1:1916832/47‑1 (MQ=255)
               gAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGcaac                           <  1:17877/48‑1 (MQ=255)
               gAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGcaac                           <  1:27966/48‑1 (MQ=255)
                 aaaCGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGcaaca                          >  2:2269249/1‑47 (MQ=255)
                  aaCGCCCGTTGCGCTGGCGCCGCCCGCACCAAAGACGTCGGcaaca                          >  2:2624784/1‑46 (MQ=39)
                   aCGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAg                       <  1:402658/48‑1 (MQ=255)
                       ccGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGa                    >  2:1951618/1‑47 (MQ=255)
                        cGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGaaa                  <  1:740546/48‑1 (MQ=255)
                          ttGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGaaaaa                <  1:2590622/48‑1 (MQ=255)
                           tGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGaaaaa                >  2:92621/1‑47 (MQ=255)
                           tGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGaaaaa                >  2:216359/1‑47 (MQ=255)
                            gcgcTGGTGCCGCCCGCTCCAAAGACGTCGGCAACAAAGCGaaaaaa               <  1:1991767/47‑1 (MQ=255)
                             cgcTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGaaaaaa               >  1:1579432/1‑46 (MQ=255)
                                   tGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAg       >  2:2438194/1‑48 (MQ=255)
                                   tGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAg       <  1:1345633/48‑1 (MQ=255)
                                     ccgcccgcACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAgg      <  1:2275928/47‑1 (MQ=255)
                                     ccgcccgcACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAgg      <  2:464892/47‑1 (MQ=255)
                                      cgccCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAGAGAAGGaaa   <  1:1402725/48‑1 (MQ=255)
                                       gccCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGaaa   <  1:2529689/48‑1 (MQ=255)
                                       gccCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGaaa   <  2:1397557/48‑1 (MQ=255)
                                       gccCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGaaa   <  1:157500/48‑1 (MQ=255)
                                       gccCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGaaa   <  1:1567168/48‑1 (MQ=255)
                                       gccCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGaaa   <  2:369332/48‑1 (MQ=255)
                                       gccCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGaa    >  1:2252120/1‑47 (MQ=255)
                                        cccGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGaaaa  <  2:2292857/48‑1 (MQ=255)
                                           |                                            
GTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAA  >  NC_017635/2753627‑2753714

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGGTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGAT  >  NC_017635/2753625‑2753720
                                             |                                                  
AGGTTCCCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCA                                               <  M01420:84:000000000‑ABVEH:1:1110:26680:21776/51‑1 (MQ=60)
AGGTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACC                                                >  M01420:84:000000000‑ABVEH:1:1107:24087:4983/1‑50 (MQ=60)
 GGTGACCTGCCAGTCGGAAAACGCCCGTTGCGATGGTGCCGCCCGCACCA                                               >  M01420:84:000000000‑ABVEH:1:2103:27161:21723/1‑50 (MQ=60)
     CCCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGAC                                          <  M01420:84:000000000‑ABVEH:1:2106:9613:3244/51‑1 (MQ=60)
         GCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGT                                        <  M01420:84:000000000‑ABVEH:1:1113:10944:4267/49‑1 (MQ=60)
               CGGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAAC                                 <  M01420:84:000000000‑ABVEH:1:2103:23283:11549/50‑1 (MQ=60)
                GGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAA                               <  M01420:84:000000000‑ABVEH:1:1101:8834:7240/51‑1 (MQ=60)
                GGAAAACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAA                               <  M01420:84:000000000‑ABVEH:1:1101:22743:9447/51‑1 (MQ=60)
                    AACGCCCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCG                           <  M01420:84:000000000‑ABVEH:1:1104:25659:14562/51‑1 (MQ=60)
                         CCGTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAA                      <  M01420:84:000000000‑ABVEH:1:1107:27012:12150/51‑1 (MQ=60)
                           GTTGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAA                    <  M01420:84:000000000‑ABVEH:1:2109:16126:6106/51‑1 (MQ=60)
                             TGCGCTGGTGCCGCCCGCTCCAAAGACGTCGGCAACAAAGCGAAAAAAAA                   <  M01420:84:000000000‑ABVEH:1:2103:21458:27864/50‑1 (MQ=60)
                             TGCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAA                    >  M01420:84:000000000‑ABVEH:1:1114:23410:13672/1‑49 (MQ=60)
                              GCGCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAG                  <  M01420:84:000000000‑ABVEH:1:2112:17906:24196/50‑1 (MQ=60)
                                GCTGGTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAG                  >  M01420:84:000000000‑ABVEH:1:2112:22773:22117/1‑48 (MQ=60)
                                    GTGCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGA           <  M01420:84:000000000‑ABVEH:1:1112:10649:14731/51‑1 (MQ=60)
                                      GCCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAA          <  M01420:84:000000000‑ABVEH:1:2106:20964:14182/50‑1 (MQ=60)
                                       CCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAA         >  M01420:84:000000000‑ABVEH:1:2106:13070:9555/1‑50 (MQ=60)
                                       CCGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAGAGAAGGAAAAT       <  M01420:84:000000000‑ABVEH:1:1112:19907:26697/51‑1 (MQ=57)
                                        CGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAAT       <  M01420:84:000000000‑ABVEH:1:1114:10491:11351/51‑1 (MQ=60)
                                        CGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAAT       <  M01420:84:000000000‑ABVEH:1:2108:24228:17519/51‑1 (MQ=60)
                                        CGCCCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAAT       <  M01420:84:000000000‑ABVEH:1:1102:20158:12154/51‑1 (MQ=60)
                                           CCGCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCG    <  M01420:84:000000000‑ABVEH:1:2103:27395:8273/51‑1 (MQ=60)
                                             GCACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGAT  <  M01420:84:000000000‑ABVEH:1:2113:25502:23933/51‑1 (MQ=60)
                                             |                                                  
AGGTGACCTGCCAGTCGGAAAACGCCCGTTGCGCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGAT  >  NC_017635/2753625‑2753720

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