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 1,071,935 +G coding (626/627 nt) h16_A0979 → Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,071,9341.G100.0% 81.6 / NA 22*209G (TGA→GGA) h16_A0979Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (13/9);  total (13/9)

GAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGG  >  AM260479/1071829‑1072061
                                                                                                          |                                                                                                                                
gAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCt                                                                                                       <  1:558831/134‑1 (MQ=255)
     ggggTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCg                                                                                                  >  2:798055/1‑134 (MQ=255)
                 cGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAg                                                                                      >  1:237370/1‑134 (MQ=255)
                     ggCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGt                                                                                  >  1:307800/1‑134 (MQ=255)
                                           cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGt                                                             >  1:874712/1‑133 (MQ=255)
                                           cgcTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGt                                                             >  2:1119632/1‑133 (MQ=255)
                                             cTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGAt                                                          <  1:798055/134‑1 (MQ=255)
                                                tgGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGa                                                                                                                         >  2:589921/1‑68 (MQ=255)
                                                tgGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGa                                                                                                                         <  1:589921/68‑1 (MQ=255)
                                                     gggCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACcg                                                  <  2:307800/134‑1 (MQ=255)
                                                            aCCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGAc                                           >  2:175958/1‑134 (MQ=255)
                                                                               aCCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAggg                                                                  >  2:352192/1‑90 (MQ=255)
                                                                               aCCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAggg                                                                  <  1:352192/91‑2 (MQ=255)
                                                                                   cacCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCTGGc                    <  2:237370/134‑1 (MQ=255)
                                                                                   cacCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGc                    <  1:175958/134‑1 (MQ=255)
                                                                                       gTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCgg                 >  2:62068/1‑133 (MQ=255)
                                                                                       gTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCgg                 >  2:784803/1‑133 (MQ=255)
                                                                                          ctgctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAAt             >  2:246108/1‑134 (MQ=255)
                                                                                            gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAg                                                                                      <  2:214792/59‑1 (MQ=39)
                                                                                            gctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAg                                                                                      >  1:214792/1‑59 (MQ=39)
                                                                                             ctgACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCt          >  2:586563/1‑134 (MQ=255)
                                                                                                     cGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTgg  <  1:784803/134‑1 (MQ=255)
                                                                                                          |                                                                                                                                
GAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGG  >  AM260479/1071829‑1072061

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
GCGCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGCCGCATCC  >  AM260479/1071810‑1072076
                                                                                                                              |                                                                                                                                              
catacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCAC                                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:2101:10403:45089/99‑1 (MQ=60)
                  CGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGA                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1214:18356:82182/150‑1 (MQ=60)
                     GGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGC                                                                                                    >  GWNJ‑0478:712:GW2002102894th:2:1115:20535:46368/1‑150 (MQ=60)
                         GGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTG                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:1203:18332:33070/1‑150 (MQ=60)
                                               ACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTG                                                                           >  GWNJ‑0478:712:GW2002102894th:2:2201:21063:32619/1‑149 (MQ=60)
                                                      gtgggctcggagatgtgtataagagacaGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACAC                                                                   <  GWNJ‑0478:712:GW2002102894th:2:1204:18791:3709/122‑1 (MQ=60)
                                                               GCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGAC                                                          <  GWNJ‑0478:712:GW2002102894th:2:2112:17825:42518/150‑1 (MQ=60)
                                                                                                GGACCACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCActgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttct                         >  GWNJ‑0478:712:GW2002102894th:2:1112:12587:34001/1‑89 (MQ=60)
                                                                                                     ACACCGTGCTGCTGACCAACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAAC                    <  GWNJ‑0478:712:GW2002102894th:2:1110:3087:32237/150‑1 (MQ=60)
                                                                                                                       ACGCCTGGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGCCGCATCC  <  GWNJ‑0478:712:GW2002102894th:2:2113:2666:73328/150‑1 (MQ=60)
                                                                                                                              |                                                                                                                                              
GCGCAACCTCTACCGGGCCGAGGTGGGGTTGTCCGTCGCCGGCGACAACGCGGCCCTGCGCGCGCTGTGGCCGGGCTGGACCGGCTTCGGCATGCAGGACCACACCGTGCTGCTGACCAACGCCTGACGCCCGAGGGGCGCAACGGCACGCTCGCCGTGATCAAGGAAGCGGTGCCCACCAGCAAGGCGGCGTGATCGACACCGCGCTGACCGTCGGGCTGGCCGCCTTCGGGCCGGAAATCCTGAAACTGGACCTGGGCCGCATCC  >  AM260479/1071810‑1072076

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: