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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,396,099 G→C A413G (GCC→GGC)  nuoC ← NADH‑quinone oxidoreductase subunit C/D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,396,0990GC100.0% 83.5 / NA 24A413G (GCC→GGC) nuoCNADH‑quinone oxidoreductase subunit C/D
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (13/11);  total (13/11)

GCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTC  >  NZ_CP009273/2396013‑2396184
                                                                                      |                                                                                     
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                    >  1:780956/1‑90 (MQ=255)
gcgcCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAgccgg                                                                                    >  2:190305/1‑90 (MQ=255)
      tACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGc                                                                              >  2:460703/1‑90 (MQ=255)
            cGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTggg                                                                        >  1:747958/1‑90 (MQ=255)
               cGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAAc                                                                     <  1:543607/90‑1 (MQ=255)
               cGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAAc                                                                     <  2:279232/90‑1 (MQ=255)
                 aaGTCGTGCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGa                                                                   <  2:552498/90‑1 (MQ=255)
                  aGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCCCTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGAc                                                                  <  2:123112/90‑1 (MQ=255)
                      gATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCttt                                                              >  1:179734/1‑90 (MQ=255)
                      gATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCttt                                                              >  1:7570/1‑90 (MQ=255)
                            ggTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAAt                                                        <  1:432885/90‑1 (MQ=255)
                                 cACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAAt                                                        <  1:429835/85‑1 (MQ=255)
                                 cACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAAt                                                        >  2:429835/1‑85 (MQ=255)
                                  aCGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGtt                                                  >  2:174636/1‑90 (MQ=255)
                                        gCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAg                                            <  1:891686/90‑1 (MQ=255)
                                             cgcCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGcgccg                                       >  1:797233/1‑90 (MQ=255)
                                                           cACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGAc                         >  1:48803/1‑90 (MQ=255)
                                                            aCTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACg                        >  1:857397/1‑90 (MQ=255)
                                                                    gcCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGt                >  2:221500/1‑90 (MQ=255)
                                                                    gcCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGt                >  1:316008/1‑90 (MQ=255)
                                                                       tctTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTc             >  1:101263/1‑90 (MQ=255)
                                                                       tctTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTc             <  1:638622/90‑1 (MQ=255)
                                                                            cgcgCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAATACGCCAGACGTTTCGGCAt        <  2:179734/90‑1 (MQ=255)
                                                                            cgcgCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCAt        <  2:780956/90‑1 (MQ=255)
                                                                                  aTAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTc  <  2:10301/90‑1 (MQ=255)
                                                                                  aTAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTc  >  2:849820/1‑90 (MQ=255)
                                                                                      |                                                                                     
GCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTC  >  NZ_CP009273/2396013‑2396184

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ 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
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTCGAGGAACTCACGC  >  NZ_CP009273/2396003‑2396197
                                                                                                |                                                                                                  
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGG                                                                                                 >  SRR3722072.794653/1‑100 (MQ=60)
            GCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGG                                                                                     >  SRR3722072.761183/1‑100 (MQ=60)
                      CGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTT                                                                           >  SRR3722072.182846/1‑100 (MQ=60)
                      CGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTT                                                                           >  SRR3722072.7721/1‑100 (MQ=60)
                         CGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAG                                                                        <  SRR3722072.553391/100‑1 (MQ=60)
                                      GGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGC                                                           <  SRR3722072.437660/100‑1 (MQ=60)
                                      GGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGC                                                           <  SRR3722072.440763/100‑1 (MQ=60)
                                             CGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCG                                                    >  SRR3722072.811178/1‑100 (MQ=60)
                                                  GCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTC                                               <  SRR3722072.907037/100‑1 (MQ=60)
                                                           AGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGAC                                      >  SRR3722072.49692/1‑100 (MQ=60)
                                                            GTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACG                                     >  SRR3722072.872201/1‑100 (MQ=60)
                                                                    CCACTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGT                             >  SRR3722072.321534/1‑100 (MQ=60)
                                                                       CTCCAGCGCCTCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTC                          >  SRR3722072.103133/1‑100 (MQ=60)
                                                                                 TCTTTCGCGCCATAGCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGT                <  SRR3722072.650014/100‑1 (MQ=60)
                                                                                               GCCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTCGAGGAACTCACGC  <  SRR3722072.350856/100‑1 (MQ=60)
                                                                                                |                                                                                                  
GAATAAGGACGCGCCTTACGCACGTCGAAGTCGATCCCGGTAGCACGCAGGCCCGCGCCAGTGGTGCCCCACTCCAGCGCCTCTTTCGCGCCATAGGCGGCAACGCCCTGGGAACGACCTTTCAGAATGGTGTTTTGCAGCGCCGCTTTCTCGTAAGACGCCAGACGTTTCGGCATCCAGTCGAGGAACTCACGC  >  NZ_CP009273/2396003‑2396197

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