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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

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

TTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACA  >  NZ_CP009273/1980334‑1980506
                                                                                      |                                                                                      
ttCTACAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATg                                                                                     <  1:442823/90‑1 (MQ=255)
    gCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTtc                                                                                 >  2:518006/1‑90 (MQ=255)
       aaCAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCa                                                                              <  1:518006/90‑1 (MQ=255)
         cAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCaaa                                                                            >  1:388155/1‑90 (MQ=255)
               ttATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAaactcaac                                                                      <  2:222243/90‑1 (MQ=255)
                 aTGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAaactcaactc                                                                    <  2:833754/90‑1 (MQ=255)
                   ggTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACtctc                                                                  >  2:856677/1‑90 (MQ=255)
                   ggTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACtctc                                                                  <  1:846863/90‑1 (MQ=255)
                              ttATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACcgc                                                       <  2:691574/90‑1 (MQ=255)
                              ttATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACcgc                                                       >  2:269134/1‑90 (MQ=255)
                               tATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACcgcg                                                      <  2:420462/90‑1 (MQ=255)
                                     gCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGtttt                                                >  2:916497/1‑90 (MQ=255)
                                          tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                           >  2:512619/1‑90 (MQ=255)
                                              cACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCa                                       <  1:329868/90‑1 (MQ=255)
                                                  tCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtct                                   >  1:795679/1‑90 (MQ=255)
                                                  tCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtct                                   >  1:762823/1‑90 (MQ=255)
                                                   cTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                                  <  2:488014/90‑1 (MQ=255)
                                                     ggATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTa                                >  1:690612/1‑90 (MQ=255)
                                                          taAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTc                           >  1:761580/1‑90 (MQ=255)
                                                               atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  1:858362/1‑90 (MQ=255)
                                                               atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  2:505894/1‑90 (MQ=255)
                                                               atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  2:360253/1‑90 (MQ=255)
                                                               atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  1:763531/1‑90 (MQ=255)
                                                               atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  2:13387/1‑90 (MQ=255)
                                                                  ttAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGt                   >  2:67518/1‑90 (MQ=255)
                                                                       tATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTg              >  1:481233/1‑90 (MQ=255)
                                                                        aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             >  2:386420/1‑90 (MQ=255)
                                                                        aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             <  2:657664/90‑1 (MQ=255)
                                                                        aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             >  1:62792/1‑90 (MQ=255)
                                                                        aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             >  1:560830/1‑90 (MQ=255)
                                                                               ccGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTg      >  2:284932/1‑90 (MQ=255)
                                                                                cGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGa     >  1:151742/1‑90 (MQ=255)
                                                                                   tGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACa  <  1:897739/90‑1 (MQ=255)
                                                                                      |                                                                                      
TTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACA  >  NZ_CP009273/1980334‑1980506

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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
ATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAAC  >  NZ_CP009273/1980333‑1980516
                                                                                       |                                                                                                
ATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAA                                                                                      >  SRR3722072.395106/1‑100 (MQ=60)
 TTCTACAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAAC                                                                                     <  SRR3722072.450927/100‑1 (MQ=60)
        AACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTC                                                                              <  SRR3722072.527330/100‑1 (MQ=60)
                    GGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCG                                                                  <  SRR3722072.861509/100‑1 (MQ=60)
                                         CCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCT                                             >  SRR3722072.776267/1‑100 (MQ=60)
                                         CCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCT                                             >  SRR3722072.809602/1‑100 (MQ=60)
                                            CTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTA                                          >  SRR3722072.702962/1‑100 (MQ=60)
                                               CACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCT                                       <  SRR3722072.335658/100‑1 (MQ=60)
                                                 CTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTC                                     >  SRR3722072.775005/1‑100 (MQ=60)
                                                      GGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTG                                >  SRR3722072.776987/1‑100 (MQ=60)
                                                      GGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTG                                >  SRR3722072.873179/1‑100 (MQ=60)
                                                              GGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTG                        >  SRR3722072.489962/1‑100 (MQ=60)
                                                               GATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGG                       >  SRR3722072.570893/1‑100 (MQ=60)
                                                               GATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGG                       >  SRR3722072.63946/1‑100 (MQ=60)
                                                                       GTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGA               >  SRR3722072.154413/1‑100 (MQ=60)
                                                                                    TGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAAC  <  SRR3722072.913190/100‑1 (MQ=60)
                                                                                       |                                                                                                
ATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAAC  >  NZ_CP009273/1980333‑1980516

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