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,882,413 A→C V349G (GTT→GGT)  cysI ← assimilatory sulfite reductase (NADPH) hemoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,882,4130AC100.0% 82.4 / NA 24V349G (GTT→GGT) cysIassimilatory sulfite reductase (NADPH) hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (12/12);  total (12/12)

GGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCG  >  NZ_CP009273/2882328‑2882497
                                                                                     |                                                                                    
ggTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCa                                                                                  >  1:433140/1‑90 (MQ=255)
                       aaTCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCgg                                                           <  1:866680/90‑1 (MQ=255)
                          cAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGa                                                        <  1:617766/90‑1 (MQ=255)
                            aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                      <  1:318703/90‑1 (MQ=255)
                            aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                      <  2:547307/90‑1 (MQ=255)
                            aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                      <  1:547416/90‑1 (MQ=255)
                            aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                      <  1:537155/90‑1 (MQ=255)
                                 gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                      >  2:646257/1‑85 (MQ=255)
                                 gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                      <  1:646257/85‑1 (MQ=255)
                                  cGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCatat                                                >  1:238147/1‑90 (MQ=255)
                                             gATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCgg                                     >  1:66056/1‑90 (MQ=255)
                                                     gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAAtt                             >  1:634311/1‑90 (MQ=255)
                                                     gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAAtt                             >  2:402289/1‑90 (MQ=255)
                                                      cGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAAttt                            >  2:543761/1‑90 (MQ=255)
                                                       gTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTg                           >  1:396650/1‑90 (MQ=255)
                                                             tGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATcccc                     >  1:383007/1‑90 (MQ=255)
                                                                 aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc                 <  2:816428/90‑1 (MQ=255)
                                                                    tATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGAc              <  2:383007/90‑1 (MQ=255)
                                                                            aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGCCGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt      <  2:537198/90‑1 (MQ=255)
                                                                            aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt      >  2:740212/1‑90 (MQ=255)
                                                                               cccTTACCCCATCCAATACGAACGCCCCGTCCGGTTAACTCATATTGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCCCTTcc   <  2:694015/90‑1 (MQ=255)
                                                                                ccTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCg  >  2:655331/1‑90 (MQ=255)
                                                                                 cTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc                                       >  2:532962/1‑52 (MQ=255)
                                                                                 cTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc                                       <  1:532962/52‑1 (MQ=255)
                                                                                     |                                                                                    
GGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCG  >  NZ_CP009273/2882328‑2882497

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 30 ≤ 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
CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG  >  NZ_CP009273/2882316‑2882506
                                                                                                 |                                                                                             
CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                             <  SRR3722072.241131/100‑1 (MQ=60)
CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                             <  SRR3722072.302085/100‑1 (MQ=60)
CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                             <  SRR3722072.881806/100‑1 (MQ=60)
 tataagagacaGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCC                                                                                            <  SRR3722072.6526/89‑1 (MQ=60)
  CCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCA                                                                                           >  SRR3722072.441024/1‑100 (MQ=60)
                                   AATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATA                                                          <  SRR3722072.881625/100‑1 (MQ=60)
                                    ATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATAT                                                         >  SRR3722072.242128/1‑100 (MQ=60)
                                      CAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGG                                                       <  SRR3722072.628785/100‑1 (MQ=60)
                                        AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                     <  SRR3722072.557259/100‑1 (MQ=60)
                                        AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                     <  SRR3722072.657791/100‑1 (MQ=60)
                                        AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                     <  SRR3722072.546825/100‑1 (MQ=60)
                                        AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                     <  SRR3722072.324282/100‑1 (MQ=60)
                                               GACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGG                                              >  SRR3722072.67276/1‑100 (MQ=60)
                                                       tgggctcggagatgtgtataagagacagCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATT                                      <  SRR3722072.542565/72‑1 (MQ=60)
                                                       TCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATT                                      >  SRR3722072.645631/1‑100 (MQ=60)
                                                         GATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTG                                    >  SRR3722072.403808/1‑100 (MQ=60)
                                                               CAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCC                              >  SRR3722072.389853/1‑100 (MQ=60)
                                                                            CAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCtgtctcttatacacatctccgagcccacgagacgtccgattatctcgtatgccgtcttctgc                 >  SRR3722072.624553/1‑38 (MQ=60)
                                                                                     TCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTT        >  SRR3722072.93533/1‑100 (MQ=60)
                                                                                         TGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGGTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAA    >  SRR3722072.528054/1‑100 (MQ=60)
                                                                                          GCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCTACGCTCCACTTCCGCTTTAAAC   >  SRR3722072.599457/1‑100 (MQ=60)
                                                                                          GCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAAC   >  SRR3722072.439705/1‑100 (MQ=60)
                                                                                          GCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAAC   >  SRR3722072.793797/1‑100 (MQ=60)
                                                                                           CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCctgtctcttatacacatctccgagcccacgagacgtccgattatct  >  SRR3722072.745657/1‑54 (MQ=60)
                                                                                           CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG  >  SRR3722072.11339/1‑100 (MQ=60)
                                                                                           CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG  >  SRR3722072.305483/1‑100 (MQ=60)
                                                                                                 |                                                                                             
CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG  >  NZ_CP009273/2882316‑2882506

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