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 3,247,959 2 bp→CA coding (723‑724/813 nt) h16_A3001 → transcriptional regulator, LysR‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,247,9590AC100.0% 91.2 / NA 24G241G (GGA→GGCh16_A3001transcriptional regulator, LysR‑family
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (13/11);  total (13/11)
*AM2604793,247,9600TA100.0% 89.0 / NA 24F242I (TTC→ATC) h16_A3001transcriptional regulator, LysR‑family
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (13/11);  total (13/11)

CTTTGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCT  >  AM260479/3247836‑3248087
                                                                                                                           ||                                                                                                                                
cTTTGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGc                                                                                                                           >  2:936713/1‑132 (MQ=255)
cTTTGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGc                                                                                                                           <  1:936713/132‑1 (MQ=255)
  ttGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGc                                                                                                                       >  2:414447/1‑134 (MQ=255)
   tGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCgg                                                                                                                        <  1:359754/132‑1 (MQ=255)
   tGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCgg                                                                                                                        <  2:25822/132‑1 (MQ=255)
                                          ccACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCa                                                                           >  1:812775/1‑134 (MQ=255)
                                                                        cATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAAcg                                             >  2:212241/1‑133 (MQ=255)
                                                                          tCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAAcgt                                            >  2:542547/1‑131 (MQ=255)
                                                                                  gACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCt                                  <  1:212241/134‑1 (MQ=255)
                                                                                                 tacgtGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAggg                   >  1:344636/1‑134 (MQ=255)
                                                                                                   cgtGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTg                 >  2:283091/1‑134 (MQ=255)
                                                                                                      ggAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCAcc                                                            <  2:45539/89‑1 (MQ=255)
                                                                                                      ggAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCAcc                                                            >  1:45539/1‑89 (MQ=255)
                                                                                                        agaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTAc            <  1:414447/134‑1 (MQ=255)
                                                                                                         gaTCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCTGTCAAGCAGGGTGCCTAcc           <  1:542547/134‑1 (MQ=255)
                                                                                                              ggCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTgcgc      <  2:344636/134‑1 (MQ=255)
                                                                                                               gCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGc                               >  2:99224/1‑108 (MQ=255)
                                                                                                               gCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGc                               <  1:99224/108‑1 (MQ=255)
                                                                                                                cTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGc                >  2:201208/1‑122 (MQ=255)
                                                                                                                cTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGc                <  1:201208/122‑1 (MQ=255)
                                                                                                                cTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAg    >  1:36601/1‑134 (MQ=255)
                                                                                                                  gggCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGc                <  2:270530/120‑1 (MQ=255)
                                                                                                                  gggCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGc                >  1:270530/1‑120 (MQ=255)
                                                                                                                  gggCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCt  >  2:78555/1‑134 (MQ=255)
                                                                                                                           ||                                                                                                                                
CTTTGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCT  >  AM260479/3247836‑3248087

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 38 ≤ 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
CAACTTTGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCG  >  AM260479/3247833‑3248105
                                                                                                                              |                                                                                                                                                   
CAACTTTGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGA                                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:2210:12493:28069/150‑1 (MQ=60)
     TTGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGC                                                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:1206:5890:84906/148‑1 (MQ=60)
                              GGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCgagcgcgaccgcaacctgcgcggcat                                                                                                >  GWNJ‑0478:712:GW2002102894th:2:2202:14296:44154/1‑124 (MQ=60)
                                                                                    GGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGG                                     <  GWNJ‑0478:712:GW2002102894th:2:1115:18786:40220/150‑1 (MQ=60)
                                                                                     GACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGT                                    >  GWNJ‑0478:712:GW2002102894th:2:1207:11502:90060/1‑150 (MQ=60)
                                                                                        ataagagacagGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGC                                  <  GWNJ‑0478:712:GW2002102894th:2:1106:6908:74389/138‑1 (MQ=60)
                                                                                          GATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACgctgtctcttatacacatctccgagcccacg                                 >  GWNJ‑0478:712:GW2002102894th:2:1102:3360:14516/1‑118 (MQ=60)
                                                                                                    TACGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGC                     >  GWNJ‑0478:712:GW2002102894th:2:1105:9069:55427/1‑150 (MQ=60)
                                                                                                      CGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTT                   <  GWNJ‑0478:712:GW2002102894th:2:1112:11740:79660/150‑1 (MQ=60)
                                                                                                      CGTGGAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTT                   >  GWNJ‑0478:712:GW2002102894th:2:1201:1934:90663/1‑150 (MQ=60)
                                                                                                          GAGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTC               <  GWNJ‑0478:712:GW2002102894th:2:1210:9864:9194/150‑1 (MQ=60)
                                                                                                           AGATCGGGCTGGGCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCTGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCT              <  GWNJ‑0478:712:GW2002102894th:2:2103:17909:39301/150‑1 (MQ=60)
                                                                                                                       GCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCG  <  GWNJ‑0478:712:GW2002102894th:2:1203:11249:81056/150‑1 (MQ=60)
                                                                                                                       GCGTGGGCATCGTCGCCGGCGTCGCCTATGATGCCGAGCGCGACCGCAACCTGCGCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACGTGACGCACCTGGCCGTCAAGCAGGGTGCCTACCTGCGCAGCTTTGTCTACACCTTCATCG  <  GWNJ‑0478:712:GW2002102894th:2:1208:17714:66442/150‑1 (MQ=60)
                                                                                                                              |                                                                                                                                                   
CAACTTTGCCGGCCGGACCAAGATCGACAAGGCCTTCGCGCTGCGCCACCTGGAGCCGGACATCGTGCTGGAGGCCATCGACGCGGACGTGATCAAGACGTACGTGGAGATCGGGCTGGGCGTGGGATTCGTCGCCGGCGTCGCCTATGATGCCCGAGCCGCGACCCGCAACCTTGCCCGGCATCTCGGCCGGCCACCTGTTTGGCAGCAACCGTGACGCACCTGGCCGTCAGCAGGGTGCCTACCTGCGCAGCTTTGTGTACACCTTCATCG  >  AM260479/3247833‑3248105

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: