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 1,499,508 +C intergenic (+85/‑79) trpS → / → zapE tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC1,499,5061.C94.4% 53.2 / ‑2.8 18intergenic (+83/‑81)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (11/6);  total (12/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

CGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCG  >  NC_002947_CJ‑RC/1499360‑1499649
                                                                                                                                                     |                                                                                                                                               
cGCTGACGCCCGTGATGCGGCCATTCAGAGCCTGAAGGAAGGCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgc                                                                                                                                                <  4:181102/149‑1 (MQ=25)
 gcggAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgcc                                                                                                                                               <  5:281157/146‑1 (MQ=37)
                                                                                      gCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACt                                                           >  7:272131/1‑149 (MQ=25)
                                                                                           cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                       >  5:376560/1‑148 (MQ=37)
                                                                                            ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                       >  8:442517/1‑147 (MQ=37)
                                                                                            ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                       >  7:356960/1‑147 (MQ=37)
                                                                                            ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATAAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                       >  5:382995/1‑147 (MQ=25)
                                                                                                   aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATAAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt                                              <  6:382995/149‑1 (MQ=37)
                                                                                                    cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc                                               >  1:1129/1‑147 (MQ=37)
                                                                                                          ccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                                       <  7:442517/149‑1 (MQ=37)
                                                                                                           cgcgGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                        >  7:352543/1‑147 (MQ=37)
                                                                                                              ggTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                                >  2:85359/1‑120 (MQ=25)
                                                                                                              ggTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                                <  1:85359/120‑1 (MQ=25)
                                                                                                                     ggcaccGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                                       >  7:275586/6‑138 (MQ=38)
                                                                                                                     ggcaccGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                                       <  8:275586/133‑1 (MQ=38)
                                                                                                                         ccGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTc                         <  2:1129/147‑1 (MQ=39)
                                                                                                                                     gtgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCa            >  7:382525/8‑149 (MQ=255)
                                                                                                                                             gggCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                                                      >  4:313804/1‑99 (MQ=255)
                                                                                                                                             gggCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                                                      <  3:313804/99‑1 (MQ=255)
                                                                                                                                               gCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACgcggcg  >  7:355747/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCG  >  NC_002947_CJ‑RC/1499360‑1499649

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGC  >  NC_002947_CJ‑RC/1499393‑1499650
                                                                                                                      |                                                                                                                                              
gggggggtttcccgtgggggcctcctcctccctcccccccggggggggggggggggcccccccaaccccggggcgggggggccgccccgggttggggggggggggggggccccccccCCGCTCCCACAGACTGATTGCCATTAAGCGGGCC                                                                                                                <  NB501016:282:H2V7GAFX2:4:21609:12987:16950/34‑1 (MQ=60)
    ccacctttgcccccccctcccgtctggggccccaattggggggggggggggggcccccaaaaccggggcggggggttcgggccccggggttgggggggttggggggcccccccCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATC                                                                                                            <  NB501016:282:H2V7GAFX2:3:11608:20070:17245/38‑1 (MQ=60)
          cccccccccccccctcccggggcccccactcggggggggggggggcgccccccaaacccggcccggggggccccgccccggggttggggggggggggggccccccccCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACA                                                                                                      <  NB501016:282:H2V7GAFX2:4:21605:19310:16482/44‑1 (MQ=60)
                  ccccccgccccccccccccaccgggggggggggggggccccccccaacccgcccccgggggcccccccccgggtgggggggggggggccccccccccccCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGC                                                                                              <  NB501016:282:H2V7GAFX2:2:21211:5840:2962/52‑1 (MQ=60)
                        tgttcggggccccaactggggggggggggggggccccctaaatcggggggggggtctgccgcccggggtgtgggggtgtgggcgggcccccccCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGT                                                                                          <  NB501016:282:H2V7GAFX2:3:11506:16948:10785/56‑1 (MQ=60)
                                       gggggggggggggggccgccccccacccccggcgggggggccccccccgcgggggggggggggggcgcccccccccccCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTA                                                                         <  NB501016:282:H2V7GAFX2:2:11111:16063:4115/73‑1 (MQ=60)
                                                    GGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACACCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTC                                                           >  NB501016:282:H2V7GAFX2:4:21406:9027:12040/1‑151 (MQ=60)
                                                       cgcccccccgcccgggggggggggttggggggcgggggggggggggggcccccccccaacgcCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCC                                                         <  NB501016:282:H2V7GAFX2:2:21306:11468:4711/89‑1 (MQ=60)
                                                         GCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTT                                                       >  NB501016:282:H2V7GAFX2:3:21511:11509:16774/1‑150 (MQ=60)
                                                          CCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTT                                                       >  NB501016:282:H2V7GAFX2:4:21506:10703:16691/1‑149 (MQ=60)
                                                          CCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATAAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTT                                                       >  NB501016:282:H2V7GAFX2:3:21512:20089:17296/1‑149 (MQ=60)
                                                                 AACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCT                                               >  NB501016:282:H2V7GAFX2:1:11101:13911:5033/1‑150 (MQ=60)
                                                                        GCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCA                                       <  NB501016:282:H2V7GAFX2:4:21607:21443:18418/151‑1 (MQ=60)
                                                                         CCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGA                                        >  NB501016:282:H2V7GAFX2:4:21506:22802:4092/1‑149 (MQ=60)
                                                                           GCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATG                                                               <  NB501016:282:H2V7GAFX2:1:11202:9188:4425/124‑1 (MQ=60)
                                                                                   TGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCA                                       >  NB501016:282:H2V7GAFX2:4:21407:14938:2174/1‑140 (MQ=60)
                                                                                        ggggggggggggggggcgcccccccccccCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTC                          <  NB501016:282:H2V7GAFX2:2:21306:1324:15195/120‑1 (MQ=60)
                                                                                                  CGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAA                              <  NB501016:282:H2V7GAFX2:3:21510:8155:13794/134‑1 (MQ=60)
                                                                                                   ggtgttcgCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCAT            >  NB501016:282:H2V7GAFX2:4:21510:14787:13279/9‑151 (MQ=60)
                                                                                                         gggcggcaccccCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACG        <  NB501016:282:H2V7GAFX2:1:21106:2426:10940/138‑1 (MQ=60)
                                                                                                           CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTA                                                      <  NB501016:282:H2V7GAFX2:2:21205:20844:4648/101‑1 (MQ=60)
                                                                                                             GGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGC  >  NB501016:282:H2V7GAFX2:4:21506:15339:13253/1‑151 (MQ=60)
                                                                                                             GGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACG                                                 >  NB501016:282:H2V7GAFX2:1:21111:22646:13791/1‑104 (MQ=60)
                                                                                                               CACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGC  >  NB501016:282:H2V7GAFX2:2:21203:8550:10540/1‑149 (MQ=60)
                                                                                                                      |                                                                                                                                              
GAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCGC  >  NC_002947_CJ‑RC/1499393‑1499650

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: