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
A3 F53 I1 R1 62 81.2 3465250 98.6% 3416736 145.3

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% 79.5 / NA 23intergenic (+90/+20)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/4);  total (19/4)

CCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCAT  >  NC_002947_CJ‑RC/6016636‑6016841
                                                                |                                                                                                                                             
cccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCACCCCTTTTATATGGTTGCCGCCAGCCCCGCCACCAACAAGCGAGTGGAAAGTTGTTCTTTAAGACCCTGc                                                            >  5:348507/1‑148 (MQ=37)
 ccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                                           >  1:102483/1‑148 (MQ=255)
        cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                                    >  1:318709/1‑148 (MQ=255)
        cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                                    >  7:236698/1‑148 (MQ=255)
         acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                                    >  3:10731/1‑147 (MQ=255)
                cagcaCAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCa                                            >  8:347481/1‑148 (MQ=255)
                  gcaCAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATc                                          >  2:171731/1‑148 (MQ=255)
                        gTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGccc                                     >  2:202700/1‑147 (MQ=255)
                          ctcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAg                                  >  3:376040/1‑148 (MQ=255)
                            cGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttct                                                                           <  8:207731/105‑1 (MQ=255)
                            cGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttct                                                                           >  7:207731/1‑105 (MQ=255)
                              aaCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTgg                              >  1:297623/1‑148 (MQ=255)
                                     tgCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCa                        <  2:297623/147‑1 (MQ=255)
                                     tgCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCa                        <  6:403756/147‑1 (MQ=255)
                                        cGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCACAGGTGGCTTTCATAcc                    >  1:105804/1‑148 (MQ=255)
                                           aCCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCt                 >  2:204712/1‑148 (MQ=255)
                                              tgtgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCg              >  2:237700/1‑148 (MQ=255)
                                              tgtgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTATAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACACCTCCg              >  3:37301/1‑148 (MQ=255)
                                                tgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCACCAACAACAAGCGCTGGTAAGGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGc            >  3:393431/1‑148 (MQ=255)
                                                    aGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCg        >  2:33999/1‑148 (MQ=255)
                                                        ggCGTGCCCGCGAATGGCCGCACCCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCCACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGc    >  5:360443/1‑148 (MQ=255)
                                                         gCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCa   <  7:339167/148‑1 (MQ=255)
                                                          cGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCAt  >  6:216157/1‑148 (MQ=255)
                                                                |                                                                                                                                             
CCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCAT  >  NC_002947_CJ‑RC/6016636‑6016841

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCT  >  NC_002947_CJ‑RC/6016597‑6016849
                                                                                                       |                                                                                                                                                     
caaaccacaaaaccccccccccccgcctccggggcccccccccccccaccgggcccccccccgcccccccccgtgggcccgcccggggggggggGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAAC                                                                                                        <  NB501016:282:H2V7GAFX2:1:11210:16881:17610/57‑1 (MQ=60)
            aacggcaaccgggcggggggggacggcccgcccccgctgggtgcccgcgtggtgtcgcgcttggccctttgcttggggtgggggggGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGG                                                                                            <  NB501016:282:H2V7GAFX2:4:21612:14852:13281/65‑1 (MQ=60)
                             ccgggccccccccgccccaccgggcgcccaccggccccccccggggggcccccccgtggggcgggGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAG                                                                           <  NB501016:282:H2V7GAFX2:3:11509:25851:4395/86‑1 (MQ=60)
                                     CGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCACCCCTTTTATATGGTTGCCGCCAGCCCCGCCACCAACAAGCGAGTGGAAAGTTGTTCTTTAAGACCCTGCG                                                                   >  NB501016:282:H2V7GAFX2:3:21510:3594:18458/1‑151 (MQ=60)
                                     CGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGGTGGTAAAGTTCTTCTTTAATCCCCTGCC                                                                   >  NB501016:282:H2V7GAFX2:3:21601:19738:5299/1‑151 (MQ=60)
                                      GCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGG                                                                  >  NB501016:282:H2V7GAFX2:1:11204:13157:14782/1‑151 (MQ=60)
                                             CCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCC                                                           >  NB501016:282:H2V7GAFX2:1:21204:4829:10450/1‑151 (MQ=60)
                                             CCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCC                                                           >  NB501016:282:H2V7GAFX2:4:21403:8046:4866/1‑151 (MQ=60)
                                              CCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCC                                                           >  NB501016:282:H2V7GAFX2:2:11102:8247:17426/1‑150 (MQ=60)
                                                   ggtacggcacagcgttcgagacctgtgctgcACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGC                                                     >  NB501016:282:H2V7GAFX2:4:11402:12796:6592/32‑151 (MQ=60)
                                                               GTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGG                                         >  NB501016:282:H2V7GAFX2:2:21306:14884:7168/1‑151 (MQ=60)
                                                                 CTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTC                                                                                 >  NB501016:282:H2V7GAFX2:4:11610:21282:8748/1‑109 (MQ=60)
                                                                   cccccgggggccgtcagggggggggggggggGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGC                                     <  NB501016:282:H2V7GAFX2:1:11212:4988:9921/120‑1 (MQ=60)
                                                                   CGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGC                                     >  NB501016:282:H2V7GAFX2:1:21201:17674:1209/1‑151 (MQ=60)
                                                                             GCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCACAGGTGGCTTTCATACCC                           >  NB501016:282:H2V7GAFX2:1:11205:22033:6320/1‑151 (MQ=60)
                                                                               CGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTaagccccatcgggttcccccaaaacgaagaaccaacctgataacaccattaaa                          >  NB501016:282:H2V7GAFX2:4:11602:2458:4972/1‑97 (MQ=60)
                                                                                   CCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTATAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACACCTCCGG                     >  NB501016:282:H2V7GAFX2:2:11107:5082:7921/1‑151 (MQ=60)
                                                                                     TGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCACCAACAACAAGCGCTGGTAAGGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCC                   >  NB501016:282:H2V7GAFX2:2:21309:18648:2437/1‑151 (MQ=60)
                                                                                           gggGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTG             <  NB501016:282:H2V7GAFX2:3:11603:5251:14833/148‑1 (MQ=60)
                                                                                             CGGGCGTGCCCGCGAATGGCCGCACCCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCCACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCC           >  NB501016:282:H2V7GAFX2:3:21512:1863:18728/1‑151 (MQ=60)
                                                                                               GGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATC         <  NB501016:282:H2V7GAFX2:4:21507:5967:9862/151‑1 (MQ=60)
                                                                                                      CCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCT  >  NB501016:282:H2V7GAFX2:4:21508:12033:5090/1‑151 (MQ=60)
                                                                                                       |                                                                                                                                                     
GCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCT  >  NC_002947_CJ‑RC/6016597‑6016849

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: