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 149 26.7 2337422 87.5% 2045244 105.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 2,308,233 (C)5→4 intergenic (+176/‑161) h16_B2050 → / → h16_B2051 response regulator containing a GGDEF effector domain/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604802,308,2290C.100.0% 133.0 / NA 30intergenic (+172/‑165)h16_B2050/h16_B2051response regulator containing a GGDEF effector domain/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/18);  total (12/18)

GCCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATGGG  >  AM260480/2308102‑2308353
                                                                                                                               |                                                                                                                            
gcCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGccccgcg                                                                                                                        >  1:39920/1‑131 (MQ=255)
 ccACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGccccgcgc                                                                                                                       >  2:846960/1‑130 (MQ=255)
 ccACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGccccgcgc                                                                                                                       >  1:723384/1‑130 (MQ=255)
   aCCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCt                                                                                                                    <  2:39920/134‑1 (MQ=255)
   aCCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCt                                                                                                                    <  1:391975/134‑1 (MQ=255)
   aCCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCt                                                                                                                    <  2:243739/134‑1 (MQ=255)
   aCCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCt                                                                                                                    <  1:846960/134‑1 (MQ=255)
         cAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCg                                                                                                              <  1:879092/134‑1 (MQ=255)
          aaCCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGc                                                                                                             <  2:963215/134‑1 (MQ=255)
             cGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTCGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGcccc                                                                                                          <  1:415765/134‑1 (MQ=255)
              gggCCACCAGCCCGGCCAACGGCCCGCCAGGCCCCGCGCCTACGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCTCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCt                                                                                                         <  2:34924/134‑1 (MQ=255)
                                       gCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCg                                                                                                              >  1:1134679/1‑104 (MQ=255)
                                       gCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCg                                                                                                              <  2:1134679/104‑1 (MQ=255)
                                                          ggCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCt                                                                                                                    <  1:430798/79‑1 (MQ=255)
                                                          ggCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCt                                                                                                                    >  2:430798/1‑79 (MQ=255)
                                                          ggCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTgcg                                                             <  2:66863/134‑1 (MQ=255)
                                                           gCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGt                                                                                                     <  1:151077/93‑1 (MQ=255)
                                                           gCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGt                                                                                                     >  2:151077/1‑93 (MQ=255)
                                                           gCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATc                                                                                  <  1:398179/112‑1 (MQ=255)
                                                           gCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATc                                                                                  >  2:398179/1‑112 (MQ=255)
                                                                gccggccCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGc                                                        >  2:992856/1‑133 (MQ=255)
                                                                gccggccCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGc                                                        >  1:26723/1‑133 (MQ=255)
                                                                 ccggccCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGcaa                                                      >  2:98929/1‑134 (MQ=255)
                                                                         gCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCg                                                                                                              <  2:1026943/70‑1 (MQ=255)
                                                                         gCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCg                                                                                                              >  1:1026943/1‑70 (MQ=255)
                                                                                         tCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGc                                                                                                             <  1:269675/55‑1 (MQ=255)
                                                                                         tCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGc                                                                                                             >  2:269675/1‑55 (MQ=255)
                                                                                             cGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGAcc                          <  1:98929/134‑1 (MQ=255)
                                                                                                      gAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAacaaca                                                 <  2:933773/102‑1 (MQ=255)
                                                                                                      gAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAacaaca                                                 >  1:933773/1‑102 (MQ=255)
                                                                                                      gAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGccc                 >  2:398620/1‑134 (MQ=255)
                                                                                                            tCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCcgcgc           <  2:26723/134‑1 (MQ=255)
                                                                                                                     cGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATggg  >  2:158246/1‑134 (MQ=255)
                                                                                                                               |                                                                                                                            
GCCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATGGG  >  AM260480/2308102‑2308353

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
AGTTACACACTGGCGGCCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATGGGCCAGATCGTCTCTGACTC  >  AM260480/2308087‑2308371
                                                                                                                                                  |                                                                                                                                          
AGTTACACACTGGCGGCCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGC                                                                                                                                        >  GWNJ‑0478:712:GW2002102894th:2:1105:7590:55180/1‑150 (MQ=60)
 GTTACACACTGGCGGCCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCC                                                                                                                                       >  GWNJ‑0478:712:GW2002102894th:2:2109:9852:85292/1‑150 (MQ=60)
                 agcgagtgtctcgtgggctcggagatgtgtataagagacagGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGT                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1209:12180:91861/109‑1 (MQ=60)
                 CACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGT                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:1206:21154:71554/150‑1 (MQ=60)
                 CACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGT                                                                                                                       <  GWNJ‑0478:712:GW2002102894th:2:2114:18421:95733/150‑1 (MQ=60)
                       CCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATAC                                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2201:18249:62487/150‑1 (MQ=60)
                        aagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACC                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1116:9308:27833/86‑1 (MQ=60)
                          ttttgatttgtttttgttcttttttttttttttttttcaagcagaagacgtcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATAC                                                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2207:2800:81395/44‑1 (MQ=60)
                           CCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTCGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTAT                                                                                                             <  GWNJ‑0478:712:GW2002102894th:2:1209:10363:7146/150‑1 (MQ=60)
                                gggctcggagatgtgtataagagaCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTG                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1111:5902:17730/126‑1 (MQ=60)
                                       tgcgtggcggggccgGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTgtctcttatacac                                                                                                  >  GWNJ‑0478:712:GW2002102894th:2:2215:16625:53033/16‑136 (MQ=60)
                                                   agagaCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATC                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:1208:13326:8636/145‑1 (MQ=60)
                                                     ttttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCG                                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2215:2317:97064/71‑1 (MQ=60)
                                                                GCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCA                                                                         >  GWNJ‑0478:712:GW2002102894th:2:1103:2033:72615/1‑149 (MQ=60)
                                                                         GGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTgtctcttatacacatctccgagcccacgagacactcgctaatctcgt                                                                >  GWNJ‑0478:712:GW2002102894th:2:2211:18605:27118/1‑102 (MQ=60)
                                                                                                      CGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCctgtctcttatacacatc                                  >  GWNJ‑0478:712:GW2002102894th:2:2204:3614:99235/1‑132 (MQ=60)
                                                                                                           CCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGC                             <  GWNJ‑0478:712:GW2002102894th:2:1106:18602:18107/150‑1 (MQ=60)
                                                                                                                                      ATTCGGAGCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgct  >  GWNJ‑0478:712:GW2002102894th:2:2114:9085:23467/1‑86 (MQ=60)
                                                                                                                                                  |                                                                                                                                          
AGTTACACACTGGCGGCCACCGGCCAACCGGGCCACCAGCCCGGCCAACGTCCCGCCAGGCCCCGCACCTGCGGGCCTGGCCGGCCCCGCCACGCACGCACGCGTCACCGCTTGTCGGAATTGTCCGACATGCGATTCGGAGCCCCCGCGCCTGTCTCGCCCCTGGGTTGATACCTATCGTTGCATCACGACAGGTGCAATCGCTGCGCCTGCAACAACAGACAGCCAACTGGCCACGCGACCGGCGCGCCCCCGCGCGGGCATGGGCCAGATCGTCTCTGACTC  >  AM260480/2308087‑2308371

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