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 154 26.6 2069536 90.2% 1866721 108.9

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% 91.8 / NA 24L40P (CTG→CCG) h16_A1354conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/12);  total (12/12)

CATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTG  >  AM260479/1464316‑1464526
                                                                                                                       |                                                                                           
cATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgg                                                                              <  2:304309‑M1/134‑11 (MQ=255)
cATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgg                                                                              <  2:390557‑M1/134‑11 (MQ=255)
            tgACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcct                                                                  >  2:957442‑M1/1‑112 (MQ=255)
                cTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                              >  1:89234‑M1/1‑108 (MQ=255)
                       atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccag                                                       <  1:901836‑M1/134‑34 (MQ=255)
                       atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccag                                                       <  1:957442‑M1/134‑34 (MQ=255)
                       atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtacca                                                        <  2:89234‑M1/133‑33 (MQ=255)
                       atcaGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtac                                                          >  1:838934‑M1/1‑101 (MQ=255)
                            cTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcac                                                  >  1:113229‑M1/1‑96 (MQ=255)
                                   gCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGTCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggct                                                                           <  2:986947‑M1/102‑14 (MQ=255)
                                   gCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGTCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggct                                                                           >  1:986947‑M1/1‑89 (MQ=255)
                                               agcgtgcCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgcc                                <  2:986332‑M1/131‑57 (MQ=255)
                                                         cacaCGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtgg                      >  2:805119‑M1/1‑68 (MQ=255)
                                                          acacGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggc                     >  2:331366‑M1/1‑67 (MQ=255)
                                                          acacGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggc                     <  2:566912‑M1/134‑68 (MQ=255)
                                                                          gCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgccc     >  2:664172‑M1/1‑51 (MQ=255)
                                                                                  cATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcc                                                                   >  2:696887‑M1/1‑43 (MQ=255)
                                                                                  cATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcc                                                                   <  1:696887‑M1/64‑22 (MQ=255)
                                                                                           cATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggt  >  2:759275‑M1/1‑34 (MQ=255)
                                                                                           cATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggt  <  1:759275‑M1/120‑87 (MQ=255)
                                                                                                           aaGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccgg   <  2:239491‑M1/103‑86 (MQ=255)
                                                                                                           aaGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccgg   >  1:239491‑M1/1‑18 (MQ=255)
                                                                                                            aGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                              >  2:425306‑M1/1‑17 (MQ=255)
                                                                                                            aGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                              <  1:425306‑M1/43‑27 (MQ=255)
                                                                                                                       |                                                                                           
CATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTG  >  AM260479/1464316‑1464526

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
GACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGAC  >  AM260479/1464304‑1464567
                                                                                                                                   |                                                                                                                                    
gagatgtgtataagagacaGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTG                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2101:7469:63797/131‑1 (MQ=60)
             ATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATG                                                                                                      >  GWNJ‑0478:712:GW2002102894th:2:1106:1634:8122/1‑150 (MQ=60)
                    CCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACC                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2205:18402:48831/1‑147 (MQ=60)
                       acggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagCCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGC                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:2113:18604:62230/94‑1 (MQ=60)
                         GACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACG                                                                                          >  GWNJ‑0478:712:GW2002102894th:2:1108:7005:61705/1‑150 (MQ=60)
                             ttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagTGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgc                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1210:6210:86554/73‑3 (MQ=38)
                                GTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGTCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTctgtctcttatacacatc                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2214:3563:38645/1‑132 (MQ=60)
                                  CATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattgggt                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2208:11971:45971/150‑10 (MQ=60)
                                  CATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGgcattgggt                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2212:11302:74673/150‑10 (MQ=60)
                                           GTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTctgtctcttatacacatctccgagcccacgagacggagctacat                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1109:14344:50570/1‑106 (MQ=60)
                                                                                                                   GTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaaaagaaaa  >  GWNJ‑0478:712:GW2002102894th:2:1207:18615:34373/1‑58 (MQ=16)
                                                                                                                                   |                                                                                                                                    
GACAGGAGATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGAC  >  AM260479/1464304‑1464567

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