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 4,180,586 (A)7→6 intergenic (+333/+77) yjaZ → / ← thiH heat‑shock protein C/2‑iminoacetate synthase ThiH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,180,5800A.100.0% 101.1 / NA 25intergenic (+327/+83)yjaZ/thiHheat‑shock protein C/2‑iminoacetate synthase ThiH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (10/15);  total (10/15)

TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG  >  NZ_CP009273/4180492‑4180660
                                                                                        |                                                                                
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  1:17083/1‑90 (MQ=255)
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  2:253392/1‑90 (MQ=255)
  gCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaaaa                                                                               <  1:426696/90‑1 (MQ=255)
          gAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAAc                                                                      >  1:683616/1‑90 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  1:307849/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:770897/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  1:539359/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:201431/90‑1 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  1:169208/90‑1 (MQ=255)
              ttttAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCtt                                                                  <  1:253392/90‑1 (MQ=255)
                       aGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAcc                                                         >  1:61226/1‑90 (MQ=255)
                        ggCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCt                                                        >  2:226281/1‑90 (MQ=255)
                                 aaaCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGttttt                                               >  2:580613/1‑90 (MQ=255)
                                 aaaCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGttttt                                               >  1:582077/1‑90 (MQ=255)
                                  aaCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    >  1:538336/1‑68 (MQ=255)
                                  aaCGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  2:538336/68‑1 (MQ=255)
                                             agCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAg                                   >  2:610775/1‑90 (MQ=255)
                                                  ttCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTGTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCt                              >  2:300646/1‑90 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                             <  1:858951/90‑1 (MQ=255)
                                                                caTTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc                <  2:292982/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt              <  2:149192/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt              <  2:432836/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt              <  1:887702/90‑1 (MQ=255)
                                                                  ttCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt              <  1:569084/90‑1 (MQ=255)
                                                                     ctgctgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCt           >  2:637139/1‑90 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    >  2:841348/1‑90 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  2:683616/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  1:380121/90‑1 (MQ=255)
                                                                                        |                                                                                
TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG  >  NZ_CP009273/4180492‑4180660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAGGC  >  NZ_CP009273/4180492‑4180680
                                                                                              |                                                                                              
TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAAC                                                                                          >  SRR3722072.695854/1‑100 (MQ=60)
  GCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGC                                                                                        <  SRR3722072.434457/100‑1 (MQ=60)
            ATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAC                                                                              <  SRR3722072.172134/100‑1 (MQ=60)
            ATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAC                                                                              <  SRR3722072.313191/100‑1 (MQ=60)
            ATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAAC                                                                              <  SRR3722072.549066/100‑1 (MQ=60)
             TTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACC                                                                             >  SRR3722072.62345/1‑100 (MQ=60)
              TTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCT                                                                            <  SRR3722072.257689/100‑1 (MQ=60)
                       AGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTT                                                                   >  SRR3722072.592470/1‑100 (MQ=60)
                        GGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTgtctcttata                                                                  >  SRR3722072.548028/1‑90 (MQ=60)
                                                   TCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATT                                       <  SRR3722072.873782/100‑1 (MQ=60)
                                                                 gtgggctcggagatgtgtataagagacagACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAT                          <  SRR3722072.695136/71‑1 (MQ=60)
                                                                  TTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTC                        <  SRR3722072.579251/100‑1 (MQ=60)
                                                                  TTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTC                        <  SRR3722072.902967/100‑1 (MQ=60)
                                                                              GCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC            <  SRR3722072.386912/100‑1 (MQ=60)
                                                                              GCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTC            >  SRR3722072.539822/1‑100 (MQ=60)
                                                                                        TAAAAAACACAACGCTTTCAACAACCTTGCGTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAGGC  >  SRR3722072.7/1‑100 (MQ=60)
                                                                                              |                                                                                              
TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTCATAGTCTTTGCGAGGC  >  NZ_CP009273/4180492‑4180680

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