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 152 24.9 2083062 91.6% 1908084 109.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,111,013 +G coding (553/1089 nt) h16_A1942 → Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,111,0121.G100.0% 39.4 / NA 12D184E (GAT→GAGh16_A1942Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (5/7);  total (5/7)

GGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCC  >  AM260479/2110890‑2111097
                                                                                                                           |                                                                                       
ggACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTcgacga                                                                              <  2:442531/134‑1 (MQ=255)
      aTCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACgaggatgc                                                                       >  2:194416/1‑129 (MQ=255)
                               aaGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTg                                               <  1:376689/134‑1 (MQ=255)
                                  cgcAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCcgccgcg                                                        <  1:291255/122‑1 (MQ=255)
                                  cgcAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCcgccgcg                                                        >  2:291255/1‑122 (MQ=255)
                                             gggCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGc                                  <  2:295076/133‑1 (MQ=255)
                                                                    gcgcATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATCGCACTCCGCCCCGCCCCCCCCTCGGTcg          >  2:758152/1‑134 (MQ=255)
                                                                              gCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGcgcgcc  <  1:419459/132‑1 (MQ=255)
                                                                                  aCGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGAcg                                                                <  1:110601/66‑1 (MQ=255)
                                                                                  aCGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGAcg                                                                >  2:110601/1‑66 (MQ=255)
                                                                                  aCGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCAt                                       >  1:391197/1‑91 (MQ=255)
                                                                                  aCGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCAt                                       <  2:391197/91‑1 (MQ=255)
                                                                                                                           |                                                                                       
GGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCC  >  AM260479/2110890‑2111097

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ 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
GCCGCGCGTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCGTGGCCCAGGCGGTGCGAACGGCCACCCGC  >  AM260479/2110875‑2111160
                                                                                                                                           |                                                                                                                                                    
GCCGCGCGTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAG                                                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:1210:11214:73206/150‑1 (MQ=60)
                            ggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGC                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1108:9301:47007/97‑1 (MQ=60)
                                    CCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCctg                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1202:6190:37919/150‑4 (MQ=60)
                                     CTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGctgtctcttatacacatctccgagcccacgagacggagc                                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:1114:16281:38239/1‑111 (MQ=60)
                                             CAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTC                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1208:13310:88168/150‑1 (MQ=60)
                                                                                  GGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCtgtctcttatacacatctccgagcccac                                                          >  GWNJ‑0478:712:GW2002102894th:2:1209:3197:77544/1‑122 (MQ=60)
                                                                                            CGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGT                                                  <  GWNJ‑0478:712:GW2002102894th:2:1210:14987:37195/148‑1 (MQ=60)
                                                                                                    gggctcggagatgtgtataagagacaGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGG                                        <  GWNJ‑0478:712:GW2002102894th:2:1216:17141:32674/124‑1 (MQ=60)
                                                                                                              GTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGgactgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaagagaagaaaaaagtaaagaattaac                             >  GWNJ‑0478:712:GW2002102894th:2:2211:14948:70287/1‑40 (MQ=60)
                                                                                                                                        atgGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCGTGGCCCAGGCGGTGCGAACGGCCACCC    <  GWNJ‑0478:712:GW2002102894th:2:1112:16412:42165/147‑1 (MQ=60)
                                                                                                                                         GGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCGTGGCCCAGGCGGTGCGAACGGCCACCCGC  >  GWNJ‑0478:712:GW2002102894th:2:2215:19021:37473/1‑150 (MQ=60)
                                                                                                                                           |                                                                                                                                                    
GCCGCGCGTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCGTGGCCCAGGCGGTGCGAACGGCCACCCGC  >  AM260479/2110875‑2111160

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