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 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

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% 80.0 / 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)

CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC  >  NZ_CP009273/2882326‑2882498
                                                                                       |                                                                                     
ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     <  1:1048836/90‑1 (MQ=255)
ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     >  1:323154/1‑90 (MQ=255)
ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     >  2:61534/1‑90 (MQ=255)
             gACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGAt                                                                        >  2:1043145/1‑90 (MQ=255)
             gACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGGGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGAt                                                                        >  2:52262/1‑90 (MQ=255)
                 cgcCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGcc                                                                    <  1:487213/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTGCCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  1:1118646/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  1:34913/90‑1 (MQ=255)
                                gATGCGACCATTTTAGATAAACAGCGTCAGGTGCCAGTTATCATCAATTCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTc                                                     >  2:943979/1‑90 (MQ=255)
                                  tGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCat                                                   >  1:510579/1‑90 (MQ=255)
                                   gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                  >  1:1126934/1‑90 (MQ=255)
                                    cGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCatat                                                 >  2:125664/1‑90 (MQ=255)
                                         aTTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  2:916484/79‑1 (MQ=255)
                                         aTTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       >  1:916484/1‑79 (MQ=255)
                                                aTAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                  >  1:316126/1‑77 (MQ=255)
                                                aTAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                  <  2:316126/77‑1 (MQ=255)
                                                    aCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCaa                                 <  1:478372/90‑1 (MQ=255)
                                                                   aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc                  <  2:998199/90‑1 (MQ=255)
                                                                   aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc                  <  2:510579/90‑1 (MQ=255)
                                                                      tATCATCAATGCCCTTCCCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGAc               <  1:515843/90‑1 (MQ=255)
                                                                      tATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGAc               <  2:1126934/90‑1 (MQ=255)
                                                                      tATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGAc               <  1:641337/90‑1 (MQ=255)
                                                                           tcaATGCCCTTCCCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTcc          <  1:747636/90‑1 (MQ=255)
                                                                            caATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCa         >  1:617042/1‑90 (MQ=255)
                                                                              aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt       >  1:925115/1‑90 (MQ=255)
                                                                                   cTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGc  >  2:961623/1‑90 (MQ=255)
                                                                                   cTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGc  >  1:234440/1‑90 (MQ=255)
                                                                                       |                                                                                     
CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC  >  NZ_CP009273/2882326‑2882498

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG  >  NZ_CP009273/2882314‑2882506
                                                                                                   |                                                                                             
ATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAC                                                                                               >  SRR3722089.417483/1‑100 (MQ=60)
 TCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACC                                                                                              <  SRR3722089.963347/100‑1 (MQ=60)
 TCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACC                                                                                              <  SRR3722089.758633/100‑1 (MQ=60)
 TCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACC                                                                                              <  SRR3722089.489275/100‑1 (MQ=60)
  CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                             <  SRR3722089.1256962/100‑1 (MQ=60)
  CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                             <  SRR3722089.429766/100‑1 (MQ=60)
  CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                             >  SRR3722089.327221/1‑100 (MQ=60)
       GCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCC                                                                                        <  SRR3722089.1048529/100‑1 (MQ=60)
       GCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCC                                                                                        <  SRR3722089.248591/100‑1 (MQ=60)
            CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATAC                                                                                   <  SRR3722089.1064554/100‑1 (MQ=60)
               GTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCctgtctct                                                                                >  SRR3722089.97299/1‑92 (MQ=60)
                             CGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTG                                                                  <  SRR3722089.492999/100‑1 (MQ=60)
                                    TAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCAT                                                           >  SRR3722089.516591/1‑100 (MQ=60)
                                     AATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATA                                                          >  SRR3722089.1144044/1‑100 (MQ=60)
                                          AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTGCCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                     <  SRR3722089.1135608/100‑1 (MQ=60)
                                          AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                     <  SRR3722089.35349/100‑1 (MQ=60)
                                           GGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACC                                                    >  SRR3722089.929838/1‑100 (MQ=60)
                                                  ACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCctg                                             >  SRR3722089.320132/1‑97 (MQ=60)
                                                         atggtcgtcggcagcgtcagatgtgtataagagacagCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATT                                      <  SRR3722089.771104/63‑1 (MQ=60)
                                                                ACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCC                               <  SRR3722089.484055/100‑1 (MQ=60)
                                                                              CAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCA                 >  SRR3722089.624773/1‑100 (MQ=60)
                                                                                GTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACT               >  SRR3722089.938613/1‑100 (MQ=60)
                                                                                  TATCATCAATGCCCTTCCCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTC             <  SRR3722089.521898/100‑1 (MQ=60)
                                                                                  TATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTC             <  SRR3722089.649525/100‑1 (MQ=60)
                                                                                     CATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC          >  SRR3722089.237417/1‑100 (MQ=60)
                                                                                       TCAATGCCCTTCCCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTT        <  SRR3722089.757503/100‑1 (MQ=60)
                                                                                             CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACctgtctcttatacacatctgacgctgccaacgaccatcttagtgtagatct  >  SRR3722089.726775/1‑49 (MQ=60)
                                                                                             CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCtgtctcttatacacatctgacgctgccg  >  SRR3722089.436238/1‑72 (MQ=60)
                                                                                             CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG  >  SRR3722089.218908/1‑100 (MQ=60)
                                                                                             CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG  >  SRR3722089.331800/1‑100 (MQ=60)
                                                                                             CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATACGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG  >  SRR3722089.451738/1‑100 (MQ=60)
                                                                                                   |                                                                                             
ATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG  >  NZ_CP009273/2882314‑2882506

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: