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 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% 100.7 / NA 23coding (49/552 nt)h16_A1354conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (11/12);  total (11/12)

GCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCA  >  AM260479/1464253‑1464472
                                                                                                                |                                                                                                           
gCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTATGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTgg                                                                                       <  1:104221/134‑1 (MQ=255)
                      tGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCAt                                                                 >  2:996499/1‑134 (MQ=255)
                      tGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCAt                                                                 >  2:779626/1‑134 (MQ=255)
                      tGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCAt                                                                 >  1:353121/1‑134 (MQ=255)
                         gggCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAg                                                              <  1:779626/134‑1 (MQ=255)
                         gggCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAg                                                              <  1:836245/134‑1 (MQ=255)
                           gcgcCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTgccgcc                                                                                                     >  2:955912/1‑93 (MQ=255)
                           gcgcCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTgccgcc                                                                                                     <  1:955912/93‑1 (MQ=255)
                                gccgcGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAg                                                              >  2:46084/1‑127 (MQ=255)
                                gccgcGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAg                                                              <  1:46084/127‑1 (MQ=255)
                                  cgcgCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTg                                                     >  2:294448/1‑134 (MQ=255)
                                                                     cGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcg                  <  1:873641‑M1/134‑17 (MQ=255)
                                                                           tgACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgc              <  1:294448‑M1/132‑21 (MQ=255)
                                                                            gACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctg           >  1:426762‑M1/1‑111 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCgg                                                       >  2:172058/1‑80 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCgg                                                       <  1:172058/80‑1 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtacca  <  2:353121‑M1/133‑33 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggc                <  2:148835‑M1/119‑19 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggc                >  1:148835‑M1/1‑101 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgc                   >  2:394075‑M1/1‑101 (MQ=255)
                                                                                      atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgc                   <  1:394075‑M1/116‑16 (MQ=255)
                                                                                                  gCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcg                                >  1:635033‑M1/1‑89 (MQ=255)
                                                                                                  gCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcg                                <  2:635033‑M1/91‑3 (MQ=255)
                                                                                                                |                                                                                                           
GCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCA  >  AM260479/1464253‑1464472

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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
GATTTCGTGGGCGCGGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACG  >  AM260479/1464238‑1464486
                                                                                                                               |                                                                                                                         
gggctcggagatgtgtataagagacagGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGC                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2212:20012:39341/123‑1 (MQ=60)
              GGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTATGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTG                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1108:15496:7716/150‑1 (MQ=60)
                      GTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATG                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1206:2231:44297/1‑150 (MQ=60)
                                       AGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGT                                                             <  GWNJ‑0478:712:GW2002102894th:2:1102:12834:18116/150‑1 (MQ=60)
                                       AGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGT                                                             <  GWNJ‑0478:712:GW2002102894th:2:2201:4789:11389/150‑1 (MQ=60)
                                       AGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGT                                                             <  GWNJ‑0478:712:GW2002102894th:2:2205:12565:5553/150‑1 (MQ=60)
                                              agctccgtctcgtgggctcggagatgtgtataagagacaGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGC                                                      <  GWNJ‑0478:712:GW2002102894th:2:1113:10613:15277/111‑1 (MQ=60)
                                                                            GCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGG                        >  GWNJ‑0478:712:GW2002102894th:2:1210:19172:79180/1‑150 (MQ=60)
                                                                                  acaGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTAC                  <  GWNJ‑0478:712:GW2002102894th:2:1209:6475:94074/147‑1 (MQ=60)
                                                                                   ACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACC                 <  GWNJ‑0478:712:GW2002102894th:2:2206:17113:40693/150‑1 (MQ=60)
                                                                                      CCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGC              >  GWNJ‑0478:712:GW2002102894th:2:1111:20505:77502/1‑150 (MQ=60)
                                                                                         GTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAG             <  GWNJ‑0478:712:GW2002102894th:2:1202:13515:59984/148‑1 (MQ=60)
                                                                                                  GTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCtgtctcttatacacatctccgagcccac  >  GWNJ‑0478:712:GW2002102894th:2:2108:13895:70704/1‑122 (MQ=60)
                                                                                                         GCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTTGCATTGCGGCCATG                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1109:5816:81725/1‑67 (MQ=60)
                                                                                                                               |                                                                                                                         
GATTTCGTGGGCGCGGCAATGGGTGCAGCACAATGGCTGAGGGCGCCGCCGCGCAATCATCGCCGTGACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACG  >  AM260479/1464238‑1464486

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: