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,464,365 Δ1 bp coding (49/552 nt) h16_A1354 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,464,3650A.100.0% 123.2 / NA 28coding (49/552 nt)h16_A1354conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (16/12);  total (16/12)

GCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGC  >  AM260479/1464253‑1464492
                                                                                                                |                                                                                                                               
gCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCAc                                                                                                               >  1:403139/1‑130 (MQ=255)
gCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCAc                                                                                                               <  2:403139/130‑1 (MQ=255)
              acaATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCAtt                                                                                             >  1:808880/1‑134 (MQ=255)
                      tGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCAcc                                                                                                              >  2:758223/1‑109 (MQ=255)
                      tGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCAcc                                                                                                              <  1:758223/109‑1 (MQ=255)
                      tGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGc                                                                                                <  1:18480/123‑1 (MQ=255)
                      tGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGc                                                                                                >  2:18480/1‑123 (MQ=255)
                           gcgcCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGc                                                                                 >  1:48122/1‑133 (MQ=255)
                                  cgcgCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTg                                                                         >  1:779017/1‑134 (MQ=255)
                                      cAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAg                                                                     >  1:131608/1‑134 (MQ=255)
                                                       ggAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcg                                                    >  2:781416‑M1/1‑132 (MQ=255)
                                                                 tGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGc                                                                                                <  2:472820/80‑1 (MQ=255)
                                                                 tGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGc                                                                                                >  1:472820/1‑80 (MQ=255)
                                                                 tGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggc                                          <  2:48122‑M1/134‑13 (MQ=255)
                                                                     cGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcg                                      <  2:779017‑M1/134‑17 (MQ=255)
                                                                               cTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCt                                                                                                             <  2:590927/53‑1 (MQ=255)
                                                                               cTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCt                                                                                                             >  1:590927/1‑53 (MQ=255)
                                                                                    tcatcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCAc                                                                                                               <  2:658494/46‑1 (MQ=39)
                                                                                    tcatcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCAc                                                                                                               >  1:658494/1‑46 (MQ=255)
                                                                                    tcatcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtac                        >  1:491925‑M1/1‑103 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGccgct                                                        <  1:586684/99‑1 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGccgct                                                        >  2:586684/1‑99 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccag                     <  2:146716‑M1/134‑34 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtac                        >  2:629468‑M1/1‑101 (MQ=255)
                                                                                          gCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagca                 <  2:808880‑M1/134‑38 (MQ=255)
                                                                                           cTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcac                >  2:685341‑M1/1‑96 (MQ=255)
                                                                                                 ggCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcatt          <  1:629468‑M1/134‑45 (MQ=255)
                                                                                                         gTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagc  >  2:660099‑M1/1‑82 (MQ=255)
                                                                                                                |                                                                                                                               
GCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGC  >  AM260479/1464253‑1464492

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGCGATTTCGTGGGCGCGGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCT  >  AM260479/1464235‑1464499
                                                                                                                                  |                                                                                                                                      
cgtgggctcggagatgtgtataagagacagGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCAC                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2208:20257:22564/118‑1 (MQ=60)
   GATTTCGTGGGCGCGGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGC                                                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1211:1942:70636/1‑150 (MQ=60)
              ataagagacaGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCA                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:2116:6256:49226/140‑1 (MQ=60)
                 GGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTG                                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2202:12271:89305/1‑150 (MQ=60)
                            CAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGC                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1103:8773:30180/150‑1 (MQ=60)
                              GCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCA                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:1102:6841:28194/1‑149 (MQ=60)
                                     GGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGC                                                                               >  GWNJ‑0478:712:GW2002102894th:2:2201:13593:89294/1‑150 (MQ=60)
                                         gcacccattGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGG                                                                           >  GWNJ‑0478:712:GW2002102894th:2:1109:7455:69071/10‑150 (MQ=60)
                                                                    ggagatgtgtataagagacaGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCT                                                <  GWNJ‑0478:712:GW2002102894th:2:2107:11581:2949/130‑1 (MQ=60)
                                                                    TGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCctgtctcttatacacatctccgagcccacgagacggagct                                                >  GWNJ‑0478:712:GW2002102894th:2:1215:20609:62882/1‑110 (MQ=60)
                                                                                  ATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttg                                   >  GWNJ‑0478:712:GW2002102894th:2:2107:9074:30913/1‑83 (MQ=60)
                                                                                       CGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaa                              >  GWNJ‑0478:712:GW2002102894th:2:2111:16694:57471/1‑76 (MQ=60)
                                                                                       CGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACC                              >  GWNJ‑0478:712:GW2002102894th:2:1214:10458:68008/1‑149 (MQ=60)
                                                                                         CCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaatacagtcaat                           >  GWNJ‑0478:712:GW2002102894th:2:2214:16766:93126/1‑62 (MQ=60)
                                                                                                       CATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattgggt             <  GWNJ‑0478:712:GW2002102894th:2:2115:5808:101192/150‑10 (MQ=60)
                                                                                                                  GGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattgggtgagccgccgaa  <  GWNJ‑0478:712:GW2002102894th:2:2108:11258:79792/150‑21 (MQ=60)
                                                                                                                                  |                                                                                                                                      
GGCGATTTCGTGGGCGCGGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCT  >  AM260479/1464235‑1464499

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