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 F64 I1 R1 60 70.4 2989830 98.5% 2944982 146.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 4,586,164 +C intergenic (+140/+75) PP_4061 → / ← PP_4063 hypothetical protein/long‑chain fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC4,586,1631.C95.8% 79.6 / ‑2.9 24intergenic (+139/+76)PP_4061/PP_4063hypothetical protein/long‑chain fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (9/14);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGC  >  NC_002947_CJ‑RC/4586018‑4586306
                                                                                                                                                    |                                                                                                                                                 
ggCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGAcc                                                                                                                                                 <  4:159102/149‑2 (MQ=255)
            gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCtt                                                                                                                                      <  7:151296/148‑1 (MQ=39)
            gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCtt                                                                                                                                      <  6:87904/148‑1 (MQ=39)
            gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCtt                                                                                                                                      <  4:98813/148‑1 (MQ=39)
              ccACGGCACCGGCGTGAACAGCTTATCTGCGCTGATCCATTGAAAGTTTGAGGATGACCGGGTGCTTTGCAAGTTGAAGATCTTGCATTGATCGGGCCTGCCTCTTCGCGTGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGa                                                                                                                                    <  5:10460/148‑1 (MQ=255)
                aCGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATg                                                                                                                                 <  3:138156/149‑1 (MQ=39)
                               aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTg                                                                                                                  <  2:299606/148‑1 (MQ=39)
                                 aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                               <  1:169729/149‑1 (MQ=39)
                                 aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa                                                                                                               <  5:93977/149‑1 (MQ=39)
                                      aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGa                                                                                                            <  2:333196/147‑1 (MQ=39)
                                           cgcgGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGtt                                                                                                     <  1:122808/149‑1 (MQ=38)
                                           cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGc                                                                                                          <  7:380244/144‑1 (MQ=38)
                                           cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCggg                                                                                                       <  3:30001/147‑1 (MQ=38)
                                           cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg                                                                                                        <  7:276777/146‑1 (MQ=38)
                                                ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc                                                                                                  <  4:163924/147‑1 (MQ=38)
                                                                                                                           gggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCCGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACCCTGATCTCCCGCATg                     >  6:161386/2‑149 (MQ=38)
                                                                                                                            ggTAAACCCGCTCCTACCAGGGGCAGCGCCGCATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGc                    >  5:298231/1‑149 (MQ=18)
                                                                                                                            ggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCGGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGc                    >  3:238626/1‑149 (MQ=255)
                                                                                                                               tAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg                   >  5:34619/1‑148 (MQ=255)
                                                                                                                                aaaCCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg                  >  3:32329/1‑148 (MQ=255)
                                                                                                                                   cccGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt              >  7:147571/1‑149 (MQ=255)
                                                                                                                                   cccGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACt              >  3:199982/1‑149 (MQ=255)
                                                                                                                                          cTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAc       >  6:54044/1‑149 (MQ=255)
                                                                                                                                          cTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAc       >  4:35561/1‑149 (MQ=255)
                                                                                                                                               aGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGc  >  4:124519/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGC  >  NC_002947_CJ‑RC/4586018‑4586306

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
AGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTT  >  NC_002947_CJ‑RC/4586029‑4586295
                                                                                                                                          |                                                                                                                                   
AGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTG                                                                                                                         <  NB501016:282:H2V7GAFX2:4:11605:21098:5378/150‑1 (MQ=60)
  ACCACGGCACCGGCGTGAACAGCTTATCTGCGCTGATCCATTGAAAGTTTGAGGATGACCGGGTGCTTTGCAAGTTGAAGATCTTGCATTGATCGGGCCTGCCTCTTCGCGTGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAA                                                                                                                       <  NB501016:282:H2V7GAFX2:3:11403:25329:1361/150‑1 (MQ=60)
    CACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGC                                                                                                                    <  NB501016:282:H2V7GAFX2:2:11303:13010:12625/151‑1 (MQ=60)
                     CAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAG                                                                                                   <  NB501016:282:H2V7GAFX2:1:11308:25376:10014/151‑1 (MQ=60)
                     CAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAG                                                                                                   <  NB501016:282:H2V7GAFX2:3:11507:16744:5172/151‑1 (MQ=60)
                               GCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGT                                                                                           <  NB501016:282:H2V7GAFX2:2:11106:25812:18205/149‑1 (MQ=60)
                               GCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGG                                                                                            <  NB501016:282:H2V7GAFX2:4:21504:4805:14137/148‑1 (MQ=60)
                               GCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCG                                                                                              <  NB501016:282:H2V7GAFX2:4:21612:10586:6072/146‑1 (MQ=60)
                               GCGCGGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTT                                                                                         <  NB501016:282:H2V7GAFX2:1:11211:13162:6393/151‑1 (MQ=60)
                                                                                                               GGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCGGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCG        >  NB501016:282:H2V7GAFX2:2:21111:14664:2324/1‑151 (MQ=60)
                                                                                                                 GGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCC                                                                                     <  NB501016:282:H2V7GAFX2:3:21510:13057:10567/73‑1 (MQ=60)
                                                                                                                   GTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGG       >  NB501016:282:H2V7GAFX2:3:11407:12962:16341/1‑150 (MQ=60)
                                                                                                                    TAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGA      >  NB501016:282:H2V7GAFX2:2:11107:23947:7679/1‑150 (MQ=60)
                                                                                                                       ACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTT  >  NB501016:282:H2V7GAFX2:2:21103:23075:3276/1‑151 (MQ=60)
                                                                                                                       ACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTT  >  NB501016:282:H2V7GAFX2:4:11604:16274:10475/1‑151 (MQ=60)
                                                                                                                                          |                                                                                                                                   
AGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTT  >  NC_002947_CJ‑RC/4586029‑4586295

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