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,778,934 (T)6→5 intergenic (‑564/+187) BW25113_RS25705 ← / ← BW25113_RS13825 DUF5507 domain‑containing protein/tRNA‑Ile

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,778,9290T.100.0% 213.6 / NA 49intergenic (‑559/+192)BW25113_RS25705/BW25113_RS13825DUF5507 domain‑containing protein/tRNA‑Ile
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (25/24);  total (25/24)

AGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTC  >  NZ_CP009273/2778841‑2779012
                                                                                        |                                                                                   
aGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAtt                                                                                    >  2:256981/1‑90 (MQ=255)
aGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAtt                                                                                    >  2:76145/1‑90 (MQ=255)
aGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCAtt                                                                                    >  1:343589/1‑90 (MQ=255)
            aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACt                                                                       >  1:510531/1‑90 (MQ=255)
            aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACt                                                                       >  2:314293/1‑90 (MQ=255)
               gATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACtttt                                                                    <  1:314293/90‑1 (MQ=255)
                     aCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGa                                                              >  2:421889/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCAGATGTAAGCAATTTGCCCGCTTGTCTTAGCATGTATTTTTCAAGGTACTCTTGAATGAGtc                                                           >  1:815658/1‑89 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                           >  1:56251/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                           >  2:718581/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                           >  1:1151702/1‑90 (MQ=255)
                        ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                           >  1:1011305/1‑90 (MQ=255)
                          ggCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACt                                                         <  2:397196/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      <  1:1199928/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      <  2:122237/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      <  1:421889/90‑1 (MQ=255)
                             tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      <  1:1073907/90‑1 (MQ=255)
                              ttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATg                                                     >  2:1196740/1‑90 (MQ=255)
                               tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGt                                                            >  1:561193/1‑82 (MQ=255)
                               tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGt                                                            <  2:561193/82‑1 (MQ=255)
                               tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATgg                                                    >  2:362259/1‑90 (MQ=255)
                                aTCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGaa                                                                 >  1:288224/1‑76 (MQ=255)
                                aTCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGaa                                                                 <  2:288224/76‑1 (MQ=255)
                                  cGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa                                                 <  2:345747/90‑1 (MQ=255)
                                    ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTACAGGTACTTTTGAATGAGTACTGATGGATa                                                 >  1:355238/1‑88 (MQ=255)
                                    ttAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTACAGGTACTTTTGAATGAGTACTGATGGATa                                                 <  2:355238/88‑1 (MQ=255)
                                     tAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAAt                                              >  2:1006596/1‑90 (MQ=255)
                                         tACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACAt                                          <  1:668703/90‑1 (MQ=255)
                                             tAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACAt                                          >  2:411580/1‑86 (MQ=255)
                                             tAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACAt                                          <  1:411580/86‑1 (MQ=255)
                                              aaTTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAg                                     <  2:506928/90‑1 (MQ=255)
                                                    acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGc                                       >  2:1068561/1‑82 (MQ=255)
                                                    acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGc                                       <  1:1068561/82‑1 (MQ=255)
                                                    acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGt                               >  2:872824/1‑90 (MQ=255)
                                                    acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGt                               >  2:377808/1‑90 (MQ=255)
                                                         tAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCAc                          <  1:1139319/90‑1 (MQ=255)
                                                         tAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCAc                          <  2:1098902/90‑1 (MQ=255)
                                                         tAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCAc                          <  1:868339/90‑1 (MQ=255)
                                                          aaGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACg                         >  2:26841/1‑90 (MQ=255)
                                                           aGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGt                        <  2:25839/90‑1 (MQ=255)
                                                           aGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGt                        <  2:198756/90‑1 (MQ=255)
                                                              aaTTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGACGGATAAATACATTGCAg                                     >  1:1046806/1‑74 (MQ=255)
                                                              aaTTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGACGGATAAATACATTGCAg                                     <  2:1046806/74‑1 (MQ=255)
                                                               aTTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCa                    >  1:442573/1‑90 (MQ=255)
                                                               aTTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCa                    >  2:251477/1‑90 (MQ=255)
                                                                 ttGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCaaa                  <  1:718581/90‑1 (MQ=255)
                                                                 ttGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCaaa                  >  2:502928/1‑90 (MQ=255)
                                                                 ttGCCCGCTTGGCATAGCAGGCATGTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCaaa                  <  1:343877/90‑1 (MQ=255)
                                                                       gCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCa            >  1:663766/1‑90 (MQ=255)
                                                                       gCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCa            >  1:249880/1‑90 (MQ=255)
                                                                            gCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCt       <  2:86907/90‑1 (MQ=255)
                                                                                 gCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTc  <  2:172391/90‑1 (MQ=255)
                                                                                        |                                                                                   
AGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTC  >  NZ_CP009273/2778841‑2779012

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTCGAA  >  NZ_CP009273/2778843‑2779015
                                                                                           |                                                                                 
GGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACT                                                                          >  SRR3722089.516543/1‑100 (MQ=60)
            TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCAGATGTAAGCAATTTGCCCGCTTGTCTTAGCATGTATTTTTCAAGGTACTCTTGAATGAGTC                                                              >  SRR3722089.827078/1‑100 (MQ=60)
            TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTA                                                              >  SRR3722089.1026369/1‑100 (MQ=60)
            TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTA                                                              >  SRR3722089.1169235/1‑100 (MQ=60)
            TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTA                                                              >  SRR3722089.56946/1‑100 (MQ=60)
             GATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTAC                                                             <  SRR3722089.318276/100‑1 (MQ=60)
                   ACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGG                                                       >  SRR3722089.567745/1‑100 (MQ=60)
                    CTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGActgt                                                      >  SRR3722089.291913/1‑96 (MQ=60)
                        GGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTACAGGTACTTTTGAATGAGTACTGATGGATAAA                                                  >  SRR3722089.359626/1‑100 (MQ=60)
                           TTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                                               <  SRR3722089.1090049/100‑1 (MQ=60)
                           TTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                                               <  SRR3722089.1218171/100‑1 (MQ=60)
                           TTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                                               <  SRR3722089.426978/100‑1 (MQ=60)
                                       TACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCG                                   <  SRR3722089.677355/100‑1 (MQ=60)
                                       TACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCG                                   <  SRR3722089.416557/100‑1 (MQ=60)
                                          TTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGC                                <  SRR3722089.1084601/100‑1 (MQ=60)
                                                  GCATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGACGGATAAATACATTGCAGTGGCGTGCCACtgtct                        >  SRR3722089.1062494/1‑95 (MQ=60)
                                                   CATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCA                       >  SRR3722089.447878/1‑100 (MQ=60)
                                                       TAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAAC                   <  SRR3722089.1156645/100‑1 (MQ=60)
                                                       TAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAAC                   <  SRR3722089.880790/100‑1 (MQ=60)
                                                           CAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCA               >  SRR3722089.672324/1‑100 (MQ=60)
                                                           CAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCA               >  SRR3722089.253069/1‑100 (MQ=60)
                                                               TTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCC           <  SRR3722089.727922/100‑1 (MQ=60)
                                                               TTGCCCGCTTGGCATAGCAGGCATGTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCC           <  SRR3722089.348159/100‑1 (MQ=60)
                                                                        TGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTCGAA  >  SRR3722089.750571/1‑100 (MQ=60)
                                                                                           |                                                                                 
GGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTCGAA  >  NZ_CP009273/2778843‑2779015

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