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 I225 R1 227 21.4 1179702 97.1% 1145490 86.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,678,723 A→C G48G (GGT→GGG glyA ← serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,678,7230AC100.0% 63.3 / NA 21G48G (GGT→GGGglyAserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (11/10);  total (11/10)

AGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTG  >  NZ_CP009273/2678637‑2678806
                                                                                      |                                                                                   
aGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAcccc                                                                                  >  1:52729/1‑90 (MQ=255)
            atacCAACATCCCCCCACCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCAtt                                                                      <  2:529942/90‑1 (MQ=255)
              atCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAc                                                                    >  1:125171/1‑90 (MQ=255)
              atCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAc                                                                    >  1:220855/1‑90 (MQ=255)
                     tACTCGCAATCGCCGTAGTAGCGTTTGCCCGGATAGCCTTCAGCTGGTTTGTTGGTCAGCTTAGACCCCTGCGCCTGCATTACGCGCggg                                                             >  1:468843/1‑90 (MQ=255)
                          gCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTggg                                                        >  2:544772/1‑89 (MQ=255)
                           cacccgccgTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGgtg                                                       <  1:126913/87‑1 (MQ=255)
                              ccgccgTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAg                                                    <  1:296749/90‑1 (MQ=255)
                                           gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAgg                                          >  2:517225/1‑87 (MQ=255)
                                           gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAgg                                          <  1:517225/87‑1 (MQ=255)
                                           gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                       >  1:1396/1‑90 (MQ=255)
                                           gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                       >  1:560924/1‑90 (MQ=255)
                                            tttGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGAt                                      <  1:148509/90‑1 (MQ=255)
                                             ttGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATc                                     >  1:96255/1‑90 (MQ=255)
                                                               gCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTcc                   <  2:387813/90‑1 (MQ=255)
                                                                     ttgttgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGt             >  1:142740/1‑90 (MQ=255)
                                                                         tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACtt         >  1:501061/1‑90 (MQ=255)
                                                                              aGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTcc    <  2:517325/90‑1 (MQ=255)
                                                                              aGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTcc    <  1:12665/90‑1 (MQ=255)
                                                                               gCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTGGTTTTCGGAGGCGGTCAGTTCGATGTGCTCTTCCTGAGGTACTGTTTCCt   >  1:469126/1‑90 (MQ=255)
                                                                               gCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCt   <  1:262289/90‑1 (MQ=255)
                                                                                cTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTg  <  2:122191/90‑1 (MQ=255)
                                                                                      |                                                                                   
AGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTG  >  NZ_CP009273/2678637‑2678806

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 29 ≤ 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
ATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  >  NZ_CP009273/2678627‑2678822
                                                                                                |                                                                                                   
ATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCC                                                                                                  >  SRR3722109.53358/1‑100 (MQ=60)
              GTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAC                                                                                    >  SRR3722109.223495/1‑100 (MQ=60)
              GTTCAACCATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAC                                                                                    >  SRR3722109.126711/1‑100 (MQ=60)
                     GATATCAACATACTCGCAATCGCCGTAGTAGCGTTTGCCCGGATAGCCTTCAGCTGGTTTGTTGGTCAGCTTAGACCCCTGCGCCTGCATTACGCGCGGG                                                                             >  SRR3722109.475902/1‑100 (MQ=60)
                      atacaaacatccccccacCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGC                                                                            <  SRR3722109.229588/82‑1 (MQ=60)
                                     CACCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGG                                                             <  SRR3722109.128476/100‑1 (MQ=60)
                                        CCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGG                                                          <  SRR3722109.300580/100‑1 (MQ=60)
                                           CCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGA                                                       >  SRR3722109.569434/1‑100 (MQ=60)
                                           CCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGA                                                       >  SRR3722109.1412/1‑100 (MQ=60)
                                             GTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATC                                                     >  SRR3722109.97454/1‑100 (MQ=60)
                                                  AGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTC                                                <  SRR3722109.525090/100‑1 (MQ=60)
                                                      TTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATG                                            <  SRR3722109.150341/100‑1 (MQ=60)
                                                                     TTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGT                             >  SRR3722109.144513/1‑100 (MQ=60)
                                                                         GCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTT                         >  SRR3722109.508622/1‑100 (MQ=60)
                                                                               TTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTGGTTTTCGGAGGCGGTCAGTTCGATGTGCTCTTCCTGAGGTACTGTTTCCT                   >  SRR3722109.476189/1‑100 (MQ=60)
                                                                                        AGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAG          <  SRR3722109.12843/100‑1 (MQ=60)
                                                                                         GCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGC         <  SRR3722109.265489/100‑1 (MQ=60)
                                                                                                CCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  >  SRR3722109.543026/1‑100 (MQ=60)
                                                                                                CCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  <  SRR3722109.11860/100‑1 (MQ=60)
                                                                                                |                                                                                                   
ATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  >  NZ_CP009273/2678627‑2678822

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