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,354,645 +C intergenic (+98/‑15) h16_A1248 → / → h16_A1249 predicted signal transduction protein containingCHASE, PAS, GGDEF and EAL domains/Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,354,6431.C100.0% 112.4 / NA 30intergenic (+96/‑17)h16_A1248/h16_A1249predicted signal transduction protein containingCHASE, PAS, GGDEF and EAL domains/Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (22/8);  total (22/8)

CAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATC  >  AM260479/1354530‑1354765
                                                                                                                  |                                                                                                                          
cAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGAt                                                                                                          >  1:820002/1‑133 (MQ=255)
   ccGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAg                                                                                                      <  2:820002/134‑1 (MQ=255)
    cGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAgg                                                                                                     >  1:797842/1‑134 (MQ=255)
          gcgcCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAgc                                                                                               >  1:691151/1‑134 (MQ=255)
          gcgcCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAgc                                                                                               >  2:160065/1‑134 (MQ=255)
                  cTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCAt                                                                                       <  2:691151/134‑1 (MQ=255)
                                      tgcTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCa                                                                                                              <  2:729332/91‑1 (MQ=255)
                                      tgcTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCa                                                                                                              >  1:729332/1‑91 (MQ=255)
                                                     aTAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTg                                                      >  1:324637/1‑132 (MQ=255)
                                                               gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                <  1:117673/80‑1 (MQ=255)
                                                               gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                                >  2:117673/1‑80 (MQ=255)
                                                               gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAAt                                                           <  2:77802/117‑1 (MQ=255)
                                                               gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAAt                                                           >  1:77802/1‑117 (MQ=255)
                                                               gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCg                                             >  2:376682/1‑131 (MQ=255)
                                                               gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCg                                             >  1:196448/1‑131 (MQ=255)
                                                                            ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACt                                                                                                                    <  1:259023/47‑1 (MQ=38)
                                                                            ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACt                                                                                                                    >  2:259023/1‑47 (MQ=38)
                                                                            ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCgg                              >  2:659813/1‑133 (MQ=255)
                                                                            ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCgg                              >  2:671453/1‑133 (MQ=255)
                                                                                      attcaaaTCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTc                     >  2:313279/5‑132 (MQ=255)
                                                                                      aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTc                     >  2:256417/1‑132 (MQ=255)
                                                                                      aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTc                     >  2:927154/1‑132 (MQ=255)
                                                                                           aaTCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTgg              >  2:943331/1‑134 (MQ=255)
                                                                                             tCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCAt                                                                                       <  1:92698/59‑1 (MQ=39)
                                                                                             tCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCAt                                                                                       >  2:92698/1‑59 (MQ=39)
                                                                                             tCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTggcg            >  2:96102/1‑134 (MQ=255)
                                                                                                tgtgGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTggcggc          >  2:190243/1‑133 (MQ=255)
                                                                                                tgtgGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTggcggc          >  2:1049453/1‑133 (MQ=255)
                                                                                                     cAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCa    >  2:889911/1‑134 (MQ=255)
                                                                                                       ggCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATc  <  1:256417/134‑1 (MQ=255)
                                                                                                                  |                                                                                                                          
CAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATC  >  AM260479/1354530‑1354765

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTCGAGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCGGACACGCCCGTCGTCCCCAGCGAG  >  AM260479/1354515‑1354791
                                                                                                                                   |                                                                                                                                                  
caagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATG                                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2109:14740:76209/83‑1 (MQ=60)
CTCGAGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATG                                                                                                                                   >  GWNJ‑0478:712:GW2002102894th:2:2115:21220:53370/1‑149 (MQ=60)
    ttttttttttttttttttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGG                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:2111:21315:76018/57‑1 (MQ=60)
    tttttttcaagcagaagacggcatacgagattagcgagtgtctcgtggctcggagatgtgtataagagacaGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGG                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1114:17672:54463/78‑1 (MQ=60)
    AGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGC                                                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:2112:10688:41594/1‑150 (MQ=60)
          GACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCT                                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:2106:7499:68975/1‑150 (MQ=60)
                       gcgagtgtctcgtgggctcggagatgtgtataagagacagGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAG                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:1108:20379:13816/110‑1 (MQ=60)
                                agaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCAT                                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:1107:18263:68857/90‑1 (MQ=60)
                                      CGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCctgtctcttatacacatctccgagcccac                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2109:2269:82183/1‑121 (MQ=60)
                                                     TGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGG                                                                               >  GWNJ‑0478:712:GW2002102894th:2:1202:17039:37156/1‑148 (MQ=60)
                                                               GTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGctg                                                                   >  GWNJ‑0478:712:GW2002102894th:2:1104:16656:82458/1‑147 (MQ=60)
                                                               GTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGG                                                                      >  GWNJ‑0478:712:GW2002102894th:2:1113:13065:39858/1‑147 (MQ=60)
                                                                                                                GTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCGGACACGCC                  >  GWNJ‑0478:712:GW2002102894th:2:2208:5423:76204/1‑150 (MQ=60)
                                                                                                                  GGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCctgtctcttatacacatctccgagccc                >  GWNJ‑0478:712:GW2002102894th:2:1102:13350:8758/1‑123 (MQ=60)
                                                                                                                     AGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCGGACACGCCCGTCG             <  GWNJ‑0478:712:GW2002102894th:2:1114:19123:41659/150‑1 (MQ=60)
                                                                                                                                gcCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCGGACctgtctcttatacacatctcc  >  GWNJ‑0478:712:GW2002102894th:2:2116:7944:39690/3‑129 (MQ=60)
                                                                                                                                   |                                                                                                                                                  
CTCGAGCATGGACCGCAACCGATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAAGCATCCGGACACGCCCGTCGTCCCCAGCGAG  >  AM260479/1354515‑1354791

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