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 I206 R1 206 23.7 1314362 96.8% 1272302 86.2

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

GCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACA  >  NZ_CP009273/1980338‑1980506
                                                                                  |                                                                                      
gCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTtc                                                                                 >  2:218593/1‑90 (MQ=255)
      aGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAaac                                                                           >  2:549091/1‑90 (MQ=255)
        tAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAaactc                                                                         <  2:160236/90‑1 (MQ=255)
               ggTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAaactcaac                                                                      >  1:245618/1‑86 (MQ=255)
               ggTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAaactcaac                                                                      <  2:245618/86‑1 (MQ=255)
                  gtTTTGATTTATCTTGCACCTCACCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCaaa                                                               >  1:169616/1‑90 (MQ=255)
                            aTCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGa                                                     >  1:252820/1‑90 (MQ=255)
                            aTCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGa                                                     >  1:655241/1‑90 (MQ=255)
                                 gCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGtttt                                                >  1:368916/1‑90 (MQ=255)
                                      tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                           >  1:97082/1‑90 (MQ=255)
                                              tCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtct                                   <  1:145885/90‑1 (MQ=255)
                                               cTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                                  <  2:579379/90‑1 (MQ=255)
                                                    tataAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCt                             <  1:204328/90‑1 (MQ=255)
                                                           atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  1:65536/1‑90 (MQ=255)
                                                           atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  1:46061/1‑90 (MQ=255)
                                                            tatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGa                     >  1:374871/1‑90 (MQ=255)
                                                            tatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGa                     >  2:485481/1‑90 (MQ=255)
                                                                   tATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCt                             >  1:440828/1‑75 (MQ=255)
                                                                   tATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCt                             <  2:440828/75‑1 (MQ=255)
                                                                    aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             <  1:549091/90‑1 (MQ=255)
                                                                    aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             >  1:655836/1‑90 (MQ=255)
                                                                    aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             >  1:185450/1‑90 (MQ=255)
                                                                     tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt            >  1:301001/1‑90 (MQ=255)
                                                                     tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt            >  2:326236/1‑90 (MQ=255)
                                                                          aCCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCt       >  2:208986/1‑90 (MQ=255)
                                                                               tGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACa  <  2:46061/90‑1 (MQ=255)
                                                                                  |                                                                                      
GCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACA  >  NZ_CP009273/1980338‑1980506

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ 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
CAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGC  >  NZ_CP009273/1980343‑1980519
                                                                             |                                                                                                   
CAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTC                                                                               >  SRR3722088.248026/1‑100 (MQ=60)
   TAATTATGGTGTTTTGATTTATCTTGCACCTCACCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAA                                                                            >  SRR3722088.171323/1‑100 (MQ=60)
             GTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGA                                                                  >  SRR3722088.663359/1‑100 (MQ=60)
             GTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGA                                                                  >  SRR3722088.255284/1‑100 (MQ=60)
                  GATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTT                                                             >  SRR3722088.372837/1‑100 (MQ=60)
                       ATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCA                                                        >  SRR3722088.98077/1‑100 (MQ=60)
                                         TCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCAC                                      <  SRR3722088.147375/100‑1 (MQ=60)
                                            GGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTG                                   >  SRR3722088.46534/1‑100 (MQ=60)
                                            GGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTG                                   >  SRR3722088.66219/1‑100 (MQ=60)
                                             GATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGA                                  >  SRR3722088.378860/1‑100 (MQ=60)
                                               TATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGT                                <  SRR3722088.206375/100‑1 (MQ=60)
                                                    GGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTctgtc                           >  SRR3722088.445904/1‑95 (MQ=60)
                                                     GATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGG                          >  SRR3722088.187309/1‑100 (MQ=60)
                                                     GATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGG                          >  SRR3722088.663961/1‑100 (MQ=60)
                                                      ATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGT                         >  SRR3722088.303957/1‑100 (MQ=60)
                                                               ATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAAC                <  SRR3722088.555802/100‑1 (MQ=60)
                                                                             GATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGC  <  SRR3722088.589742/100‑1 (MQ=60)
                                                                             |                                                                                                   
CAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACGGC  >  NZ_CP009273/1980343‑1980519

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