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 I228 R1 214 26.5 1448192 97.1% 1406194 86.4

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% 87.8 / NA 27G48G (GGT→GGGglyAserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (20/7);  total (20/7)

TTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCC  >  NZ_CP009273/2678639‑2678804
                                                                                    |                                                                                 
tttttCAACAATATCAACATCCCCGCACCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTg                                                                              <  2:653787/87‑1 (MQ=255)
            atCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAc                                                                  >  1:492734/1‑90 (MQ=255)
                  aTACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCgg                                                            >  1:708982/1‑90 (MQ=255)
                    aCTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGc                                                          >  1:112647/1‑90 (MQ=255)
                         cAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGgtg                                                     >  2:59132/1‑90 (MQ=255)
                            ccgccgTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAg                                                  <  2:425724/90‑1 (MQ=255)
                              gccgTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGtt                                                >  2:19199/1‑90 (MQ=255)
                               ccgTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGttt                                               >  2:400729/1‑90 (MQ=255)
                                      agCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAgg                                        >  1:44615/1‑90 (MQ=255)
                                      agCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAgg                                        <  1:692223/90‑1 (MQ=255)
                                      agCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAgg                                        <  1:694832/90‑1 (MQ=255)
                                         gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                     >  2:96335/1‑90 (MQ=255)
                                         gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                     >  2:648790/1‑90 (MQ=255)
                                         gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                     >  2:635682/1‑90 (MQ=255)
                                         gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                     >  2:598753/1‑90 (MQ=255)
                                         gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                     >  2:405338/1‑90 (MQ=255)
                                           ttGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATc                                   >  1:304642/1‑90 (MQ=255)
                                           ttGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATc                                   >  1:222130/1‑90 (MQ=255)
                                                 ggATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCg                             >  1:249123/1‑90 (MQ=255)
                                                       ccTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGc                       <  1:598753/90‑1 (MQ=255)
                                                            aGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTc                  >  2:110141/1‑90 (MQ=255)
                                                                  tttGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACg            >  2:223049/1‑90 (MQ=255)
                                                                   ttgttgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGt           >  1:345436/1‑90 (MQ=255)
                                                                       tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACtt       >  1:163680/1‑90 (MQ=255)
                                                                           cAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTc   >  1:578229/1‑90 (MQ=255)
                                                                            aGCTGAGACCCCTGCGCCTGCATTTCGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTcc  <  1:375216/90‑1 (MQ=255)
                                                                            aGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTcc  <  1:208705/90‑1 (MQ=255)
                                                                                    |                                                                                 
TTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCC  >  NZ_CP009273/2678639‑2678804

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
GTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  >  NZ_CP009273/2678641‑2678822
                                                                                  |                                                                                                   
GTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAC                                                                                    >  SRR3722112.499854/1‑100 (MQ=60)
      CGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGG                                                                              >  SRR3722112.720067/1‑100 (MQ=60)
        ATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGC                                                                            >  SRR3722112.114044/1‑100 (MQ=60)
              caacacccgcatCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGT                                                                      <  SRR3722112.387166/88‑1 (MQ=60)
              acctcccctcacCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGT                                                                      <  SRR3722112.420990/88‑1 (MQ=60)
                          CCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGG                                                          >  SRR3722112.45156/1‑100 (MQ=60)
                               GTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATC                                                     >  SRR3722112.308141/1‑100 (MQ=60)
                               GTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATC                                                     >  SRR3722112.224805/1‑100 (MQ=60)
                                    AGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTC                                                <  SRR3722112.703034/100‑1 (MQ=60)
                                    AGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTC                                                <  SRR3722112.705683/100‑1 (MQ=60)
                                     GCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCG                                               >  SRR3722112.252071/1‑100 (MQ=60)
                                                     CCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGAC                               <  SRR3722112.607926/100‑1 (MQ=60)
                                                       TTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGT                             >  SRR3722112.349749/1‑100 (MQ=60)
                                                           GCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTT                         >  SRR3722112.165697/1‑100 (MQ=60)
                                                               ATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTC                     >  SRR3722112.587056/1‑100 (MQ=60)
                                                                      GGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCC              >  SRR3722112.46794/1‑100 (MQ=60)
                                                                          AGCTGAGACCCCTGCGCCTGCATTTCGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAG          <  SRR3722112.380070/100‑1 (MQ=60)
                                                                          AGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAG          <  SRR3722112.211232/100‑1 (MQ=60)
                                                                                  CCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  >  SRR3722112.422183/1‑100 (MQ=60)
                                                                                  CCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  <  SRR3722112.211226/100‑1 (MQ=60)
                                                                                  |                                                                                                   
GTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  >  NZ_CP009273/2678641‑2678822

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: