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 2,907,305 2 bp→AT coding (181‑182/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,3050GA100.0% 44.1 / NA 14D61I (GAC→ATC) h16_A2693conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/6);  total (8/6)
*AM2604792,907,3060AT100.0% 41.6 / NA 14D61I (GAC→ATC) h16_A2693conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (8/6);  total (8/6)

GCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGG  >  AM260479/2907212‑2907431
                                                                                             ||                                                                                                                             
gccgccGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATg                                                                                                             <  1:764835/112‑1 (MQ=255)
gccgccGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATg                                                                                                             >  2:764835/1‑112 (MQ=255)
                                    gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCaa                                                                      <  1:747058/113‑1 (MQ=255)
                                    gCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCaa                                                                      >  2:747058/1‑113 (MQ=255)
                                                               gagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGcc                                                >  1:770886/1‑108 (MQ=255)
                                                               gagCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGcc                                                <  2:770886/108‑1 (MQ=255)
                                                                 gCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAAATATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                               >  1:303976/1‑91 (MQ=255)
                                                                 gCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAAATATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGc                                                               <  2:303976/91‑1 (MQ=255)
                                                                              agctTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCAcc       >  1:462449/1‑134 (MQ=255)
                                                                                   cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGc                                                     >  1:1030784/1‑83 (MQ=255)
                                                                                   cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGc                                                     <  2:1030784/83‑1 (MQ=255)
                                                                                   cTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAgg  >  1:33568/1‑134 (MQ=255)
                                                                                    tGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTa                                                                                        >  1:830875/1‑47 (MQ=12)
                                                                                    tGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTa                                                                                        <  2:830875/47‑1 (MQ=12)
                                                                                             ||                                                                                                                             
GCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGG  >  AM260479/2907212‑2907431

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
CGCCGCCGCGCCGCTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCC  >  AM260479/2907175‑2907442
                                                                                                                                  |                                                                                                                                         
gagacagGCGCCGCTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGC                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2115:7147:85885/143‑1 (MQ=60)
              agagacaggttcttgggcatgaTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGG                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2116:8295:23172/128‑1 (MQ=60)
                          gagatgtagctccgtctcgtgggctcggagatgtgtataagagacagGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGC                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:2213:6771:66332/103‑1 (MQ=60)
                                     cctgcgcgggcaGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAAC                                                                                <  GWNJ‑0478:712:GW2002102894th:2:2211:9769:57262/138‑1 (MQ=60)
                                                   gagacaGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGT                                                                  <  GWNJ‑0478:712:GW2002102894th:2:2114:14751:22252/144‑1 (MQ=60)
                                                                                     CCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCctgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtctt                                >  GWNJ‑0478:712:GW2002102894th:2:1116:20207:67077/1‑91 (MQ=60)
                                                                                     CCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCtgtctcttat                                 >  GWNJ‑0478:712:GW2002102894th:2:2116:16690:57338/1‑139 (MQ=60)
                                                                                       CAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAAATATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCctgtctcttatacacatctccgagcccac                              >  GWNJ‑0478:712:GW2002102894th:2:1205:4614:24433/1‑121 (MQ=60)
                                                                                                    GAGCAGTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCC                 >  GWNJ‑0478:712:GW2002102894th:2:1212:7782:74366/1‑150 (MQ=60)
                                                                                                         GTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCctgtctcttatacacatctccgagcccacgagacgg             >  GWNJ‑0478:712:GW2002102894th:2:2216:13442:43221/1‑113 (MQ=60)
                                                                                                         GTACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGC            >  GWNJ‑0478:712:GW2002102894th:2:1105:14920:35728/1‑150 (MQ=60)
                                                                                                          TACGAGCTGAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCtgtctcttatacacatctccgagcccacgagacggagctacatctcgtatgccgtcttctgcttgaaaaaa            >  GWNJ‑0478:712:GW2002102894th:2:2204:14764:68859/1‑78 (MQ=55)
                                                                                                                  GAGCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCctgtctcttatacacatctccgagcccacgagacggagctacatct   >  GWNJ‑0478:712:GW2002102894th:2:2207:13926:8641/1‑104 (MQ=60)
                                                                                                                    GCTTCTGGGAATCGATCAAGAACAGCAATTATGCCGGCGACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCtgtctcttatacacatctccgagcccacgagacggagctac  >  GWNJ‑0478:712:GW2002102894th:2:1103:11048:91678/1‑108 (MQ=60)
                                                                                                                                  |                                                                                                                                         
CGCCGCCGCGCCGCTGCCCGCGCAGGCGGCCGCCGCTGCCGCCGGCCCCGCGGCCACGGAACCGCTGCTGCTCGCCGGCATCATGCCCAAGAACACCGAGGAGCAGTACGAGCTGAGCTTCTGGGAATCGGACAAGAAACAGCAATTATGCCGGGCGAACTACGAGGCCTACCTGAAGCAGTACCCCAACGGACGCTTCGCGGTGCTGGCCAGGGCGCGCCTGGAGCGGCTGGCGGCGTCCGCCTCCGCACCCAAGGCCCAGGCGCCC  >  AM260479/2907175‑2907442

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