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 F2 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,849,816 T→C E117G (GAA→GGA)  ygbA ← nitrous oxide‑stimulated promoter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,849,8160TC100.0% 85.4 / NA 25E117G (GAA→GGA) ygbAnitrous oxide‑stimulated promoter family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/15);  total (10/15)

TCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGG  >  NZ_CP009273/2849731‑2849901
                                                                                     |                                                                                     
tCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATg                                                                                   <  1:232622/90‑1 (MQ=255)
 cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                                  >  2:27037/1‑90 (MQ=255)
 cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                                  >  1:534144/1‑90 (MQ=255)
             ggCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGa                                                                      <  1:604632/90‑1 (MQ=255)
                           aaaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCgg                                                        >  2:691103/1‑90 (MQ=255)
                           aaaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCgg                                                        <  2:579198/90‑1 (MQ=255)
                            aaTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGc                                                       <  2:296048/90‑1 (MQ=255)
                             aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                                      <  2:478538/90‑1 (MQ=255)
                                 gCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTc                                                  <  1:74319/90‑1 (MQ=255)
                                 gCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTc                                                  >  1:443210/1‑90 (MQ=255)
                                 gCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTc                                                  <  1:254726/90‑1 (MQ=255)
                                              tGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGc                                     >  1:262816/1‑90 (MQ=255)
                                                cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCtt                                   >  1:473136/1‑90 (MQ=255)
                                                          gACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                         <  1:774642/90‑1 (MQ=255)
                                                              gATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagat                     >  1:407664/1‑90 (MQ=255)
                                                               aTCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagatg                    <  2:573906/90‑1 (MQ=255)
                                                                  tGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACg                 <  1:736123/90‑1 (MQ=255)
                                                                   gCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                <  2:676971/90‑1 (MQ=255)
                                                                    ccTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa               <  2:262816/90‑1 (MQ=255)
                                                                    ccTGTCCGGCATACTATCCATGAGGTTTCTTCGGACGATATTTTTCCGGCAGTTTTGGCACCGGACGCTTGTCATCGATGAGATGACGCa               <  1:122135/90‑1 (MQ=255)
                                                                        tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                >  1:656250/1‑85 (MQ=255)
                                                                        tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc                <  2:656250/85‑1 (MQ=255)
                                                                         ccGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGtt          <  2:200897/90‑1 (MQ=255)
                                                                           ggCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTaa        >  2:84463/1‑90 (MQ=255)
                                                                                 cTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTAGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCgg  >  1:42007/1‑90 (MQ=255)
                                                                                     |                                                                                     
TCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGG  >  NZ_CP009273/2849731‑2849901

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGC  >  NZ_CP009273/2849717‑2849914
                                                                                                   |                                                                                                  
GCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATC                                                                                                    <  SRR3722073.127588/100‑1 (MQ=60)
 CACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCC                                                                                                   <  SRR3722073.48652/100‑1 (MQ=60)
 CACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCC                                                                                                   <  SRR3722073.630312/100‑1 (MQ=60)
     ACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGA                                                                                               >  SRR3722073.545289/1‑100 (MQ=60)
      ataagagacagCGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAG                                                                                              <  SRR3722073.13956/89‑1 (MQ=60)
              TCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTC                                                                                      <  SRR3722073.237484/100‑1 (MQ=60)
                           GGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTC                                                                         <  SRR3722073.617218/100‑1 (MQ=60)
                                     TGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTC                                                               >  SRR3722073.452692/1‑100 (MQ=60)
                                               GCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGA                                                     <  SRR3722073.260075/100‑1 (MQ=60)
                                               GCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGA                                                     <  SRR3722073.76043/100‑1 (MQ=60)
                                                  TTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGC                                                  >  SRR3722073.268361/1‑100 (MQ=60)
                                                    GATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTT                                                >  SRR3722073.483203/1‑100 (MQ=60)
                                                                  TTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGAT                                  >  SRR3722073.416372/1‑100 (MQ=60)
                                                                        GACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCA                            <  SRR3722073.790250/100‑1 (MQ=60)
                                                                            GATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGT                        >  SRR3722073.669711/1‑100 (MQ=60)
                                                                                TGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAG                    <  SRR3722073.751016/100‑1 (MQ=60)
                                                                                  CCTGTCCGGCATACTATCCATGAGGTTTCTTCGGACGATATTTTTCCGGCAGTTTTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGAT                  <  SRR3722073.124724/100‑1 (MQ=60)
                                                                                     GTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTAGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGG               >  SRR3722073.43023/1‑100 (MQ=60)
                                                                                                 ATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAG   <  SRR3722073.772346/100‑1 (MQ=60)
                                                                                                  TCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGC  <  SRR3722073.663054/100‑1 (MQ=60)
                                                                                                   |                                                                                                  
GCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGC  >  NZ_CP009273/2849717‑2849914

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