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 146 27.9 2327082 89.7% 2087392 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 2,907,313 Δ1 bp coding (189/1221 nt) h16_A2693 → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604792,907,3110A.100.0% 95.3 / NA 23coding (187/1221 nt)h16_A2693conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (12/11);  total (12/11)

CCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGG  >  AM260479/2907222‑2907431
                                                                                         |                                                                                                                        
ccGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTCCGAGGCCTACCt                                                                                     >  1:537598/1‑124 (MQ=255)
ccGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTCCGAGGCCTACCt                                                                                     <  2:537598/124‑1 (MQ=255)
  gcgGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTAccc                                                                         >  1:1101359/1‑134 (MQ=255)
            aaCCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                               <  1:213493/134‑1 (MQ=255)
                    ctgctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGt                                                       >  1:118273/1‑134 (MQ=255)
                       ctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGc                                                     <  2:118273/133‑1 (MQ=255)
                       ctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGc                                                     <  2:769768/133‑1 (MQ=255)
                          gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGcc                                                                                          >  1:107529/1‑93 (MQ=255)
                          gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGcc                                                                                          <  2:107529/93‑1 (MQ=255)
                          gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGc                                                 >  2:1150240/1‑134 (MQ=255)
                                atcatGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCggcg                                                                                                       <  1:22521/76‑3 (MQ=255)
                                atcatGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCggcg                                                                                                       >  2:22521/1‑74 (MQ=255)
                                                   aggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGaa                                                                                  >  2:1089421/1‑76 (MQ=255)
                                                   aggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGaa                                                                                  <  1:1089421/76‑1 (MQ=255)
                                                   aggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcg                        >  2:1026422/1‑134 (MQ=255)
                                                     gagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGACCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcggc                      >  1:512743/1‑134 (MQ=255)
                                                       gCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGt                    >  2:1050647/1‑134 (MQ=255)
                                                             cGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCt              <  1:531909/134‑1 (MQ=255)
                                                               agctgagctTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcgg                       >  1:427803/1‑123 (MQ=255)
                                                               agctgagctTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTggcgg                       <  2:427803/123‑1 (MQ=255)
                                                                         cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGcc                                                                                          >  2:804116/1‑46 (MQ=12)
                                                                         cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGcc                                                                                          <  1:804116/46‑1 (MQ=12)
                                                                            ggAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAgg  <  1:1150240/131‑1 (MQ=255)
                                                                                         |                                                                                                                        
CCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGG  >  AM260479/2907222‑2907431

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
CGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCGCTCCG  >  AM260479/2907195‑2907461
                                                                                                                      |                                                                                                                                                    
attagcgagtgtctcgtgggctcggagatgtgtataagagacaGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACC                                                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:1103:9260:45115/107‑1 (MQ=60)
            tttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTAC                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:2112:15702:74081/76‑1 (MQ=51)
            ggcatgatGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTCCGAGGCCTACCTGAAGCAGTACC                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1215:9364:82882/9‑150 (MQ=60)
              GCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCC                                                                                                      >  GWNJ‑0478:712:GW2002102894th:2:2214:7113:6629/1‑150 (MQ=60)
                   attagcgagtgtctcgtgggctcggagatgtgtataagagacagTGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGG                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2213:10909:48511/106‑1 (MQ=60)
                                GCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTG                                                                                    >  GWNJ‑0478:712:GW2002102894th:2:1108:11215:17402/1‑150 (MQ=60)
                                      GAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCtgtctcttatacacatctccgagccc                                                                              >  GWNJ‑0478:712:GW2002102894th:2:1106:17476:55529/1‑124 (MQ=60)
                                      GAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCC                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1112:11134:27700/150‑1 (MQ=60)
                                               CTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCctgtctcttatacacatctccgagcccacgagacactcgct                                                                     >  GWNJ‑0478:712:GW2002102894th:2:1110:13201:54193/1‑109 (MQ=60)
                                                                 CCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGACCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCG                                                   >  GWNJ‑0478:712:GW2002102894th:2:1212:10112:65071/1‑150 (MQ=60)
                                                                           CCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCG                                         >  GWNJ‑0478:712:GW2002102894th:2:1209:15241:33554/1‑150 (MQ=60)
                                                                                       ACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCC                             <  GWNJ‑0478:712:GW2002102894th:2:1215:6895:56878/150‑1 (MQ=60)
                                                                                          AGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGggcggcggtgcgggcgcctgggccttgggctgtc                          >  GWNJ‑0478:712:GW2002102894th:2:2112:18627:37719/1‑116 (MQ=60)
                                                                                                   TCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCtgtctcttatac                 >  GWNJ‑0478:712:GW2002102894th:2:1110:20403:39619/1‑138 (MQ=60)
                                                                                                      GGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAC                 <  GWNJ‑0478:712:GW2002102894th:2:2216:5885:38211/147‑1 (MQ=60)
                                                                                                             cgatcaagAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCctgtctcttatacacatctcc        >  GWNJ‑0478:712:GW2002102894th:2:1108:18808:42610/9‑129 (MQ=60)
                                                                                                                    GAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaa  >  GWNJ‑0478:712:GW2002102894th:2:2112:4955:33687/1‑78 (MQ=60)
                                                                                                                      |                                                                                                                                                    
CGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCGCTCCG  >  AM260479/2907195‑2907461

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: