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 I229 R1 214 17.4 943020 97.5% 919444 87.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.97.0% 124.6 / ‑1.1 33intergenic (‑559/+192)BW25113_RS25705/BW25113_RS13825DUF5507 domain‑containing protein/tRNA‑Ile
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (18/14);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

AATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTC  >  NZ_CP009273/2778853‑2779012
                                                                            |                                                                                   
aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACt                                                                       >  2:342884/1‑90 (MQ=255)
aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACt                                                                       >  1:188765/1‑90 (MQ=255)
aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACt                                                                       >  2:134954/1‑90 (MQ=255)
aaTGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACt                                                                       >  1:398027/1‑90 (MQ=255)
aaTGATAATACACTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACt                                                                       >  1:218616/1‑90 (MQ=255)
       ataCTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAAt                                                                >  1:383438/1‑90 (MQ=255)
            ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGt                                                                          <  2:372255/75‑1 (MQ=255)
            ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGt                                                                          >  1:372255/1‑75 (MQ=255)
            ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                           >  2:467171/1‑90 (MQ=255)
            ctGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTa                                                           >  2:444940/1‑90 (MQ=255)
              ggCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACt                                                         >  1:339677/1‑90 (MQ=255)
                 tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      <  2:328744/90‑1 (MQ=255)
                 tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      <  1:399513/90‑1 (MQ=255)
                 tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      <  1:391437/90‑1 (MQ=255)
                 tttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      <  2:446904/90‑1 (MQ=255)
                  ttATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      <  1:257705/90‑1 (MQ=255)
                   tATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATgg                                                    >  1:191428/1‑90 (MQ=255)
                      cGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATa                                                 <  2:191428/90‑1 (MQ=255)
                         tAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAAt                                              <  1:165921/90‑1 (MQ=255)
                           aTTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAc                                            <  1:107211/90‑1 (MQ=255)
                             tACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACAt                                          <  2:383438/90‑1 (MQ=255)
                                     tCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTgg                                  <  2:218616/90‑1 (MQ=255)
                                        acaTGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGt                               <  1:467171/90‑1 (MQ=255)
                                             tAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATg                                                               <  2:84468/53‑1 (MQ=255)
                                             tAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATg                                                               >  1:84468/1‑53 (MQ=255)
                                              aaGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACg                         >  2:467364/1‑90 (MQ=255)
                                                  aaTTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      <  2:32823/57‑1 (MQ=255)
                                                  aaTTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGAt                                                      >  1:32823/1‑57 (MQ=255)
                                                   aTTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCa                    >  2:359103/1‑90 (MQ=255)
                                                     ttGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCaaa                  <  1:359103/90‑1 (MQ=255)
                                                           gCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCa            >  2:303330/1‑90 (MQ=255)
                                                           gCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCa            >  1:141208/1‑90 (MQ=255)
                                                                     gCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTc  >  1:170866/1‑90 (MQ=255)
                                                                            |                                                                                   
AATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTC  >  NZ_CP009273/2778853‑2779012

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AATGAAGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTC  >  NZ_CP009273/2778836‑2779012
                                                                                                  |                                                                              
AATGAAGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCC                                                                               <  SRR3722113.81728/100‑1 (MQ=60)
       GGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACT                                                                       >  SRR3722113.190245/1‑100 (MQ=60)
       GGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACT                                                                       >  SRR3722113.402194/1‑100 (MQ=60)
       GGGTATAATAAATGATAATACACTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACT                                                                       >  SRR3722113.220365/1‑100 (MQ=60)
              ATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAAT                                                                >  SRR3722113.387446/1‑100 (MQ=60)
                   TGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATctgtc                                                           >  SRR3722113.376115/1‑95 (MQ=60)
                     ATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACT                                                         >  SRR3722113.343148/1‑100 (MQ=60)
                          ACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGG                                                    >  SRR3722113.192923/1‑100 (MQ=60)
                                  TTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                                            <  SRR3722113.403693/100‑1 (MQ=60)
                                  TTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                                            <  SRR3722113.395540/100‑1 (MQ=60)
                                   TTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATAC                                            <  SRR3722113.259964/100‑1 (MQ=60)
                                          TAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGT                                    <  SRR3722113.167271/100‑1 (MQ=60)
                                            ATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGG                                  <  SRR3722113.108161/100‑1 (MQ=60)
                                                    ATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGctgtctcttatacacatctccgagccc                          >  SRR3722113.85209/1‑73 (MQ=60)
                                                         ACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACctgtctcttatacacatctccga                     >  SRR3722113.33099/1‑77 (MQ=60)
                                                         ACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACC                     <  SRR3722113.472120/100‑1 (MQ=60)
                                                                  CAATTTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCA            >  SRR3722113.142379/1‑100 (MQ=60)
                                                                      TTGCCCGCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCC        <  SRR3722113.362786/100‑1 (MQ=60)
                                                                            GCTTGGCATAGCAGGCATTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTC  >  SRR3722113.172241/1‑100 (MQ=60)
                                                                                                  |                                                                              
AATGAAGGGGTATAATAAATGATAATACTCTGGCTTTATCGTTAATTACTTAATTCCACATGTAAGCAATTTGCCCGCTTGGCATAGCAGGCATTTTTTCCAGGTACTTTTGAATGAGTACTGATGGATAAATACATTGCAGTGGCGTGCCACGTACCAAAACACCAGCCCTCATTC  >  NZ_CP009273/2778836‑2779012

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