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 I206 R1 206 23.7 1314362 96.8% 1272302 86.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,995,567 A→G T635A (ACC→GCC)  xdhA → xanthine dehydrogenase molybdenum‑binding subunit XdhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,995,5670AG100.0% 67.7 / NA 22T635A (ACC→GCC) xdhAxanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/11);  total (11/11)

GAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGG  >  NZ_CP009273/2995480‑2995646
                                                                                       |                                                                               
gAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGcc                                                                               >  1:419602/1‑90 (MQ=255)
gAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGcc                                                                               >  1:158107/1‑90 (MQ=255)
  aaGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCAt                                                                             >  2:220902/1‑90 (MQ=255)
         aTCAAAACCACCTCTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGc                                                                      <  2:163270/90‑1 (MQ=255)
                   ccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACg                                                            >  1:208542/1‑90 (MQ=255)
                        aaCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                       <  2:562795/90‑1 (MQ=255)
                              ccGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCa                                                 <  2:69450/90‑1 (MQ=255)
                                   gTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCa                                            <  2:278840/90‑1 (MQ=255)
                                   gTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCa                                            <  2:390845/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                          <  2:643069/90‑1 (MQ=255)
                                           gTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTa                                    >  2:610998/1‑90 (MQ=255)
                                           gTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTa                                    >  2:593936/1‑90 (MQ=255)
                                              ccTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATc                                 <  1:593936/90‑1 (MQ=255)
                                                 ttgttgATCTGACGGTCTGTATTGCGCGGTGCAAGGTCGGTATCAAACGCATCCTTAAGGGTCATGGTTCAGGGCATATTCTTAATCCAc                              >  1:94489/1‑90 (MQ=255)
                                                           gACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGa                    >  2:214024/1‑90 (MQ=255)
                                                           gACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGa                    <  1:90497/90‑1 (MQ=255)
                                                              ggTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAgg                 <  2:208542/90‑1 (MQ=255)
                                                                 cGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCa              <  2:525369/90‑1 (MQ=255)
                                                                 cGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCa              <  1:370716/90‑1 (MQ=255)
                                                                        gcgcTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTacac       >  2:414565/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg  >  1:459320/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg  >  1:205792/1‑90 (MQ=255)
                                                                                       |                                                                               
GAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGG  >  NZ_CP009273/2995480‑2995646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ 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
GCTCTCTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCT  >  NZ_CP009273/2995470‑2995666
                                                                                                 |                                                                                                   
GCTCTCTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCC                                                                                                   >  SRR3722088.159691/1‑100 (MQ=60)
GCTCTCTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCC                                                                                                   >  SRR3722088.424268/1‑100 (MQ=60)
      ctagtacggtctcgtgggctcggagatgtgtataagagacagTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAAC                                                                                             <  SRR3722088.117905/58‑1 (MQ=60)
      TGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAAC                                                                                             <  SRR3722088.318059/100‑1 (MQ=60)
                   ATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACG                                                                                >  SRR3722088.210623/1‑100 (MQ=60)
                                                 GGGTGTGCCTTTGTTGATCTGACGGTCTGTATTGCGCGGTGCAAGGTCGGTATCAAACGCATCCTTAAGGGTCATGGTTCAGGGCATATTCTTAATCCAC                                                  >  SRR3722088.95441/1‑100 (MQ=20)
                                                        CCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGC                                           <  SRR3722088.601242/100‑1 (MQ=60)
                                                                     GACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTA                              <  SRR3722088.91414/100‑1 (MQ=60)
                                                                           CGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgt                        >  SRR3722088.191028/1‑44 (MQ=60)
                                                                           CGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGC                        <  SRR3722088.374648/100‑1 (MQ=60)
                                                                             ATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGG                      >  SRR3722088.207853/1‑100 (MQ=60)
                                                                             ATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGG                      >  SRR3722088.464692/1‑100 (MQ=60)
                                                                                       GTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATG            >  SRR3722088.284805/1‑100 (MQ=60)
                                                                                                 GCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCT  <  SRR3722088.644107/100‑1 (MQ=60)
                                                                                                 |                                                                                                   
GCTCTCTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCT  >  NZ_CP009273/2995470‑2995666

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