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
A0 F0 I1 R1 149 26.7 2337422 87.5% 2045244 105.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
JC AM260479 1,464,441 Δ37 bp coding (125‑161/552 nt) h16_A1354 → conserved hypothetical protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AM260479 = 14644402 (0.080)23 (0.880) 17/208 0.4 95.8% coding (124/552 nt) h16_A1354 conserved hypothetical protein
?AM260479 1464478 = 0 (0.000)coding (162/552 nt) h16_A1354 conserved hypothetical protein

GGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AM260479/1464308‑1464440
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGACTGGCCTGATAATGACGGCTGAGCAGGCATGGCCGCGCCGCCAC  >  AM260479/1464478‑1464610
                                                                                                                                                                                                                                                                          
GGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCG                                                                                                                                     >  2:42572/1‑134
 GAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGC                                                                                                                                    <  2:559925/134‑1
          TGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGC                                                                                                                           >  1:486839/1‑134
              CGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCG                                                                                                                       <  1:362612/134‑1
                    TGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGC                                                                                                                   <  1:509573/132‑1
                                TCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGA                                                                                                     <  2:577823/134‑1
                                        TGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCA                                                                                                       <  2:284647/124‑1
                                        TGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCA                                                                                                       >  1:284647/1‑124
                                                  GTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGC                                                                                   >  1:387850/1‑134
                                                  GTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGAT                                                                                                             >  2:1031781/1‑108
                                                  GTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGAT                                                                                                             <  1:1031781/108‑1
                                                                  ACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGC                                                                    >  2:278627/1‑134
                                                                                  GCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTTAGCCGCCGAATATGTGGC                                                                    >  1:121716/1‑118
                                                                                  GCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTTAGCCGCCGAATATGTGGC                                                                    <  2:121716/118‑1
                                                                                  GCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCC                                                    >  1:1102518/1‑134
                                                                                         GCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGAC                                             <  2:1102518/134‑1
                                                                                                 GCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGAC                                             <  2:894944/126‑1
                                                                                                 GCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGAC                                             >  1:894944/1‑126
                                                                                                   CATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCA                                                                                                       <  2:220191/66‑1
                                                                                                   CATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCA                                                                                                       >  1:677024/1‑66
                                                                                                   CATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCA                                                                                                       <  2:677024/66‑1
                                                                                                   CATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCA                                                                                                       >  1:220191/1‑66
                                                                                                                                    GCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGACTGGCCTGATAATGACGGCTGAGCAGGCATGGCCGCGCCGCCAC  <  1:42572/134‑1
                                                                                                                                                                                                                                                                          
GGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AM260479/1464308‑1464440
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGACTGGCCTGATAATGACGGCTGAGCAGGCATGGCCGCGCCGCCAC  >  AM260479/1464478‑1464610

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGACTGGCCTGATAA  >  AM260479/1464303‑1464578
                                                                                                                                          |                                                                                                                                         
TGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCT                                                                                                                               >  GWNJ‑0478:712:GW2002102894th:2:1213:8705:82300/1‑150 (MQ=60)
                  ACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACC                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1204:12546:96763/150‑1 (MQ=60)
                        GTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAG                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:1212:6698:97536/148‑1 (MQ=60)
                             ataagagacagGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgc                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2211:7513:52768/138‑3 (MQ=60)
                              TGTGTCATCAGCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcat                                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1201:17187:4157/1‑146 (MQ=60)
                                        GCTGGTGGGCGAACCGTCCAGC‑GTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattgggtgagcc                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:1206:20207:46547/1‑136 (MQ=60)
                                                                 GCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCagcggctgtcgaccttgcaccggttgctcatggccgcaatgccagaggcagccagg                                                               <  GWNJ‑0478:712:GW2002102894th:2:1104:4279:6126/150‑57 (MQ=60)
                                                                        CACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTTAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCT                   >  GWNJ‑0478:712:GW2002102894th:2:1108:8823:52466/1‑150 (MQ=60)
                                                                        CACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCG                   >  GWNJ‑0478:712:GW2002102894th:2:2214:6464:6702/1‑150 (MQ=60)
                                                                                       GCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttga                                          >  GWNJ‑0478:712:GW2002102894th:2:2212:11960:5224/1‑82 (MQ=60)
                                                                                       GCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGACTGGCCTGAT    >  GWNJ‑0478:712:GW2002102894th:2:2203:5718:61484/1‑150 (MQ=60)
                                                                                         CTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgcc  >  GWNJ‑0478:712:GW2002102894th:2:1112:11214:83559/1‑96 (MQ=36)
                                                                                         CTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgcc  >  GWNJ‑0478:712:GW2002102894th:2:2106:3427:22927/1‑96 (MQ=36)
                                                                                                                                          |                                                                                                                                         
TGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGACTGGCCTGATAA  >  AM260479/1464303‑1464578

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑