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,329 Δ1 bp coding (205/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,3270G.100.0% 128.7 / NA 30coding (203/1221 nt)h16_A2693conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (15/15);  total (15/15)

CCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCG  >  AM260479/2907222‑2907456
                                                                                                         |                                                                                                                                 
ccGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTCCGAGGCCTACCt                                                                                                              <  2:537598/124‑1 (MQ=255)
ccGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTCCGAGGCCTACCt                                                                                                              >  1:537598/1‑124 (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:769768/133‑1 (MQ=255)
                       ctCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGc                                                                              <  2:118273/133‑1 (MQ=255)
                          gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGcc                                                                                                                   <  2:107529/93‑1 (MQ=255)
                          gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGcc                                                                                                                   >  1:107529/1‑93 (MQ=255)
                          gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGc                                                                          >  2:1150240/1‑134 (MQ=255)
                                atcatGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCggcg                                                                                                                                >  2:22521/1‑74 (MQ=255)
                                atcatGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCggcg                                                                                                                                <  1:22521/76‑3 (MQ=255)
                                                   aggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGaa                                                                                                           <  1:1089421/76‑1 (MQ=255)
                                                   aggagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGaa                                                                                                           >  2:1089421/1‑76 (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                                                                                                                   <  1:804116/46‑1 (MQ=12)
                                                                         cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGcc                                                                                                                   >  2:804116/1‑46 (MQ=12)
                                                                            ggAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAgg                           <  1:1150240/131‑1 (MQ=255)
                                                                                       aagaaCAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGc              <  2:1101359/131‑1 (MQ=255)
                                                                                        agaaCAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGcc                                        >  1:177368/3‑107 (MQ=255)
                                                                                        agaaCAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGcc                                        <  2:177368/105‑1 (MQ=255)
                                                                                        agaaCAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCa             >  2:198746/3‑134 (MQ=255)
                                                                                             aGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAccgcc        <  1:642558/134‑1 (MQ=255)
                                                                                             aGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCAccgcc        <  1:111672/134‑1 (MQ=255)
                                                                                                 aTTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCt                                       <  2:122857/99‑1 (MQ=255)
                                                                                                 aTTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCt                                       >  1:122857/1‑99 (MQ=255)
                                                                                                   tATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCg  >  2:964483/1‑134 (MQ=255)
                                                                                                         |                                                                                                                                 
CCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCG  >  AM260479/2907222‑2907456

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
CGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCGCTCCGGCAC  >  AM260479/2907195‑2907465
                                                                                                                                      |                                                                                                                                        
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)
                                                                                                                       CAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCGCTCCGGCAC  <  GWNJ‑0478:712:GW2002102894th:2:2105:7228:76699/150‑1 (MQ=60)
                                                                                                                       CAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCGCTCCGGCAC  <  GWNJ‑0478:712:GW2002102894th:2:1106:2895:79411/150‑1 (MQ=60)
                                                                                                                                      |                                                                                                                                        
CGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCCGCACCGCCGCCCCGCTCCGGCAC  >  AM260479/2907195‑2907465

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: