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 150 26.0 2053212 90.0% 1847890 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,056,947 +G intergenic (+211/‑15) h16_A0964 → / → h16_A0965 conserved hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,056,9451.G100.0% 51.2 / NA 15intergenic (+209/‑17)h16_A0964/h16_A0965conserved hypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (8/7);  total (8/7)

GCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGG  >  AM260479/1056825‑1057043
                                                                                                                         |                                                                                                  
gCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGccgcccgc                                                                                          >  2:212106/1‑132 (MQ=255)
                   cAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCcgcg                                                                            <  1:319745/127‑1 (MQ=255)
                   cAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCcgcg                                                                            >  2:319745/1‑127 (MQ=255)
                      ccgccgGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCcgcg                                                                            <  1:46002/124‑1 (MQ=255)
                      ccgccgGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCcgcg                                                                            >  2:46002/1‑124 (MQ=255)
                                                    tggCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGa                                                                <  1:858935/106‑1 (MQ=255)
                                                    tggCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGa                                                                >  2:858935/1‑106 (MQ=255)
                                                    tggCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGAcccc                                    >  2:807211/1‑134 (MQ=255)
                                                                           cAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTc                                 <  1:716358/114‑1 (MQ=255)
                                                                           cAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTc                                 >  2:716358/1‑114 (MQ=255)
                                                                                      cAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCggg  <  1:807211/134‑1 (MQ=255)
                                                                                                       cccGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                <  1:170842/39‑1 (MQ=25)
                                                                                                       cccGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGc                                                                                >  2:170842/1‑39 (MQ=25)
                                                                                                        ccGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCa                                                                   <  1:857565/51‑1 (MQ=38)
                                                                                                        ccGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCa                                                                   >  2:857565/1‑51 (MQ=38)
                                                                                                                         |                                                                                                  
GCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGG  >  AM260479/1056825‑1057043

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ 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
GTGCATTAGCTTCTGTTGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTACGGCCGCATGATCCAG  >  AM260479/1056808‑1057084
                                                                                                                                            |                                                                                                                                         
GTGCATTAGCTTCTGTTGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCC                                                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2107:6871:27098/1‑150 (MQ=60)
                cgtatatttttgctttttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGG                                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:2205:5223:57042/59‑1 (MQ=60)
                         tattttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACG                                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:1113:9057:6589/69‑1 (MQ=60)
                            ACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGC                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1102:8162:13945/150‑1 (MQ=60)
                            ACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGC                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1204:11904:54979/150‑1 (MQ=60)
                                      ttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAA                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:2209:2624:100476/81‑1 (MQ=60)
                                         gtataagagacaGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGT                                                                                          <  GWNJ‑0478:712:GW2002102894th:2:2206:5529:9137/137‑1 (MQ=60)
                                                                       agacagCGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATC                                                           <  GWNJ‑0478:712:GW2002102894th:2:2115:16167:3927/144‑1 (MQ=60)
                                                                                                      CCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGC                            <  GWNJ‑0478:712:GW2002102894th:2:2202:16913:78265/150‑1 (MQ=60)
                                                                                                                                gacagcggcagGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTACGGCCGCATGATCCAG  <  GWNJ‑0478:712:GW2002102894th:2:1106:8125:60538/139‑1 (MQ=60)
                                                                                                                                            |                                                                                                                                         
GTGCATTAGCTTCTGTTGCAGATGCTGCACAAATGGCAACCGCCGGCGGCGCCGGCATTGGGGGCATGGTGGCGACACGTCAGGAGGGCGGGCAAGCCGTGGCCAAGGGTACTGTCGAGGCCCGGCGCGAGCGCGGCAGGCCGCCCGCCAGCAGGTGCCGCGCGGTGCGCACGAGGCCATCGGCAACGTCGACCGCGACCCCGTCGAACTGCTGGAAATCAGCAGCCAGGGCCGGGTTGAGCGGCTGGTGCCGCTGCGCTACGGCCGCATGATCCAG  >  AM260479/1056808‑1057084

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