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 F1 I1 R1 154 26.6 2069536 90.2% 1866721 108.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,907,313 Δ1 bp coding (189/1221 nt) h16_A2693 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,907,3110A.100.0% 44.7 / NA 11coding (187/1221 nt)h16_A2693conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (5/6);  total (5/6)

GCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAG  >  AM260479/2907194‑2907430
                                                                                                                     |                                                                                                                       
gcgcAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCCgg                                                                                                        >  2:235580/1‑134 (MQ=255)
                                        aaCCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                              <  1:177228/134‑1 (MQ=255)
                                        aaCCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                              <  1:937494/134‑1 (MQ=255)
                                                ctgctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGt                                                      >  1:317922/1‑134 (MQ=255)
                                                                                ggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcgg                      <  2:317922/134‑1 (MQ=255)
                                                                                ggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcgg                      <  2:499535/134‑1 (MQ=255)
                                                                                   gCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGt                   >  1:103107/1‑134 (MQ=255)
                                                                                         cGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGcc              <  1:272188/133‑1 (MQ=255)
                                                                                                agctTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                              >  1:972993/1‑78 (MQ=255)
                                                                                                agctTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                              <  2:972993/78‑1 (MQ=255)
                                                                                                    tCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAg  >  2:402974/1‑134 (MQ=255)
                                                                                                                     |                                                                                                                       
GCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAG  >  AM260479/2907194‑2907430

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCGCTCCGG  >  AM260479/2907201‑2907462
                                                                                                                |                                                                                                                                                     
CGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGC                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1114:9875:80928/150‑1 (MQ=60)
CGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGC                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1205:11180:39402/150‑1 (MQ=60)
                          GCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTG                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:1204:14699:31401/1‑150 (MQ=60)
                                GAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCC                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1113:10808:45317/150‑1 (MQ=60)
                                GAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCC                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2210:2695:58858/150‑1 (MQ=60)
                                           cagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagatcaagAACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCC                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1105:20517:35267/82‑1 (MQ=60)
                                                             CAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTC                                                  >  GWNJ‑0478:712:GW2002102894th:2:1106:19576:97194/1‑150 (MQ=60)
                                                                          GAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCctgtctcttatacacatctccgagcccacgagacggagctac                                     >  GWNJ‑0478:712:GW2002102894th:2:2213:17953:61951/1‑108 (MQ=60)
                                                                               ataagagacagGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGC                                <  GWNJ‑0478:712:GW2002102894th:2:1204:19225:66548/139‑1 (MQ=60)
                                                                               GTACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaa                                 >  GWNJ‑0478:712:GW2002102894th:2:1210:18293:72434/1‑71 (MQ=42)
                                                                                 ACGAGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCC                               <  GWNJ‑0478:712:GW2002102894th:2:1203:17184:22317/149‑1 (MQ=60)
                                                                                    AGCTGAGCTTCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaa                            >  GWNJ‑0478:712:GW2002102894th:2:2202:17775:2256/1‑66 (MQ=32)
                                                                                             TCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCtgtctcttatacacatctccgagcccacgagacgg                  >  GWNJ‑0478:712:GW2002102894th:2:1203:13078:7161/1‑115 (MQ=60)
                                                                                             TCTGGGAATCGATCAAG‑AACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAC                  >  GWNJ‑0478:712:GW2002102894th:2:2210:9663:74677/1‑150 (MQ=60)
                                                                                                         atcaagAACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCctgtctcttatacacatctccgagcccacg       >  GWNJ‑0478:712:GW2002102894th:2:2110:20353:43629/7‑120 (MQ=60)
                                                                                                              GAACAGCAATTATGCC‑GGCG‑ACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCGCTCCGG  <  GWNJ‑0478:712:GW2002102894th:2:1211:9953:58655/150‑1 (MQ=60)
                                                                                                                |                                                                                                                                                     
CGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCGCTCCGG  >  AM260479/2907201‑2907462

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: ‑