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 155 25.8 2285122 88.0% 2010907 104.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,464,435 T→C L40P (CTG→CCG)  h16_A1354 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,464,4350TC100.0% 119.8 / NA 32L40P (CTG→CCG) h16_A1354conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/17);  total (15/17)

TCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGC  >  AM260479/1464337‑1464552
                                                                                                  |                                                                                                                     
tcatcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtac                                                                                    >  2:307527‑M1/1‑103 (MQ=255)
  atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccag                                                                                 <  2:96625‑M1/134‑34 (MQ=255)
  atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtac                                                                                    >  1:983639‑M1/1‑101 (MQ=255)
  atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtac                                                                                    >  1:96625‑M1/1‑101 (MQ=255)
  atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                        <  1:541166‑M1/127‑27 (MQ=255)
  atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                        >  2:810812‑M1/1‑101 (MQ=255)
  atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                        >  2:541166‑M1/1‑101 (MQ=255)
  atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                        <  1:810812‑M1/127‑27 (MQ=255)
     aGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagag                                                                               <  1:641957‑M1/133‑36 (MQ=255)
            gggCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtacca                                                                                  >  2:1118729‑M1/1‑91 (MQ=255)
            gggCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtacca                                                                                  <  1:1118729‑M1/123‑33 (MQ=255)
              gCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgt                                                                                      <  1:680041‑M1/117‑29 (MQ=255)
              gCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgt                                                                                      >  2:680041‑M1/1‑89 (MQ=255)
                gAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattggg                                                                   <  1:307527‑M1/134‑48 (MQ=255)
                          agcgtgcCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgcc                                                          >  2:133270‑M1/4‑78 (MQ=255)
                                ccgccACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaata                                                     <  1:927146‑M1/133‑62 (MQ=255)
                                ccgccACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaata                                                     <  1:133270‑M1/133‑62 (MQ=255)
                                     acacGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtacca                                                                                  >  2:105335‑M1/1‑67 (MQ=255)
                                     acacGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtacca                                                                                  <  1:105335‑M1/99‑33 (MQ=255)
                                                            gCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtc                           <  1:747949‑M1/131‑88 (MQ=255)
                                                                     ccATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgat               <  2:983639‑M1/134‑100 (MQ=255)
                                                                       aTGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataa             <  1:283077‑M1/134‑102 (MQ=255)
                                                                              aCCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgcgccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggc      >  2:206267‑M1/1‑26 (MQ=255)
                                                                                cGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGgcgc                                                                                              <  2:635517/44‑1 (MQ=18)
                                                                                cGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGgcgc                                                                                              >  1:635517/1‑44 (MQ=18)
                                                                                   tGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggctgag  >  1:841749‑M1/1‑21 (MQ=255)
                                                                                   tGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                        <  1:923232‑M1/47‑27 (MQ=255)
                                                                                   tGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                        >  2:923232‑M1/1‑21 (MQ=255)
                                                                                       aGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtacca                                                                                  >  2:62543‑M1/1‑17 (MQ=255)
                                                                                       aGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtacca                                                                                  <  1:62543‑M1/49‑33 (MQ=255)
                                                                                        ggTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccgg                             <  2:1115350‑M1/101‑86 (MQ=255)
                                                                                        ggTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccgg                             >  1:1115350‑M1/1‑16 (MQ=255)
                                                                                                  |                                                                                                                     
TCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGC  >  AM260479/1464337‑1464552

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 35 ≤ 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
CGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGG  >  AM260479/1464298‑1464550
                                                                                                                                         |                                                                                                                   
ggctcggagatgtgtataagagacaGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCAC                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1112:17019:29735/125‑1 (MQ=60)
   CGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGG                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1208:15633:46764/150‑1 (MQ=60)
                          CCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACC                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1107:3148:100132/1‑147 (MQ=60)
                          CCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACC                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:2207:11964:82580/1‑147 (MQ=60)
                                 CCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgc                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1214:7686:32376/150‑3 (MQ=60)
                                 CCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgc                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2115:1904:75520/150‑3 (MQ=60)
                                   ttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgc                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2205:20841:32750/78‑3 (MQ=49)
                                   cagGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcat                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2106:5970:69511/147‑5 (MQ=60)
                                        ggagatgtgtataagagacaggtccagcGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattggg                                                                 <  GWNJ‑0478:712:GW2002102894th:2:1106:6734:51833/122‑9 (MQ=60)
                                        CATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattggg                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2215:15980:84991/149‑9 (MQ=60)
                                         ttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagTGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattggg                                                                 <  GWNJ‑0478:712:GW2002102894th:2:1102:18946:97471/79‑9 (MQ=23)
                                           CAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattgggtga                                                              <  GWNJ‑0478:712:GW2002102894th:2:2104:5257:36844/149‑12 (MQ=60)
                                                      CGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattgggtgagccgccgaatat                                                  <  GWNJ‑0478:712:GW2002102894th:2:1203:18477:59194/150‑24 (MQ=60)
                                                                      GCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattgggtgagccgccgaatatgtggcgaatgggcc                                    <  GWNJ‑0478:712:GW2002102894th:2:1109:8981:13057/149‑38 (MQ=60)
                                                                      GCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattgggtgagccgccgaatatgtggcgaatgggcc                                    <  GWNJ‑0478:712:GW2002102894th:2:2205:3700:55691/149‑38 (MQ=60)
                                                                                                        GCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCtgtctcttatacacatcccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaa  >  GWNJ‑0478:712:GW2002102894th:2:2105:16141:96930/1‑75 (MQ=59)
                                                                                                                                         |                                                                                                                   
CGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGG  >  AM260479/1464298‑1464550

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: