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
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% 45.5 / NA 13D184E (GAT→GAGh16_A1942Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (5/8);  total (5/8)

TGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGG  >  AM260479/2110885‑2111129
                                                                                                                                |                                                                                                                      
tggtggGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATggact                                                                                                                   <  1:70301/133‑5 (MQ=255)
               ccAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATc                                                                                                    >  2:910768/1‑134 (MQ=255)
                cAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATc                                                                                                    >  2:155687/1‑133 (MQ=255)
                   cATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGAcg                                                                                                <  2:39478/134‑1 (MQ=255)
                                ccACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCcgccgcgcc                                                                                      <  1:910768/131‑1 (MQ=255)
                                                        aTCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTccgccg                                                           <  2:224849/134‑1 (MQ=255)
                                                          cGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTcg                                                                                                                  >  1:662677/1‑77 (MQ=255)
                                                          cGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTcg                                                                                                                  <  2:662677/77‑1 (MQ=255)
                                                                                                                cgacgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTgcgc    <  1:120138/133‑1 (MQ=255)
                                                                                                                cgacgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTgcg     <  1:155687/132‑1 (MQ=255)
                                                                                                                 gacgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCggg                                            >  1:224484/1‑91 (MQ=255)
                                                                                                                 gacgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCggg                                            <  2:224484/92‑2 (MQ=255)
                                                                                                                   cgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCgcgg  <  1:158722/132‑1 (MQ=255)
                                                                                                                   cgaACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCgcgg  >  2:158722/1‑132 (MQ=255)
                                                                                                                                |                                                                                                                      
TGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGG  >  AM260479/2110885‑2111129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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
AGGCCGCGCGTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCGTGGCCCAGGCGGTGCGAA  >  AM260479/2110873‑2111149
                                                                                                                                             |                                                                                                                                         
aagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACG                                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2215:13580:50507/92‑1 (MQ=60)
  aagagacagGTGGTGGGGCCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAG                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1201:11582:65440/141‑1 (MQ=60)
       tattttttctttgatgtctttttttttttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATG                                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:2203:16727:72958/46‑1 (MQ=60)
          aagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagTCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGAT                                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:1112:13867:30326/91‑1 (MQ=60)
           GTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGAT                                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:1104:2677:48901/149‑1 (MQ=60)
                                           TCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCAT                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2202:3083:59445/147‑1 (MQ=60)
                                                       GTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGctgtctcttatacacatctccgagcccacgagacactcgcta                                                                             >  GWNJ‑0478:712:GW2002102894th:2:2107:21174:41692/1‑107 (MQ=60)
                                                                                    GGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaa                                               >  GWNJ‑0478:712:GW2002102894th:2:2107:9036:91319/1‑68 (MQ=60)
                                                                                                              GAGTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGgcgcggcgcgcgcgacctgtctcttatacacatctccgagccc                      >  GWNJ‑0478:712:GW2002102894th:2:1115:18286:5791/1‑106 (MQ=60)
                                                                                                                GTATGAGACCTTCGACGAACTGGTGGATGGACTCGACGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaagaaagaatgaagcacagagactagcata                     >  GWNJ‑0478:712:GW2002102894th:2:2203:8272:41975/1‑38 (MQ=56)
                                                                                                                           TCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCGTGGCCCAGGCG         <  GWNJ‑0478:712:GW2002102894th:2:1108:10790:43857/149‑1 (MQ=60)
                                                                                                                           TCGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCGTGGCCCAGGC          <  GWNJ‑0478:712:GW2002102894th:2:1111:11759:41642/148‑1 (MQ=60)
                                                                                                                            CGACGAACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCGcgcgcgaccgagg       <  GWNJ‑0478:712:GW2002102894th:2:1111:6487:57489/150‑14 (MQ=60)
                                                                                                                                 AACTGGTGGATGGACTCGACGAGGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCGTGGCCCAGGCGGTGCGgg  >  GWNJ‑0478:712:GW2002102894th:2:2111:17275:2858/1‑148 (MQ=60)
                                                                                                                                             |                                                                                                                                         
AGGCCGCGCGTGTGGTGGGACCGATCGCCAGCATGATTCCTCTTCCACAAGCGCAGTTGATCGGGCGCATCGCCAATGTGGCCGGGCGCATGCTCGCCGACGGCGCCGACGAGTATGAGACCTTCGACGAACTGGTGGATGACTCGACGAGATGCGATCGACGCCGCCGCGCCGGTGCTGGCCGGCATGGTGCTCCGCCGCGCCCTCCCCTCGGTCGCGCGCGCCGCGCCCGCCGTGCGTCGCGCCGCGGTGCGCGGCGTGGCCCAGGCGGTGCGAA  >  AM260479/2110873‑2111149

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