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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

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% 46.3 / NA 15V349G (GTT→GGT) cysIassimilatory sulfite reductase (NADPH) hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (5/10);  total (5/10)

CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACT  >  NZ_CP009273/2882326‑2882493
                                                                                       |                                                                                
ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                <  2:81301/90‑1 (MQ=255)
ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                >  1:326763/1‑90 (MQ=255)
        cAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAAt                                                                        <  2:570281/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                  <  1:397110/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                  <  2:622109/90‑1 (MQ=255)
                               ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAcccc                                                                               >  2:661301/1‑60 (MQ=255)
                               ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAcccc                                                                               <  1:661301/60‑1 (MQ=255)
                                          ttttCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGCCCGGTAAACTCATACGGActg                                      >  1:222599/1‑88 (MQ=255)
                                             tCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTCCCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc                                   <  2:411583/90‑1 (MQ=255)
                                             tCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc                                   <  2:46149/90‑1 (MQ=255)
                                                  aaaCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTc                              <  1:16213/90‑1 (MQ=255)
                                                    aCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTCCCCCAGCCAATACGATCGCCTCGCCCGGTGAACTCATATGGACGGATCGGTTCaa                            <  1:264663/90‑1 (MQ=255)
                                                    aCAGCGTCAGGTCCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCaa                            <  1:283016/90‑1 (MQ=255)
                                                       gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGTATCGGTTCAAAtt                         >  2:18014/1‑90 (MQ=255)
                                                                 ccAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgc               <  1:719788/90‑1 (MQ=255)
                                                                 ccAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgc               <  2:407598/90‑1 (MQ=255)
                                                                             aaTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCAc   >  1:356088/1‑90 (MQ=255)
                                                                              aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt  >  2:240876/1‑90 (MQ=255)
                                                                              aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt  >  1:200455/1‑90 (MQ=255)
                                                                                       |                                                                                
CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACT  >  NZ_CP009273/2882326‑2882493

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTC  >  NZ_CP009273/2882316‑2882510
                                                                                                 |                                                                                                 
CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                                 <  SRR3722073.273479/100‑1 (MQ=60)
CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                                 >  SRR3722073.333702/1‑100 (MQ=60)
CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                                 <  SRR3722073.358218/100‑1 (MQ=60)
     GCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCC                                                                                            <  SRR3722073.522251/100‑1 (MQ=60)
           cagatgtgtataagagacagGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACG                                                                                      <  SRR3722073.674856/80‑1 (MQ=60)
                                        agatctacacgaccgccatcgtcggcagcgtcagatgtgtataagagacaGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                         <  SRR3722073.727796/50‑1 (MQ=60)
                                        AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                         <  SRR3722073.405563/100‑1 (MQ=60)
                                          GATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGCCCGGTAAACTCATACGGACTG                                                       >  SRR3722073.227215/1‑100 (MQ=60)
                                                            AAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATC                                     <  SRR3722073.16622/100‑1 (MQ=60)
                                                              ACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTCCCCCAGCCAATACGATCGCCTCGCCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCC                                   <  SRR3722073.270249/100‑1 (MQ=60)
                                                              ACAGCGTCAGGTCCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCC                                   <  SRR3722073.289010/100‑1 (MQ=60)
                                                                           CCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCC                      <  SRR3722073.734326/100‑1 (MQ=60)
                                                                             AGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCAC                    >  SRR3722073.363624/1‑100 (MQ=60)
                                                                              GTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACT                   >  SRR3722073.204510/1‑100 (MQ=60)
                                                                                     TCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTT            >  SRR3722073.672663/1‑100 (MQ=60)
                                                                                           CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG      >  SRR3722073.346373/1‑100 (MQ=60)
                                                                                           CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG      >  SRR3722073.772386/1‑100 (MQ=60)
                                                                                           CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG      >  SRR3722073.779109/1‑100 (MQ=60)
                                                                                               TACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTC  <  SRR3722073.245931/100‑1 (MQ=60)
                                                                                                 |                                                                                                 
CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTC  >  NZ_CP009273/2882316‑2882510

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