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 1,980,420 A→G G5G (GGA→GGG ftnB → non‑heme ferritin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,980,4200AG100.0% 142.8 / NA 41G5G (GGA→GGGftnBnon‑heme ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/20);  total (21/20)

CTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACA  >  NZ_CP009273/1980336‑1980506
                                                                                    |                                                                                      
cTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCt                                                                                   <  2:408672/90‑1 (MQ=255)
cTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCt                                                                                   <  1:893931/90‑1 (MQ=255)
             ttATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAaactcaac                                                                      <  1:756009/90‑1 (MQ=255)
               aTGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAaactcaactc                                                                    >  2:653938/1‑90 (MQ=255)
                 ggTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACtctc                                                                  >  1:427915/1‑90 (MQ=255)
                   tgtTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCaa                                                                >  1:955975/1‑90 (MQ=255)
                             tATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACcgcg                                                      <  2:188493/90‑1 (MQ=255)
                                   gCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGtttt                                                >  1:960091/1‑90 (MQ=255)
                                        tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTAc                                              <  2:654715/87‑1 (MQ=255)
                                        tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTAc                                              >  1:654715/1‑87 (MQ=255)
                                        tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                           >  1:239749/1‑90 (MQ=255)
                                        tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                           >  2:869348/1‑90 (MQ=255)
                                        tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                           >  2:546861/1‑90 (MQ=255)
                                        tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                           >  1:952126/1‑90 (MQ=255)
                                           ccACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATcc                                        <  1:325251/90‑1 (MQ=255)
                                             ccTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCaa                                      <  1:474035/89‑1 (MQ=255)
                                             aCTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCaa                                      <  1:529257/90‑1 (MQ=255)
                                                tCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtct                                   <  1:192443/90‑1 (MQ=255)
                                                tCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtct                                   >  2:434132/1‑90 (MQ=255)
                                                 cTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                                  <  1:86505/90‑1 (MQ=255)
                                                 cTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                                  <  1:394303/90‑1 (MQ=255)
                                                 cTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                                  <  2:399003/90‑1 (MQ=255)
                                                   ggATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTa                                >  2:807585/1‑90 (MQ=255)
                                                         aaGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCa                          <  2:21341/90‑1 (MQ=255)
                                                             atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  2:866585/1‑90 (MQ=255)
                                                             atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  1:928055/1‑90 (MQ=255)
                                                             atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  2:731287/1‑90 (MQ=255)
                                                              tatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGa                     >  2:1028950/1‑90 (MQ=255)
                                                               atTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACcgcg                                                      >  2:753038/1‑56 (MQ=255)
                                                               atTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACcgcg                                                      <  1:753038/56‑1 (MQ=255)
                                                                 tAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTa                  <  2:1020683/90‑1 (MQ=255)
                                                                     tATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                                  >  2:360665/1‑70 (MQ=255)
                                                                     tATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                                  <  1:360665/70‑1 (MQ=255)
                                                                      aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             >  2:380548/1‑90 (MQ=255)
                                                                      aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             <  1:380548/90‑1 (MQ=255)
                                                                       tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt            >  1:440698/1‑90 (MQ=255)
                                                                       tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt            >  2:18261/1‑90 (MQ=255)
                                                                       tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt            >  2:988122/1‑90 (MQ=255)
                                                                           aaCCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTc        <  2:928055/90‑1 (MQ=255)
                                                                                 tGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACa  <  1:866585/90‑1 (MQ=255)
                                                                                 tGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACa  <  1:908890/90‑1 (MQ=255)
                                                                                    |                                                                                      
CTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACA  >  NZ_CP009273/1980336‑1980506

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAAC  >  NZ_CP009273/1980322‑1980516
                                                                                                  |                                                                                                
AGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGA                                                                                                 <  SRR3722089.309679/100‑1 (MQ=60)
              CTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTC                                                                                   <  SRR3722089.906864/100‑1 (MQ=60)
                     CAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTC                                                                            >  SRR3722089.433070/1‑100 (MQ=60)
                       GTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAA                                                                          >  SRR3722089.970032/1‑100 (MQ=60)
                           TTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGA                                                                      <  SRR3722089.766081/100‑1 (MQ=60)
                                       GATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTT                                                          >  SRR3722089.974237/1‑100 (MQ=60)
                                           gcgtcagatgtgtataagagacaGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGC                                                      <  SRR3722089.763053/77‑1 (MQ=60)
                                            ATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCA                                                     >  SRR3722089.663132/1‑100 (MQ=60)
                                            ATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCA                                                     >  SRR3722089.966111/1‑100 (MQ=60)
                                            ATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCA                                                     >  SRR3722089.242793/1‑100 (MQ=60)
                                                         CCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACC                                        <  SRR3722089.329338/100‑1 (MQ=60)
                                                           CCTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTT                                      <  SRR3722089.479687/100‑1 (MQ=60)
                                                           ACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTT                                      <  SRR3722089.535430/100‑1 (MQ=60)
                                                              TCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCAC                                   <  SRR3722089.194862/100‑1 (MQ=60)
                                                               taagagacagGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACC                                  <  SRR3722089.365113/90‑1 (MQ=60)
                                                               CTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACC                                  <  SRR3722089.399108/100‑1 (MQ=60)
                                                               CTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACC                                  <  SRR3722089.87572/100‑1 (MQ=60)
                                                                 GGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTG                                >  SRR3722089.941595/1‑100 (MQ=60)
                                                                           ATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGT                      >  SRR3722089.445989/1‑100 (MQ=60)
                                                                                    ATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAAC             <  SRR3722089.385195/100‑1 (MQ=60)
                                                                                           CCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCT      >  SRR3722089.82078/1‑100 (MQ=60)
                                                                                             GCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGA    >  SRR3722089.630844/1‑100 (MQ=60)
                                                                                               TGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAAC  <  SRR3722089.879001/100‑1 (MQ=60)
                                                                                               TGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAAC  <  SRR3722089.922108/100‑1 (MQ=60)
                                                                                                  |                                                                                                
AGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAAC  >  NZ_CP009273/1980322‑1980516

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