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
A122 F1 I1 R1 27 151.2 10056804 92.4% 9292486 76.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,609,130 A→C K285Q (AAG→CAG)  insF‑3 → putative transposase for IS3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,609,1300AC100.0% 194.2 / NA 47K285Q (AAG→CAG) insF‑3putative transposase for IS3
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/21);  total (26/21)

GCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGTA  >  REL606/1609104‑1609202
                          |                                                                        
gcCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCa                         <  2:4791757/76‑1 (MQ=17)
gcCTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCa                         <  1:800425/76‑1 (MQ=17)
  cTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCAcc                       >  2:3420690/1‑76 (MQ=25)
  cTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCAcc                       >  2:2902477/1‑76 (MQ=25)
  cTCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCAcc                       >  2:1392828/1‑76 (MQ=25)
   tCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCt                      >  1:1311420/1‑76 (MQ=255)
   tCAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCt                      <  1:4375739/76‑1 (MQ=255)
    cAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTa                     <  1:183965/76‑1 (MQ=255)
    cAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTa                     <  2:3523161/76‑1 (MQ=255)
    cAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTa                     <  1:1432658/76‑1 (MQ=255)
    cAGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTa                     >  2:2093543/1‑76 (MQ=255)
     aGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaa                    >  2:1390524/1‑76 (MQ=255)
     aGTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaa                    >  1:4877134/1‑76 (MQ=255)
     aGTCCGGAACAATATGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaa                    >  2:896584/1‑76 (MQ=255)
      gTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaaa                   >  1:2660266/1‑76 (MQ=255)
      gTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaaa                   <  1:2138526/76‑1 (MQ=255)
      gTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaaa                   >  2:3719069/1‑76 (MQ=255)
      gTCCGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTaaa                   >  1:4928458/1‑76 (MQ=255)
        ccGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATa                 <  1:806855/76‑1 (MQ=255)
        ccGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATa                 <  1:1074996/76‑1 (MQ=255)
         cGGAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATAc                >  1:2866188/1‑76 (MQ=255)
           gAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGa              >  2:3228154/1‑76 (MQ=255)
           gAACAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGa              <  2:4961039/76‑1 (MQ=255)
            aaCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAg             <  1:389909/76‑1 (MQ=255)
            aaCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAg             <  2:2428078/76‑1 (MQ=255)
            aaCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAg             >  1:4303021/1‑76 (MQ=255)
             aCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAgg            <  1:4881451/76‑1 (MQ=255)
             aCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAgg            <  1:3492342/76‑1 (MQ=255)
             aCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAgg            >  2:4316468/1‑76 (MQ=255)
             aCAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAgg            <  1:1658793/76‑1 (MQ=255)
              cAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGc           <  1:2782381/76‑1 (MQ=255)
              cAATTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGc           >  2:2635261/1‑76 (MQ=255)
                aTTTGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCAt         <  1:1483732/76‑1 (MQ=255)
                 tttGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCAtt        <  2:3044173/76‑1 (MQ=255)
                 tttGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCAtt        >  1:4628460/1‑76 (MQ=255)
                  ttGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTg       >  2:770591/1‑76 (MQ=255)
                   tGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGc      >  2:2394065/1‑76 (MQ=255)
                   tGAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGc      >  1:143872/1‑76 (MQ=255)
                    gAAAACCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCt     >  1:1797303/1‑76 (MQ=255)
                     aaaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTg    >  2:1293314/1‑76 (MQ=255)
                     aaaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTg    >  1:3300429/1‑76 (MQ=255)
                      aaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGt   <  1:3061479/76‑1 (MQ=255)
                      aaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGt   >  2:4083010/1‑76 (MQ=255)
                      aaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGt   <  1:1658227/76‑1 (MQ=255)
                      aaaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGt   >  1:5016428/1‑76 (MQ=255)
                       aaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGTa  <  2:3275178/76‑1 (MQ=255)
                       aaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGTa  <  2:3396815/76‑1 (MQ=255)
                       aaCCAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGTa  >  2:2701067/1‑76 (MQ=255)
                          |                                                                        
GCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAATTCACCTAAATACGAGGCATTGCTGTA  >  REL606/1609104‑1609202

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

GATK/CNVnator alignment

N/A