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
A4 F13 I1 R1 59 72.5 3655346 98.4% 3596860 137.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 6,016,700 T→C intergenic (+90/+20) xpt → / ← PP_5266 xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC6,016,7000TC100.0% 93.1 / NA 27intergenic (+90/+20)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/12);  total (15/12)

ACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGT  >  NC_002947_CJ‑RC/6016633‑6016837
                                                                   |                                                                                                                                         
aCGCCCGCTCCCACAGGTGCAGCACACGTCTCGAACGCAGTGCCCTACCTGTGGGAGCGGGGGGGCCCCCCAATGGCCGCAAACCTTTTATATGGTTTGCGCCAGCCCCGCCAAAAACAAACCCTGGTAAAGTTATTATTTAAGccc                                                            >  4:53029/1‑147 (MQ=255)
  gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGCATGGCCGCAACCCGTTTTTATGGTTGCCGCCAGCCCCGCCAACAACAACCGCGGGTAAAGTTCTTCTTTAAGCCCCTGc                                                        >  7:417315/1‑149 (MQ=38)
  gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                                        >  6:392580/1‑149 (MQ=255)
      gCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                                       <  1:259609/146‑1 (MQ=255)
      gCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                                       <  1:266237/146‑1 (MQ=255)
      gCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                                       >  2:259609/1‑146 (MQ=255)
      gCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                                       >  2:266237/1‑146 (MQ=255)
           caccGGTGCAGCACCGGTCTCTAAGGGTGTGCGGTACGTGGGGTCGGGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                                               <  8:417315/149‑1 (MQ=255)
              aGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAAc                                            >  8:379860/1‑149 (MQ=255)
                           gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                                        >  4:152999/1‑124 (MQ=255)
                           gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                                        <  3:152999/124‑1 (MQ=255)
                              tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGt                            >  3:196026/1‑149 (MQ=255)
                                gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCa                               >  3:327989/1‑144 (MQ=255)
                                gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCa                               <  4:327989/144‑1 (MQ=255)
                                gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTgg                          >  5:168868/1‑149 (MQ=255)
                                   cGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCtt                       <  7:379860/149‑1 (MQ=255)
                                      tgtgCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGTCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCa                    >  7:61868/1‑149 (MQ=255)
                                           cGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAccc               >  8:320781/1‑149 (MQ=255)
                                                cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCg          <  4:196026/149‑1 (MQ=255)
                                                 tgtgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGAtcttct                                                                    <  6:128318/90‑1 (MQ=255)
                                                 tgtgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGAtcttct                                                                    >  5:128318/1‑90 (MQ=255)
                                                   tgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGc        <  6:286121/148‑1 (MQ=255)
                                                        gggggcgggcCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGt  <  7:258596/147‑1 (MQ=255)
                                                          gggCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCAt                                       <  2:214736/110‑1 (MQ=255)
                                                          gggCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCAt                                       >  1:214736/1‑110 (MQ=255)
                                                          gggCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCgg   <  3:262684/146‑1 (MQ=255)
                                                          gggCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCgg   >  4:262684/1‑146 (MQ=255)
                                                                   |                                                                                                                                         
ACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGT  >  NC_002947_CJ‑RC/6016633‑6016837

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGT  >  NC_002947_CJ‑RC/6016622‑6016837
                                                                              |                                                                                                                                         
gtagttacaaaaggcctctcctcttgtgtcacccattgtgtttaagggtgttcggtatgggggtgcgGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCTCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCT                                                                    <  NB501016:282:H2V7GAFX2:1:11307:2392:7519/83‑1 (MQ=60)
            CGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGCATGGCCGCAACCCGTTTTTATGGTTGCCGCCAGCCCCGCCAACAACAACCGCGGGTAAAGTTCTTCTTTAAGCCCCTGCC                                                       >  NB501016:282:H2V7GAFX2:4:21603:8967:19709/1‑151 (MQ=60)
               CCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGC                                                     <  NB501016:282:H2V7GAFX2:1:21105:1280:12814/150‑1 (MQ=60)
               CCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGC                                                     <  NB501016:282:H2V7GAFX2:1:21106:26693:12812/150‑1 (MQ=60)
                        ggcccccccccgtcgctccccgggggggcccgggcgtggggcgggGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACT                                           <  NB501016:282:H2V7GAFX2:2:11109:11550:12606/106‑1 (MQ=60)
                          ggtacggcacagcgttcgagacctgtgctgcACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGC                                         >  NB501016:282:H2V7GAFX2:3:21409:5116:9943/32‑151 (MQ=60)
                            TGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA                                                                                                                                      <  NB501016:282:H2V7GAFX2:4:21409:4867:15148/56‑1 (MQ=60)
                                     GGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCG                                                       <  NB501016:282:H2V7GAFX2:2:11301:9463:9885/126‑1 (MQ=60)
                                        CTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTG                           >  NB501016:282:H2V7GAFX2:2:11308:7610:3667/1‑151 (MQ=60)
                                          CGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGC                         >  NB501016:282:H2V7GAFX2:3:11603:5838:5659/1‑151 (MQ=60)
                                          CGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAG                              >  NB501016:282:H2V7GAFX2:2:21207:6798:9614/1‑146 (MQ=60)
                                          CGAACGCTGTGCCGTACCTGTGGGAGCGGCCGAGCCCGCGAATGGCCCCAACCCATAAATCTCCTTGCCGCCAGCCCCCCCAACAACAAGCGCAGGAAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTACATCCGCCCCAGGacgc                         >  NB501016:282:H2V7GAFX2:1:21111:1375:4407/1‑147 (MQ=60)
                                             ACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTT                      <  NB501016:282:H2V7GAFX2:4:21510:3711:5875/151‑1 (MQ=60)
                                               ggggcccctcccggggggggggGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCA                    <  NB501016:282:H2V7GAFX2:1:11112:18784:12970/129‑1 (MQ=60)
                                                CTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGTCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCAT                   >  NB501016:282:H2V7GAFX2:4:11410:13955:4770/1‑151 (MQ=60)
                                                           CTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGATCTTCTT                                                                   >  NB501016:282:H2V7GAFX2:3:11508:22780:16211/1‑92 (MQ=60)
                                                                   GCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCC                                     >  NB501016:282:H2V7GAFX2:1:11310:5603:16649/1‑114 (MQ=60)
                                                                    CGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGT  <  NB501016:282:H2V7GAFX2:2:21107:5171:1152/148‑1 (MQ=60)
                                                                              |                                                                                                                                         
GCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGT  >  NC_002947_CJ‑RC/6016622‑6016837

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