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 F5 I102 R1 278 39.1 2251254 95.8% 2156701 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 974,173 T→C S798P (TCC→CCC)  mukB → chromosome partition protein MukB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273974,1730TC100.0% 112.1 / NA 31S798P (TCC→CCC) mukBchromosome partition protein MukB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/20);  total (11/20)

TATTCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTTCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGCGTTCAGCCGCTTTATCGGC  >  NZ_CP009273/974086‑974250
                                                                                       |                                                                             
tATTCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCg                                                                            >  1:1085047/1‑91 (MQ=255)
   tCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAAc                                                                         <  1:842230/91‑1 (MQ=255)
   tCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAAc                                                                         <  1:353912/91‑1 (MQ=255)
      cGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCt                                                                      >  2:1022136/1‑91 (MQ=255)
                 aGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGctc                                                           <  2:108947/91‑1 (MQ=255)
                         tGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGa                                                   <  1:1081238/91‑1 (MQ=255)
                         tGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGa                                                   <  1:955813/91‑1 (MQ=255)
                           tttGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATg                                                 <  1:761009/91‑1 (MQ=255)
                            ttGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGt                                                >  1:277173/1‑91 (MQ=255)
                             tGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTa                                               <  2:366341/91‑1 (MQ=255)
                             tGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCAAACGCTTCGCCACGCTCTCCTTTGATGTa                                               <  1:13841/91‑1 (MQ=255)
                               gTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACa                                             <  1:613274/91‑1 (MQ=255)
                                  gTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCttt                                                     <  2:937367/80‑1 (MQ=255)
                                  gTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCttt                                                     >  1:937367/1‑80 (MQ=255)
                                  gTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGaa                                          >  1:347487/1‑91 (MQ=255)
                                  gTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGaa                                          >  2:568614/1‑91 (MQ=255)
                                  gTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGaa                                          >  2:184309/1‑91 (MQ=255)
                                  gTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGaa                                          <  1:156383/91‑1 (MQ=255)
                                        cgcgTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCa                                    <  2:904112/91‑1 (MQ=255)
                                         gcgTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAg                                   <  2:241079/91‑1 (MQ=255)
                                         gcgTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAg                                   >  2:402225/1‑91 (MQ=255)
                                             gAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTc                               <  1:471491/91‑1 (MQ=255)
                                                      aTTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAgg                      <  2:266861/91‑1 (MQ=255)
                                                           aaGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGCGTTc                 >  2:892103/1‑91 (MQ=255)
                                                              ccTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGctc                                                           <  2:1025995/46‑1 (MQ=255)
                                                              ccTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGctc                                                           >  1:1025995/1‑46 (MQ=255)
                                                                tCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCa                                                                <  2:1109260/39‑1 (MQ=39)
                                                                tCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCa                                                                >  1:1109260/1‑39 (MQ=39)
                                                                tCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGCGTTCAGCCg            <  2:709934/91‑1 (MQ=255)
                                                                       cGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGCGTTCAGCCGCTTTATc     <  2:1107995/91‑1 (MQ=255)
                                                                          gCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGCGTTCAGCCGCTTTATCGGc  <  1:1093734/91‑1 (MQ=255)
                                                                                       |                                                                             
TATTCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTTCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGCGTTCAGCCGCTTTATCGGC  >  NZ_CP009273/974086‑974250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ 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
TCAGTGGCGTTATTCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTTCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGCGTTCAGCCGCTTTATCGGCAGTCATCTGGCGGTTGCGTTTGA  >  NZ_CP009273/974076‑974273
                                                                                                 |                                                                                                    
TCAGTGGCGTTATTCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCG                                                                                                   >  SRR3722125.1109472/1‑101 (MQ=60)
    TGGCGTTATTCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACG                                                                                               <  SRR3722125.791351/101‑1 (MQ=60)
             TCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCAC                                                                                      <  SRR3722125.860766/101‑1 (MQ=60)
             TCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCAC                                                                                      <  SRR3722125.361311/101‑1 (MQ=60)
                            GTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGT                                                                       >  SRR3722125.282782/1‑101 (MQ=60)
                                  CTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAC                                                                 >  SRR3722125.958174/1‑101 (MQ=60)
                                  CTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAA                                                                 >  SRR3722125.354703/1‑101 (MQ=60)
                                   TGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAA                                                                <  SRR3722125.977099/101‑1 (MQ=60)
                                   TGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAA                                                                <  SRR3722125.1105560/101‑1 (MQ=60)
                                     TTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAAC                                                              <  SRR3722125.777491/101‑1 (MQ=60)
                                       TGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCAAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTC                                                            <  SRR3722125.14107/101‑1 (MQ=60)
                                         GTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAG                                                          <  SRR3722125.626502/101‑1 (MQ=60)
                                            GTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGT                                                       <  SRR3722125.159371/101‑1 (MQ=60)
                                                       GAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGC                                            <  SRR3722125.481432/101‑1 (MQ=60)
                                                              GTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGctgtctcttatacacatctccgagcccacgagact                                     >  SRR3722125.1048993/1‑66 (MQ=60)
                                                                ATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTCCCGAGCGCTTCGCCACGCTCTCCTTctgtctcttatacacatctccgagcccacgagactaggcatg                                   >  SRR3722125.1134254/1‑59 (MQ=60)
                                                                                    GCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGCGTTCAGCCGCTTTATCGGCAGTCATCTGG               >  SRR3722125.67812/1‑101 (MQ=60)
                                                                                    GCGTGAAGTGCTTCCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGCGTTCAGCCGCTTTATCGGCAGTCATCTGG               <  SRR3722125.1118346/101‑1 (MQ=60)
                                                                                                 CCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCctgtctcttatacacatctccgagcccacgagactaggcatgatctcgtatgccgtct  >  SRR3722125.845046/1‑43 (MQ=60)
                                                                                                 |                                                                                                    
TCAGTGGCGTTATTCACGTTTCCCGGAAGTGCCGCTGTTTGGTCGTGCTGCGCGTGAAAGCCGTATTGAAAGCCTCCATGCCGAGCGTGAAGTGCTTTCCGAACGCTTCGCCACGCTCTCCTTTGATGTACAGAAAACTCAGCGTCTGCATCAGGCGTTCAGCCGCTTTATCGGCAGTCATCTGGCGGTTGCGTTTGA  >  NZ_CP009273/974076‑974273

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