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 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% 95.5 / NA 27G5G (GGA→GGGftnBnon‑heme ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/7);  total (20/7)

TGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGT  >  NZ_CP009273/1980355‑1980496
                                                                 |                                                                            
tgtTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCaa                                                      >  2:490950/1‑90 (MQ=255)
  ttttGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAAt                                                    >  2:370575/1‑90 (MQ=255)
          tATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACcgcg                                            <  1:186829/90‑1 (MQ=255)
           aTCTTGCATCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGa                                           >  2:538286/1‑90 (MQ=255)
           aTCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGa                                           >  1:722211/1‑90 (MQ=255)
                     tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                 >  1:320424/1‑90 (MQ=255)
                        ccACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATcc                              <  1:296369/90‑1 (MQ=255)
                          aCTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCaa                            >  2:230400/1‑90 (MQ=255)
                             tCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtct                         >  2:278946/1‑90 (MQ=255)
                              cTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                        <  2:722720/90‑1 (MQ=255)
                              cTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                        <  1:230400/90‑1 (MQ=255)
                                ggATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTa                      >  1:554041/1‑90 (MQ=255)
                                ggATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTa                      >  1:84838/1‑90 (MQ=255)
                                  atatAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTAcc                    >  2:179576/1‑90 (MQ=255)
                                   tataAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCCCTACCt                   <  2:413734/90‑1 (MQ=255)
                                     taAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTc                 >  2:649637/1‑90 (MQ=255)
                                     taAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTc                 >  1:202100/1‑90 (MQ=255)
                                          atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg            >  1:552124/1‑90 (MQ=255)
                                              tAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTa        <  1:597878/90‑1 (MQ=255)
                                                  tATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTg    >  1:725142/1‑90 (MQ=255)
                                                  tATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTg    >  2:663819/1‑90 (MQ=255)
                                                  tATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTg    >  1:279952/1‑90 (MQ=255)
                                                   aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg   >  2:386736/1‑90 (MQ=255)
                                                    tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt  >  2:46902/1‑90 (MQ=255)
                                                    tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt  >  2:51058/1‑90 (MQ=255)
                                                          ccGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                        >  1:83877/1‑62 (MQ=255)
                                                          ccGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                        <  2:83877/62‑1 (MQ=255)
                                                                 |                                                                            
TGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGT  >  NZ_CP009273/1980355‑1980496

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
AAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCT  >  NZ_CP009273/1980325‑1980512
                                                                                               |                                                                                            
cagatgtgtataagagacagGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGC                                                                                          <  SRR3722073.205112/80‑1 (MQ=60)
                               GTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGA                                                           >  SRR3722073.736805/1‑100 (MQ=60)
                                        TATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGC                                                  <  SRR3722073.190497/100‑1 (MQ=60)
                                         ATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCA                                                 >  SRR3722073.327214/1‑100 (MQ=60)
                                                    CTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTA                                      >  SRR3722073.565621/1‑100 (MQ=60)
                                                    CTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTA                                      >  SRR3722073.86746/1‑100 (MQ=60)
                                                      CCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACC                                    <  SRR3722073.302641/100‑1 (MQ=60)
                                                         CTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTC                                 >  SRR3722073.206195/1‑100 (MQ=60)
                                                            CTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACC                              <  SRR3722073.235211/100‑1 (MQ=60)
                                                              GGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTG                            >  SRR3722073.563661/1‑100 (MQ=60)
                                                                      GGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTG                    >  SRR3722073.285875/1‑100 (MQ=60)
                                                                      GGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTG                    >  SRR3722073.739778/1‑100 (MQ=60)
                                                                            TAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTC              <  SRR3722073.610347/100‑1 (MQ=60)
                                                                              GGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCctgtctcttatacacatc            >  SRR3722073.85765/1‑82 (MQ=60)
                                                                                        CCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCT  >  SRR3722073.406358/1‑100 (MQ=60)
                                                                                               |                                                                                            
AAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCT  >  NZ_CP009273/1980325‑1980512

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: