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 F11 I1 R1 56 91.5 3907174 98.4% 3844659 146.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 1,126,647 +C intergenic (‑101/+48) tdcG‑II ← / ← gcvP‑I L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC1,126,6451.C94.7% 60.9 / ‑3.3 19intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (13/5);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCG  >  NC_002947_CJ‑RC/1126503‑1126780
                                                                                                                                                |                                                                                                                                       
gggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCTGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgc                                                                                                                                    <  8:389548/148‑4 (MQ=255)
                            acaCGCTAAGTTACTTGGGGACCACCTGCCTGGTTGTTGGGTGCAAATTGAGCCCTTTCGGGGGGGTTTAACCCGCTAATGGGTTGTTATACCCACAATCGAATTGGGGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCa                                                                                                         <  6:414223/146‑1 (MQ=255)
                                        aCATGTGCACCTCCTGCCAGGTTGTTGGGTGCACTTTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                             <  2:86785/148‑1 (MQ=38)
                                                                                          gcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg                                          >  1:398197/1‑149 (MQ=255)
                                                                                          gcgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAgg                                          >  7:278023/1‑149 (MQ=255)
                                                                                              gTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCCCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGCCGGGCAGGCGCACACCAGGTTg                                       >  5:382241/1‑148 (MQ=255)
                                                                                                 tACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACCCCAAGTTGCGGt                                   >  1:76011/1‑149 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                >  4:299293/1‑149 (MQ=255)
                                                                                                    ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                >  7:382652/1‑149 (MQ=255)
                                                                                                      gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                               >  8:404012/1‑148 (MQ=255)
                                                                                                              ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                      >  1:126012/1‑149 (MQ=255)
                                                                                                               tttgTATACCCACCATCCCTTTCGCGGTTAACCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                      <  6:449884/145‑1 (MQ=255)
                                                                                                                     taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc               >  1:286525/1‑149 (MQ=255)
                                                                                                                     taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc                >  5:266031/1‑148 (MQ=255)
                                                                                                                       gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGt                                         >  6:133629/1‑121 (MQ=255)
                                                                                                                        ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                    <  7:191699/77‑1 (MQ=255)
                                                                                                                        ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                    >  8:191699/1‑77 (MQ=255)
                                                                                                                        ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGtt                                        <  5:133629/121‑1 (MQ=255)
                                                                                                                                gCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGctcaggtg    >  8:242973/1‑142 (MQ=255)
                                                                                                                                   ttCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg  <  6:266031/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                       
GGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCG  >  NC_002947_CJ‑RC/1126503‑1126780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCC  >  NC_002947_CJ‑RC/1126529‑1126793
                                                                                                                       |                                                                                                                                                  
cacccccccctttcacttggccccccccgcccggttttttggttgcaattcaccccctgtggggggggttaaccccgcaatggcgtcgtaaacccccccccccaatgcgggttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAG                                                                                                                     <  NB501016:282:H2V7GAFX2:4:21509:1360:7785/33‑1 (MQ=60)
           gcccctggggccccccccccgtggtttggggggtcacttcgcccggtgggggggggggttcccccgctatggggtcgctacccccaccacccctggggggtaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGG                                                                                                            <  NB501016:282:H2V7GAFX2:1:11211:9490:7097/42‑1 (MQ=60)
              ccctgggccccccccgccgtgggtggggtgcgcactgcccgccggctcgggggggtcccccccgtactggcggcgaaaaccccccaccgacggggggtaaacccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGC                                                                                                       <  NB501016:282:H2V7GAFX2:2:11307:17160:12185/47‑1 (MQ=60)
                    gcccggcggggttgtggggggggccccatggccccggtcggggggggtcccccccgggtcgggcgtttaaacccccccaccgccggggggtttaacccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAG                                                                                                 <  NB501016:282:H2V7GAFX2:3:21412:2371:12527/53‑1 (MQ=60)
                                     ggggggggtccattgccccgtgtggggggggttaccccggcttggcggtttctaaccccccctccaattgggggtaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCT                                                                                 <  NB501016:282:H2V7GAFX2:1:21209:20047:13840/69‑1 (MQ=60)
                                                               ggggggttacccggttatggggttttaagcccccctagcccttcggggttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGT                                                      <  NB501016:282:H2V7GAFX2:3:11502:11134:15742/96‑1 (MQ=60)
                                                               CGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGT                                                      >  NB501016:282:H2V7GAFX2:1:21210:21922:18643/1‑151 (MQ=60)
                                                               CGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGT                                                      >  NB501016:282:H2V7GAFX2:4:21404:9727:12352/1‑151 (MQ=60)
                                                                 gggtttcccccgctattgggttggtaacgcacccaccgatttcgcgttaaaacCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTG                                                    <  NB501016:282:H2V7GAFX2:3:11506:21198:5456/98‑1 (MQ=60)
                                                                   GGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCCCCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGCCGGGCAGGCGCACACCAGGTTGCG                                                  >  NB501016:282:H2V7GAFX2:3:21509:7901:7439/1‑151 (MQ=60)
                                                                      TTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACCCCAAGTTGCGGTC                                               >  NB501016:282:H2V7GAFX2:1:11111:23686:9574/1‑151 (MQ=60)
                                                                         CCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCC                                            >  NB501016:282:H2V7GAFX2:4:21507:2833:6045/1‑151 (MQ=60)
                                                                          ggtgttgtggggttttaccccccctatggcttgggggtttaaccCCGCTCCTACACTGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCG                                           <  NB501016:282:H2V7GAFX2:3:11507:3494:12574/107‑1 (MQ=60)
                                                                          ccgctatggggtttttaaccccacccgcccttgcgggttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCG                                           <  NB501016:282:H2V7GAFX2:4:21609:17495:15747/107‑1 (MQ=60)
                                                                           CGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACC                                                                                                                                                    >  NB501016:282:H2V7GAFX2:4:11501:17522:16663/1‑44 (MQ=60)
                                                                                   CGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTG                                  >  NB501016:282:H2V7GAFX2:1:11206:6812:12990/1‑151 (MQ=60)
                                                                                    gatggtggctatacctcCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGT                                 <  NB501016:282:H2V7GAFX2:1:21306:15468:8251/134‑1 (MQ=60)
                                                                                          ATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCC                           >  NB501016:282:H2V7GAFX2:3:21403:5356:11062/1‑151 (MQ=60)
                                                                                          ATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCC                           >  NB501016:282:H2V7GAFX2:1:21106:21894:5934/1‑151 (MQ=60)
                                                                                            AGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTG                                                    <  NB501016:282:H2V7GAFX2:3:11508:18303:12479/124‑1 (MQ=60)
                                                                                             GCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGG                                                                                                <  NB501016:282:H2V7GAFX2:4:11603:25163:20227/79‑1 (MQ=60)
                                                                                                          gggcggtaaaacCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCC            <  NB501016:282:H2V7GAFX2:4:11506:25725:13792/138‑1 (MQ=60)
                                                                                                             ggtttaaacCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTA        <  NB501016:282:H2V7GAFX2:4:21511:11942:20080/142‑1 (MQ=60)
                                                                                                                   cccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCC  <  NB501016:282:H2V7GAFX2:1:21212:5510:2972/148‑1 (MQ=60)
                                                                                                                       |                                                                                                                                                  
GAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCC  >  NC_002947_CJ‑RC/1126529‑1126793

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